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1.  INTRODUCTION 


1.1  Overview 

This  report  summarizes  the  research  conducted  in  the  second  phase  of  this 
three-year  research  and  development  program  directed  toward  the  analysis 
and  evaluation  of  myoelectric  signals  (MES)  as  indicators  of  operator 
alertness  and  piloting  workload.  The  purpose  of  the  study  was  to 
investigate  the  efficiency  of  stochastic  models  such  as  autoregressive 
(AR),  autoregressive-moving-average  (ARMA),  and  autoregressive  integrated 
moving  average  (ARIMA)  models  in  characterizing  the  MES  under  different 
levels  of  task-imposed  burden. 

The  three  year  program  built  on  the  research  performed  by  Madni  (1978, 
1981),  Graupe  and  Cline  (1975)  and  Graupe,  Magnussen  and  Beex  (1977)  who 
established  the  feasibility  of  stochastic  models  in  characterizing  sampled 
myoelectric  signal  waveforms.  In  particular,  the  work  of  Madni  has 
established  the  feasibility  of  characterizing  myoelectric  signals  under 
varying  levels  of  muscle  tension  and  physical  fatigue  via  stochastic 
models.  This  program  has  explored  methods  for  characterizing  the 
myoelectric  signal  under  varying  mental  stress  levels  of  the  human 
operator. 

The  research  conducted  provided  a  firm  analytical  and  experimental  basis 
for  model-based  feature  extraction  by  providing  techniques  for  algorithmic 
improvement  in  accuracy,  by  investigating  critical  issues  in  pattern 
variations  due  to  individual  operator's  behavioral  differences,  and  by 
providing  techniques  for  symbolic  interpretation.  The  program  culminated 
in  the  development  of  a  prototype  MES  analysis  system  and  guidelines  for 
its  application  in  assessing  operator  stress  and/or  alertness  in  various 
piloting  activities. 
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Problem  Statement 


The  definition  and  derivation  of  objectives  measures  for  assessing 
workload,  attentional  demands,  or  operator  alertness  in  specific  piloting 
tasks  has  been  an  area  of  investigation  by  several  researchers  for  more 
than  three  decades.  Myoelectric  signals  (MES)  have  been  the  subject  of  a 
search  by  some  researchers  (Kennedy  and  Travis,  1947;  Travis  and  Kennedy, 
1947;  Kennedy,  1953)  for  a  physiological  indication  of  alertness  in 
piloting  tasks.  The  results  of  these  experiments  demonstrated  that  there 
appeared  to  be  some  correlation  between  MES  properties  (e.g.,  spike 
amplitude,  zero  crossings)  and  human  alertness.  However,  the  use  of  these 
properties  as  an  indicator  of  alertness  level  was  never  successfully 
incorporated  in  a  practical  setting  primarily  because  of  the  excessively 
high  false  positives  in  certain  tasks.  That  is,  diminished  alertness  has 
been  identified  in  many  situations  when  the  subject  was  clearly  alert.  One 
plausible  explanation  for  this  unreliability  in  "answers"  extracted  from 
MES  signatures  is  that  the  information  content  of  the  original  MES  waveform 
is  underutilized.  In  other  words,  the  reliability  of  features,  the 
information  content  of  the  features,  and  the  feature  extraction  process  are 
critical  to  the  success  of  the  alertness/workload  level  discrimination 
process. 


1.3  Program  Objectives 

Several  of  the  program  objectives  were  met  in  the  first  phase  of  the 
program.  These  included: 


(1)  The  development  and  implementation  of  stochastic  model -based 
signal  processing  and  pattern  analysis  approaches  within  the 
overall  framework  of  the  data  acquisition  and  processing 
system. 
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(2)  The  derivation  of  model  structure,  feature  extraction,  and 
parameter  identification  processes  for  constructing  MES-based 
indicators  of  alertness/workload. 

(3)  The  development  of  an  experimental  plan  and  a  representative 
task  simulation  and  interface. 

(4)  A  pilot  experimental  investigation. 

(5)  The  evaluation  of  MES  features  in  terms  of  their  relevance  to 
operator  alertness  level  and/or  mental  load. 

The  second  phase  of  the  program  built  on  the  work  resulting  from  the  first 
phase  by  enhancing  model  parameter  estimates  and  algorithm  accuracy, 
investigating  experimentally  the  performance  of  the  model -based  feature 
extraction  methods,  determining  pattern  responses  to  different  task 
characteristics  and  individual  differences,  and  establishing  a  rule-based 
framework  for  interpretation  of  the  feature  parameters.  Program  objectives 
met  during  the  second  phase  of  the  program  included: 

(1)  The  expansion  of  the  ARIMA-based  feature  extraction  process  to 
provide  better  estimates  of  the  model  parameters  and  to 
improve  computational  accuracy  of  the  algorithms. 

(2)  Improved  data  collection  system  to  provide  speedy  global  and 
automatic  pick-up  of  MES  signals. 

(3)  Expansion  of  the  task  simulation  to  include  critical  air  pilot 
situations  where  the  alertness  of  the  aircrew  to  high-stress 
and  multiple-task  conditions  is  essential. 
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(4)  Identification  of  range  of  feature  variations  and  model 
accuracy/validity. 

(5)  Development  and  experimental  use  of  the  extended  model/system 
in  extracting  indications  of  operator  alertness/load  level. 

(6)  Development  of  guidelines  for  operational  application  of  the 
MES  analysis  system  in  air  piloting  tasks.  The  guidelines 
should  include:  methods  of  measurements;  model  structuring; 
parameter  estimation  and  interpretation  techniques;  and 
specification  of  the  system's  control  and  display 
requirements. 

1.4  Technical  Approach 


Stochastic  modeling  and  time  series  analysis  methods  have  been  extensively 
used  to  statistically  model  the  relationship  between  the  amplitude  of  the 
signal  at  different  points  in  time  along  the  entire  time  history.  In  this 
model,  the  amplitude  fluctuations  along  the  timeline  are  treated  as  a 
stochastic  process.  Stochastic  models  are  particularly  well  suited  as  a 
temporal  feature  extraction  tool  for  time  varying  random  signals.  Such 
features,  because  of  their  high  information  content,  have  proven  to  be 
diagnostic  indicators  in  applications  where  purely  spectral  or  ad  hoc 
feature  extraction  methods  have  failed  (Madni,  1978).  The  key  hypotheses 
underlying  the  use  of  stochastic  models  as  a  feature  extraction  method  for 
alertness  level  identification  were  that:  (1)  at  least  one  of  the  features 
will  be  relatively  constant  and  repeatable  for  the  mental  load  category  and 
task  during  which  the  signal  was  recorded;  and  (2)  at  least  one  of  the 
nearly  constant  features  for  each  alertness  category/load  condition  will  be 
distinctly  different  for  each  level  thus  enabling  identification  of  the 
category. 


The  stochastic  modeling  approach  for  characterizing  bioelectric  signals  is 
most  applicable  to  model  physiological  data  which  possess  some  or  all  of 
the  following  characteristics: 

(1)  The  data  trace  is  noisy,  i.e.,  data  points  show  random 
fluctuations  in  amplitude  and  are  thus  amenable  to  being 
modelled  as  a  random  sequence. 

(2)  The  classification  problem  is  restricted  to  a  finite, 
previously  established,  number  of  categories. 

(3)  Simpler  features  extraction  methods  such  as  power  spectrum 
analysis,  root-rnean-square-value  estimation,  and  amplitude 
level  coding  fail  to  provide  good  separation  among  the 
classes. 

Other  examples  of  physiological  signals  amenable  to  stochastic  modeling 
include:  (1)  Steady  State  Electro-encephalograms  (EEGs);  (2)  Visual  Evoked 
Response;  and  (3)  Electro-oculograms  (EOGs). 

1.5  Findings 

The  research  performed  during  this  three-year  program  revealed  several 
facts  consistent  with  past  research  in  this  field.  First,  surface 
myoelectric  signals  tend  to  be  much  too  "noisy"  to  extract  consistently 
stable  (i.e.,  invariant)  features  from  the  waveform.  Consequently, 
reliability  of  the  MES  signatures  was  found  to  be  inadequate  for  reliable 
discrimination.  Past  research  in  this  area  where  relatively  stable 
features  were  acquired  were  based  on  the  use  of  intramuscular  Basmajian 
type  electrodes  (Madni,  1978).  With  surface  electrodes,  trends  were 
observed  in  approximately  twenty  percent  of  the  subjects  tested.  However, 
the  trends  were  not  robust  enough  to  make  conclusive  assessments.  With 
respect  to  diagnosticity,  the  MES  pattern  associated  with  the  different 


levels  of  task  difficulty  could  not  be  easily  separated.  Since  great  care 
was  exercised  in  electrode  placement  and  contact  and  in  task  design,  we 
attribute  these  results  to  two  key  factors: 

(1)  The  change  of  electrodes  from  Phase  I  of  the  program  to  Phase 
II  was  motivated  by  shape  and  size  consideration.  However  it 
is  conceivable  that  the  resulting  MES  patterns  were 
contaminated  by  extraneous  signals  even  though  the  signals 
looked  clean  and  devoid  of  60  hertz  interference.  The  main 
reason  for  this  conclusion  stems  from  the  significantly  better 
results  in  terms  of  stable  features  achieved  during  Phase  I  of 
this  program.  This  argument  can  account  in  part  for  the  lack 
of  reliability  in  the  MES  signature  derived  via  ARIMA  model 
analysis. 

(2)  The  selected  muscle  sites,  i.e.,  the  frontalis  and  the 
trapezius  may  not  be  as  correlated  with  operator  alertness/ 
workload  levels  as  originally  envisioned.  This  argument 
accounts  in  a  large  part  for  the  relatively  low  diagnoscity  of 
MES  signatures  in  distinguishing  reliably  among  the  different 
task  loading  conditions. 

The  implications  from  this  three-year  research  program  are  two-fold. 
Surface  myoelectric  activity  is  not  a  reliable  measure  of  operator 
alertness.  During  Phase  I,  the  first  autoregressive  coeffocient  of  the 
ARIMA  model  revealed  a  significant  correlation  with  task  difficulty  level. 
During  Phase  III,  the  *  weights  did  not  show  the  same  trend.  Intramuscular 
electrodes,  on  the  other  hand,  that  do  pick  up  more  reliable  signatures 
have  obvious  drawbacks.  Post  hoc  analysis  of  the  experimental  data 
revealed  that  the  total  number  of  experimental  subjects  which  were 
constrained  by  program  scope  and  size  were  inadequate  in  terms  of  producing 
a  statistically  significant  difference  in  perceived  stress  between  the 
single  and  dual-task  groups.  A  detailed  description  of  findings  is 
presented  in  Section  6  of  this  report. 
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2.  BACKGROUND 


2.1  Rationale  and  Hypothesis 

The  MES,  within  the  context  of  human  performance  and  workload,  has  been 
studied  by  various  researchers  over  the  last  three  decades.  Within  the 
context  of  human  performance,  the  MES  can  be  potentially  used  to  provide  a 
measure  of  either  activity  of  the  muscles  or  the  tension  of  the  muscles. 
When  workload  estimation  is  involved,  data  processing  of  some  kind  has  to 
be  performed  on  the  raw  MES  data.  This  processing  can  range  from 
conventional  signal  processing  and  filtering  methods  to  temporal  feature 
extraction  and  pattern  analysis  methods. 

A  number  of  studies  have  been  carried  out  to  demonstrate  the  practical 
value  of  MES  as  a  measure  of  task  workload  and  performance  quality.  Among 
the  earliest  research  is  that  of  Kennedy  and  Travis  (1947,  1948,  1949)  who 
found  that  the  level  of  the  integrated  MES  recorded  over  the  supraorbital 
facial  area  was  closely  related  to  vigilance  and  tracking  performance. 
Lucaccini  (1968)  observed  similar  changes  in  the  integrated  forearm  flexor 
muscle  MES  during  simple  and  complex  visual  tasks.  He  also  reported  that 
the  average  intrasubject  correlations  between  MES  and  performance  were 
significant  in  both  tasks  (r  =  .21  and  .30  in  simple  and  complex  tasks, 
respectively).  Stern  (1966)  found  that  integrated  neck  MES  rose  initially 
and  fell  thereafter  during  easier  and  more  difficult  (lower  signal 
frequency)  versions  of  a  simple  visual  task. 

It  seems  from  these  results  that  integrated  MES  voltage  is  one  of  the 
better  predictors  of  vigilance  performance,  but  it  has  not  been  universally 
accepted  that  MES  varies  directly  with  vigilance  task  performance.  Eason, 
Beardshall  and  Jaffee  (1965)  interpreted  their  results  as  indicating  that 
sympathetic  activity  decreases  along  with  Central  nervous  system  arousal 


and  vigilance,  but  that  somatic  activity  increases  as  part  of  a 
compensatory  process.  Groll's  (1966)  conclusions  were  essentially  the 
same.  Yet,  judging  from  their  results  and  the  contradictory  findings  of 
others,  muscle  tension,  like  performance,  may  reflect  both  the  processes 
underlying  declines  in  vigilance  performance  and  those  acting  to  counteract 
it. 

For  these  reasons,  stochastic  modeling  approaches  which  examine  the  entire 
temporal  signature,  i.e.,  the  whole  latency  domain,  are  preferable  to  ad 
hoc  isolated-feature  dependent  methods.  It  is  expected  that  the  use  of 
stochastic  model-derived  features  in  conjunction  with  single-features  such 
as  P300  component  amplitude  and  latency  (e.g.,  Israel  et  al ,  1980) 
potentially  offer  a  diagnostic  and  reliable  direction  for  the  assessment  of 
internal  operator  states. 

Jex  and  Allen  (1970)  found  that  rectified  and  suitable  filtered  MES 
recorded  from  the  forearm  of  subjects  showed  a  decrease  in  amplitude  when 
subjects  changed  from  a  resting  to  a  tracking  state.  These  researchers 
also  found  that  grip  pressure  increased  with  an  increase  in  tracking 
difficulty.  Sun,  Keane  and  Stackhouse  (1976)  and  Stackhouse  (1976)  found 
that  MES  from  the  forehead  and  the  forearm  were  correlated  with  task 
loading  in  a  variety  of  aircrew  tasks.  Madni  (1978)  found  a  deterministic 
correlation  between  MES  reorded  during  isometric  contraction  of  the  deltoid 
muscle  at  various  load  levels  and  the  associated  MES  stochastic  model 
parameters. 

From  the  foregoing,  it  appears  that  MES  can  potentially  provide  reliable 
correlates  of  operator  states  such  as  alertness  levels  or  load;  however,  in 
order  to  achieve  this  goal,  computer-based  signal  analysis  and  feature 
extraction  methods  must  be  systematically  applied  to  the  MES  recordings. 
Stochastic  modeling  approaches  which  examine  the  entire  temporal  signature, 
i.e.,  the  whole  latency  domain,  are  preferable  to  ad  hoc  isolated  feature- 
dependent  methods.  It  is  expected  that  the  use  of  stochastic  model-derived 


features  in  conjuction  with  single-features  such  as  P300  component 
amplitude  and  latency  (e.g.,  Israel,  Wickens,  Chesney  and  Donchin,  1980) 
potentially  offer  a  diagnostic  and  reliable  direction  for  the  assessment  of 
internal  operator  states. 

Even  with  stochastic  modeling  approaches,  the  transient  nature  of  ERPs 
require  that  a  fairly  general  form  of  a  stochastic  model  be  used,  i.e.,  one 
that  can  at  least  accommodate  nonstationarity  in  the  mean  of  the  waveform. 
Since  ARIMA  stochastic  models  are  designed  to  handle  nonstationary  means, 
they  offer  a  particularly  promising  framework  for  extracting  informative 
features  from  ERPs. 

2.2  Stochastic  Models 

Stochastic  modeling  or  time  series  analysis  (Box  and  Jenkins,  1970)  have 
been  extensively  used  to  model  the  statistical  relationship  between  the 
amplitude  of  a  signal  at  any  point  in  time  and  the  preceeding  amplitudes 
along  the  time  history.  The  amplitude  fluctuations  along  the  time  line  are 
treated  as  a  stochastic  process.  The  future  course  of  the  process  is 
presumed  t^  be  predictable  from  information  about  its  past. 

Before  describing  these  models,  the  notation  employed  will  be  summarized. 


•^k-l^k^k+1' 


be  a  discrete  time  series  where  X^  is  the  random  variable  X  at 
time  i.  We  denote  the  series  by  [ X 1 . 


Let  u  be  the  mean  of  [ X } ,  called  the  level  of  process. 


Let  [X]  denote  the  series  of  deviations  about  ;  that  is. 


Let  [W]  be  a  series  of  outputs  from  a  white  noise  source  with  a 
mean  zero  and  variance  j  . 

Let  B  be  the  "backward"  shift  operator  for  the  deviation  series 
such  that 


Bxk  =  xk_! 

Hence,  Bmxk  =  xk_m 

Let  7  be  the  backward  difference  operator  for  the  deviation 
series  such  that 

7xk=xk-xk-l=(l-B)xk 

Hence,  7mxk=(l-B)mxk 

The  dependence  of  the  current  value  xk  on  the  past  values  of  x  and 
w  can  be  expressed  in  different  ways  giving  rise  to  several 
different  models. 


(a)  Autoregressive  (AR)  Models.  In  this  model  the  current  value 
of  x  depends  on  the  previous  p  and  x  and  on  the  current  noise 
term  w.  Thus. 


xk=alxk-l  +  a2xk-2  + 


♦  apxk-p  +  wk 


xk  =  £  aixk-2  +  wk 

i  =  l 


The  series  [X]  as  defined  above  is  known  as  the  autoregressive 
process  of  order  p.  The  name  "autoregressive"  arises  from  the 
model's  similarity  to  regression  analysis  and  the  fact  that 
the  variable  x  in  an  AR  model  is  regressed  on  previous  values 
of  itself. 

(b)  Moving  Average  (MA)  Model.  In  the  eguation  for  the  AR  model, 
xk-l  can  be  eliminated  from  the  expression  for  x|<  by 
substituting 


xk-l  *  alxk-2  +  a2xk-3  +  ••••+  apxk-p-l  +  wk-l 


expression  for  xk.  This  formulation  explicitly  treats  the 
series  as  being  observations  on  linearly  filtered  Gaussian 
noise.  A  MA  process  of  order  q  is  given  by 


P 

xk  =  I  &iwk-i  +  *k 
i=l 


(c)  Mixed  Model:  Autoregressive-Moving  Average  (ARMA)  Model.  To 
achieve  flexibility  in  the  fitting  of  actual  time  series,  this 
model  includes  both  the  AR  and  the  MA  terms.  A  (p,q)  ARMA 
model  has  the  form: 


p 

P 

*k  =  l  a1 Xk_ i  +  Wk  - 

1  biwk-i 

II 

i=l 

In  all  three  models  described  above  the  process  of  generating  the 
series  is  assumed  to  be  in  equilibrium  about  a  constant  mean 
level.  Models  characterized  by  such  an  equilibrium  condition 
are  called  stationary  models.  In  certain  time  series  data, 
the  level  u  does  not  remain  constant,  i.e.,  the  series  is 
nonstationary.  The  series  may,  nevertheless,  exhibit 
homogeneous  or  stationary  behavior  after  the  differences  due 
to  level  drift  have  been  accounted  for.  It  can  be  shown  that 
such  behavior  can  in  certain  instances  be  represented  by  an 
autoregressive-integrated-moving-average  (ARIMA)  model. 


fl  D 


(d)  Autoregressive-Inteqrated-Movinq-Averaqe  (ARIMA)  Model.  The 
general  (p,d,q)  model  has  the  form 


p 

P 

-  l  ai  vdxk_-j  +  wk  -  l  b-jwk_i 

i=l 

II 

where  xk  is  the  original  time  series 

v  is  the  backward  difference  operator 
d  is  the  number  of  differencing  operations  performed 
on  the  original  data 

p  is  the  order  of  the  autoregressive  terms 
q  is  the  order  of  the  moving  average  terms 


P 

P 

Then  Yk  =  I  aiyk_i+wk  - 

I  biwk-i 

i=l 

i=l 

This  model  is  referred  to  as  a  general  (p,d,q)  model  referring 
to  a  general  pth  order  autoregressive,  dth  data  differencing, 
qth  order  moving  average  process  (Box  et  al,  1970). 


2-7 


ARIMA  Models  in  MES  Characterization 


The  feasibility  of  ARIMA  Stochastic  Model  Identification  for  feature 
extraction  was  explored  by  Madni  (1978).  The  key  elements  of  this  study 
are  provided  in  the  following  paragraphs. 

The  experimental  data  consisted  of  MES  records  from  the  deltoid  muscle  for 
different  isometric  contraction  levels.  These  ranged  from  0%  to  100% 
(maximum),  where  100%  tension  is  defined  as  100%  of  the  force  generated  at 
maximum  effort,  not  100%  of  MES.  The  primary  assumption  in  this  experiment 
is  that  an  X%  run  corresponds  to  X%  of  muscle  tension  which  is  proportional 
to  abduction,  and  that  the  only  muscle  involved  in  abduction  is  the 
deltoid. 

The  results  of  the  spectral  analysis  performed  on  the  experimental  data 
revealed  a  gradual  but  definite  shift  of  power  to  lower  frequencies  with  an 
increase  in  muscle  contraction.  The  total  power  of  the  signal  was  found  to 
lie  below  2500  Hz.  The  most  significant  shift  of  power  to  lower 
frequencies  with  increasing  muscle  tension  was  observed  in  the  frequency 
band  that  contained  ninety  percent  of  the  total  power. 

ARIMA  models  were  fitted  to  the  MES  data  recorded  for  each  contraction 
level.  The  ARIMA  parameters  were  fitted  across  the  n  trials  for  each 
contraction  level.  It  was  shown  that  AR  terms  of  the  a  vector  do  not 

change  significantly  for  the  1%,  5%, . ,50%  tension  levels.  However, 

the  AR  coefficients  for  the  100%  tension  level  is  quite  different  from 
those  for  all  other  tension  levels  (both  in  sign  and  magnitude).  The 
contraction  level  ranges  (within  which  the  AR  coefficients  of  the  ARIMA 
pattern  vector  are  relatively  stationary)  that  resulted  from  this 
experiment  were: 


(1)  1%  through  50%  (or  'low')  contraction  range 

(2)  50%  through  100%  (or  'high')  contraction  range 

Each  of  these  two  contraction  ranges  can  be  represented  by  an  average 
pattern  vector.  A  question  that  arises  here  is  how  one  can  determine, 
online,  the  underlying  contraction  levels  from  the  MES  spectral  signature. 
The  results  from  the  spectral  signature  analysis  reveal  that  90%  of  the 
cumulative  power  was  below  400  Hz  for  the  high  tension  case,  but  was  above 
400  Hz  for  the  low  contraction  level  case.  This  fact  provided  a  useful 
criteria  for  determining  whether  a  given  MES  pattern  should  be  compared  to 
the  'low'  contraction  pattern  vector  or  the  'high'  contraction  pattern 
vector.  Details  of  the  online  recognition  rule  may  be  found  in  Madni 
(1978). 
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3.  SYSTEM  IMPLEMENTATION  AND  EXPERIMENTAL  SET-UP 


3.1  Overview 

The  system  architecture  was  designed  and  implemented  in  parallel  with  the 
first  series  of  experiments  for  assessing  the  "goodness"  of  model-derived 
features  in  terms  of  their  relevance  to  operator  alertness/workload  levels. 
The  primary  tool  used  for  performing  this  assessment  in  both  series  of 
experiments  is  a  computer-based  stochastic  signal  processing  and  pattern 
recognition  algorithm  described  in  earlier  sections. 

This  section  discusses  the  behavioral  issues  being  investigated  through  the 
use  of  model-based  feature  derivation/extraction  of  workload  correlates. 
Workload  correlates  were  primarily  gathered  through  a  task  simulation  based 
on  the  Criterion  Task  Set  (CTS)  workload  test  battery  developed  at  AFAMRL 
(Shingledecker,  1983).  The  task  simulation  was  altered  to  some  degree  to 
accommodate  the  experimental  hardware  constraints  and  to  control  the 
experimental  and  behavioral  variables  during  the  task.  Each  subject  was 
presented  with  controlled  workload  tasks  along  cognitive  processing  and 
motor  task  dimensions.  As  the  subject  performed  the  various  tasks,  the  MES 
data  was  recorded  at  rest,  at  the  beginning,  and  near  the  end  of  task 
execution.  At  the  conclusion  of  each  task,  model  outputs,  subject 
performance,  and  subjective  rating  comparisons  were  made  between  levels  of 
task  loading  and  among  the  various  task  dimensions. 

3.2  ARIMA  Model  Identification  of  MES 

The  key  research  problems  forming  the  basis  of  this  study  and  underlying 
the  use  of  stochastic  models  as  a  feature  extraction  method  are  described 
in  the  following  paragraphs. 
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The  fundamental  hypothesis  forming  the  basis  of  this  study  is  that  MES 
recorded  from  selected  muscle  groups  are  correlated  with  internal  states  of 
the  human  operator  (e.g.,  alertness  level  or  mental  load)  and  that  a 
suitably  selected/designed  feature  extraction  method  is  capable  of 
uncovering  the  underlying  operator  state  in  terms  of  invariant  features 
associated  with  the  MES.  It  was  our  hypothesis  that  stochastic  model 
characterization  of  the  MES  waveforms  is  potentially  capable  of  "capturing" 
features  that  are  both  repeatable  and  diagnostical.  Repeatability  implies 
that  there  is  at  least  one  parameter  in  the  stochastic  model 
characterization  of  MES  data  that  is  constant  or  near-constant  for  each 
underlying  level  of  alertness  or  load  in  a  given  task.  Diagnosticity 
implies  that  these  nearly  invariant  features  are  sufficiently  different  for 
each  level  of  alertness,  thereby  allowing  identification  of  the  underlying 
operator  state.  The  specific  stochastic  model  selected  for  MES 
characterization  is  the  ARIMA  model. 

As  shown  earlier,  an  ARIMA  model  for  a  general  time  series  has  d  levels  of 
differencing,  p  autoregressi ve  coefficients,  and  q  moving  average 
coefficients  as  shown  in  the  equation  below: 


;dZt-<i)i!7dZt_i  -  - 't)p7dZt_p=at-eiat-l  -  - eqat-q 


where 


dZu  is  the  dth  difference  of  the  time  series  at  time  u 

au  is  the  zero  mean,  normally-distributed  random  noise  at  time  u 

u 

°i  is  the  iu  autoregressive  coefficient 

0-j  is  the  ith  moving  average  coefficient 


Determination  of  p,  d,  and  q  is  a  three  step  procedure.  The  first  step  of 
the  ARIMA  modeling  as  provided  by  Box  &  Jenkins  (1970)  is  to  identify  p,d, 
and  q.  The  specific  software  module  associated  with  this  step  calculates 


autocorrelations  and  partial  autocorrelations  for  different  levels  of  data 
differencing.  These  autocorrelations  and  partial  autocorrelations  provide 
an  insight  in  selecting  p,d,  and  q  of  the  ARIMA  model.  To  determine  d,  one 
looks  at  the  autocorrelations  for  a  given  level  of  differencing  and 
observes  whether  or  not  they  "die  out"  rapidly.  If  they  do,  the  given 
level  of  differencing  is  adequate;  if  not,  additional  differencing 
operations  are  required  until  this  constraint  is  satisfied.  The  smallest 
level  of  differencing  for  which  the  autocorrelations  die  out  rapidly  is 
taken  as  the  optimal  level  of  differencing,  d. 

To  determine  p  and  q  we  look  at  the  autocorrelations  and  partial  auto¬ 
correlations  for  the  selected  level  of  differencing.  Box  and  Jenkins 
summarize  this  approach  as  follows: 

"Briefly,  whereas  the  autocorrelation  function  of  an 
autoregressive  process  of  order  p  tails  off,  its 
partial  autocorrelation  function  has  a  cutoff  after 
lag  p.  Conversely,  the  autocorrelation  function  of  a 
moving  average  process  of  order  q  has  a  cutoff  after 
lag  q,  while  its  partial  autocorrelation  tails  off. 

If  both  the  autocorrelations  and  partial  autocor¬ 
relations  tail  off,  a  mixed  process  is  suggested. 

Furthermore,  the  autocorrelation  function  for  a  mixed 

X-  L. 

process,  containing  a  p1,  order  autoregressive 
component  and  a  q^h  order  moving  average  component,  is 
a  mixture  of  exponentials  and  damped  sine  waves  after 
the  first  p-q  lags.  Conversely,  the  partial  auto¬ 
correlation  function  for  a  mixed  process  is  dominated 
by  a  mixture  of  exponentials  and  damped  sine  waves 
after  the  first  p-q  lags." 


Once  p,  d,  and  q  are  selected,  we  proceed  to  the  second  stage  of  the  model 
identification  process.  The  purpose  of  the  second  stage  is  to  come  up  with 
initial  estimates  of  the  autoregressi ve  and  moving  average  parameters. 
These  initial  estimates  are  then  used  by  the  third  stage  of  the  model  in 
generating  final  estimates  of  the  autoregressive  parameters.  The  result  of 
this  three-stage  process  is  a  feature  vector  consisting  of  parameters  that 
parsimoniously  characterize  the  original  MES  time  series  data. 

3.3  AR1MA  Software  (AUTOBJ)  and  Functional  Description 

The  ARIMA  modeling  of  MES  was  based  on  a  recently  released  program,  AUTOBJ 
(an  automatic  Box  &  Jenkins  modeling  procedure)  developed  by  Automatic 
Forecasting  Systems,  Inc.  for  use  on  an  IBM  PC/XT  with  an  8087 
floating-point  chip.  The  previous  system,  the  COSMOS/UNIX  system,  did  not 
support  a  floating-point  hardware.  Consequently,  it  was  necessary  in  phase 
one  of  the  effort  to  develop  a  batch  mode  program  that  allowed  the 
experimenter  to  specify  up  to  25  files  representing  25  experimental 
sessions  and  process  all  the  records  in  these  files.  AUTOBJ,  the  software 
on  the  IBM  PC/XT  with  the  supporting  floating-point  hardware  does  not 
require  this  batch  processing.  The  key  features  of  AUTOBJ  are  described  in 
Table  3-1. 

TABLE  3-1 

KEY  FEATURES  OF  AUTOBJ 


.  Computes  the  optimal  level  of  differencing. 

.  Easily  handles  a  mixed  model  (ARMA)  (the  in-house  system 
generally  set  q  to  zero). 

.  Discards  0  and  e  weights  of  low  significance. 

.  Computes  -  weights  as  opposed  to  AR  parameters  for 
comparison. 


< 
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The  models  obtained  for  the  software  used  during  phase  one  were  generally 
AR  models  whose  coefficients  were  directly  comparable,  while  the  models 
produced  by  AUTOBJ  that  have  differing  values  of  p,  d,  and  q  may  not  be 
directly  compared.  This  problem  of  incompatibility  was  solved  by  using  the 
model  to  generate  *  weights  and  using  the  T  weights  for  comparison  rather 
than  the  model  feature  vector. 

Any  ARIMA  (p,d,q)  model  can  be  expressed  as  an  infinite  series  of  weights 
as  follows: 

Take  the  general  ARIMA  (p,d,q)  model. 


Divide  both  sides  by  (1-h^(B))  to  yield. 


The  resulting  polynomial  coefficients  of  the  potentially  infinite  order 
polynomial  in  B  represent  the  ~  weights.  If  q=0,  then  the  polynomial  is  of 
order  p  and  identical  to  the  6  p(B)  polynomial.  It  is  important  to 
truncate  the  0  weights  after  some  number  such  as  (p+d+q)  to  parsimoniously 
represent  the  time  series. 
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With  its  floating-point  hardware  support  and  automated  identification 
procedure,  the  AUTOBJ  program  was  a  very  powerful  tool  for  the  purposes  of 
this  investigation. 


AUTOBJ  allows  efficient  implementation  of  Box  &  Jenkins  models  and  produce 
time  series  identifications.  Starting  with  raw  data,  the  program's 
algorithm  follows  the  model  building  methodology  described  by  Box  and 
Jenkins  —  that  of  tentative  model  identification,  parameter  estimation  and 
diagnostic  checking  (for  the  purposes  of  this  project,  only  the  model 
identification  was  employed).  While  requiring  only  minimal  input  by  the 
user,  the  program  is  quite  flexible  in  allowing  for  optimal  user  control. 
AUTOBJ  provides  a  quick,  convenient  and  powerful  means  for  developing  ARIMA 
models  for  the  MES  time  series  data. 

To  be  more  explicit,  AUTOBJ  consists  of  two  stages.  In  the  first  stage  the 
system  queries  the  user  for  information  concerning  the  data.  Stage  two 
involves  the  execution  of  the  statistical  ARIMA  modeling  analysis  based  on 
the  data  introduced  in  stage  one.  The  output  of  the  program  is  a  detailed 
printout  of  the  feature  vectors  calculated  from  the  modeling  process  of  the 
time  series  identification.  These  vectors  are  later  used  to  evaluate  the 
correlation  of  the  MES  signature  to  the  performance  measures  from  the 
various  tasks. 

3.4  ARIMA  Model  Implementation 

The  AUTOBJ  software  package  was  used  to  extract  the  time  series 
identification  from  the  experimental  data.  It  took  approximately  30 
minutes  of  run-time  to  process  each  file  consisting  of  1000  data  points. 
The  resulting  printout  for  two  subjects  is  shown  in  Appendix  D.  The 
modeling  results  include  univariate  model  parameters  (such  as  the  mean  and 
the  ARIMA  coefficients),  residual  statistics  and  pi-weights. 
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The  first  series  of  experiments  during  phase  one  did  not  incorporate 
explicit  means  of  verifying  the  integrity  of  the  signals  being  measured. 

In  order  to  alleviate  this  problem,  a  graphics  package  was  developed  on  the 
IBM  XT  to  plot  the  data  sampled  from  both  the  frontalis  and  trapezius 
muscles. 

The  graphics  package,  written  in  C  programing  language,  is  customized  to 
the  data  formats  of  each  new  series  of  experiments.  The  first  one  hundred 
characters  of  each  data  file  is  reserved  for  a  descriptive  header.  This  is 
followed  by  1000  integers  in  ASCII  format. 

Three  variables  are  presented  to  the  user  with  this  package:  vertical 
reduction,  horizontal  magni f icati jn  and  the  number  of  data  files  to  plot. 
The  maximum  screen  resolution  that  could  be  achieved  on  the  IBM  XT  using  a 
Tecmar  Graphics  Master  Card  was  640x200.  The  flexibility  resulting  from 
the  vertical  reduction  and  horizontal  magnification  allows  the  user  to 
scale  the  graph  within  these  screen  limits.  The  vertical  reduction  allows 
the  user  to  contract  the  graph  along  the  y-axis  in  order  to  ensure  that  it 
fits  within  the  screen  boundaries.  The  horizontal  magnification  allows  the 
user  to  expand  the  graph  along  the  x-axis,  where  the  number  of  data  points 
plotted  equals  640/(horizontal  magnification).  Either  one  or  two  data 
files  are  plotted  along  with  their  respective  descriptive  headers. 

3.6  Communication  Interface 

Two  separate  communication  interfaces  are  used  in  the  experimental  setup. 
The  interface  between  the  data  acquisition  system  and  the  IBM  XT  is  used 
both  for  the  transfer  of  messages  regarding  sampling  variables  to  the  data 
acquisition  system  and  the  transfer  of  the  sampled  data  itself.  The 
interface  between  the  Apple  He  and  the  IBM  XT  is  used  for  the  control  of 
the  tasks  presented  to  the  subject. 


The  line  between  the  IBM  XT  and  the  data  acquisition  system  is  an  IEEE-488 
bus.  This  line  can  transfer  data  at  rates  up  to  three  hundred  megabytes 
per  second.  The  IEEE-488  interface  to  the  data  acquisition  system  is  built 
into  the  system  while  the  interface  to  the  IBM  XT  had  to  be  purchased 
separately.  This  consists  of  a  plug-in  card  and  related  software  which  can 
be  used  with  applications  written  in  8086  assembler,  Basic  and  C. 

The  connection  between  the  Apple  lie  and  the  IBM  XT  is  a  standard  RS-232 
line  configured  for  9600  baud,  even  parity,  1  stop  bit  and  a  7  bit  word. 

3.7  Hardware  and  Instrumentation 

3.7.1  System  Selection.  The  COSMOS  computer  initially  selected  for 
system  control  and  data  analysis  functions  was  found  to  be  too  limited  for 
our  purposes  in  a  number  of  areas.  First,  the  available  A/D  cards  for  the 
COSMOS  (used  for  amplifying  and  digitizing  the  analog  myoelectric  signals 
for  input  to  the  computer)  were  too  limiting  in  the  control  of  parameters 
such  as  amplification  gain,  sampling  frequency,  number  of  data  input 
channels,  etc.  Additionally,  the  lack  of  graphics  capability  of  the  system 
prevented  us  from  observing  the  actual  signals  acquired  for  integrity.  The 
general  unavailability  of  peripheral  devices  (e.g.,  floating-point 
hardware,  graphics  hardware)  for  this  particular  machine  initially  forced 
us  to  design  the  experiment  around  hardware  limitations.  Therefore,  a 
decision  was  made  to  re-analyze  the  hardware  requirements  and  replace  the 
COSMOS  with  a  different  system  which  offered  enough  flexibility  to  handle 
most  any  contingency  which  could  occur  over  the  course  of  the  program. 
Among  the  requirements  of  this  system  were  (1)  floating-point  hardware 
support,  and  (2)  graphic  display  capability.  Initially,  an  additional 
requirement  was  the  availability  of  an  appropriate  "C"  compiler  so  that  the 
ARIMA  model  software  used  in  the  first  phase  of  the  program  could  be  easily 
transported  to  the  new  computer.  This  requirement  later  proved  unnecessary 


since  the  computer  chosen  to  replace  the  COSMOS  had  available  to  it  a 
commercial  model  package  called  AUTQBJ  which  automated  the  modeling 
procedure  and  improved  on  some  inadequacies  of  the  original  software. 
AUTOBJ  was  suitably  modified  for  Automatic  Forecasting,  Inc.  for  use  as  a 
feature  extraction  tool  for  our  purpose. 

An  IBM  PC/XT  microcomputer  was  selected  for  both  data  acquisition  and 
executing  the  application  software.  While  this  8088  based  machine,  in 
general,  is  not  as  powerful  as  the  68000  based  COSMOS  UNIX  system,  it 
offered  the  advantage  of  widely  available  peripheral  hardware  and  software 
items  which  were  added  resulting  in  an  application-specific  configuration. 
The  system  (Figure  3-1)  was  configured  as  follows: 


FIGURE  3-1. 

SYSTEM  CONFIGURATION 


-  RS-232,  IEEE-488,  and  Centronics  interfaces 

-  8087  numeric  coprocessor 

-  Trans  Era  MOAS  7000  data  acquisition  system 

-  Tecmar  graphics  card 

-  Amdek  high  resolution  monochrome  monitor 

-  Apple  lie  with  graphics  card,  serial  port,  joystick 

-  Okidata  dot  matrix  printer 

3.7.2  Acquisition  System.  A  data  acquisition  system  (DAS)  was  required 
to  sample  the  myoelectric  signals.  The  system  had  to  impose  as  few 
limitations  as  possible  regarding  sampling  speed,  resolution,  variable 
sampling  windows,  storage,  etc.  Furthermore,  cost  was  to  be  kept  to  a 
minimum  within  program  budgetary  allocation. 

Myoelectric  signals  have  a  bandwidth  ranging  from  25  Hz  to  approximately 
5Khz.  However,  all  the  necessary  information  for  our  purposes  lies  within 
the  bandwidth  of  25  Hz  to  500  Hz.  Thus,  in  accord  with  Nyquist's  sampling 
theorem  the  DAS  was  required  to  sample  at  a  minimum  of  1000  samples  per 
second. 

The  amplitude  of  myoelectric  signal  obtained  from  a  standard  electrode 
varies  from  100  V  to  approximately  5mV.  A  system  which  is  not  sensitive 
enough  for  this  range  of  signals  would  need  to  have  the  signals 
preamplified  before  they  are  digitized.  This  in  turn  would  increase  the 
cost,  the  complexity  and  the  noise  introduced  into  the  system. 

A  minimum  of  2  analog  inputs  had  to  be  present  in  the  system.  True 
simultaneous  inputs  would  be  an  advantage  in  that  any  possible  correlation 
between  the  2  channels  could  be  investigated  with  more  accuracy. 

The  system  required  a  software  package  which  made  available  to  the  user  an 
easy  and  fast  method  of  CAlling  assembly  language  subroutines  from  C  to 
handle  the  1/0  and  sampling  functions. 


Many  systems  were  investigated  for  our  purposes.  However,  only  one  system 
fulfilled  all  the  important  requirements.  TransEra's  Model  7000-MDAS  is  a 
stand-alone  data  acquisition  system  based  on  the  Motorola  68000  micro¬ 
processor.  The  basic  processing  and  Control  Unit  consists  of:  16  slot 

chassis,  basic  power  unit,  10MHz  68000  processor  card  with  timing  module 
and  real-time  clock,  2  RS-232  serial  ports,  1  IEEE-488  interface  port,  22K 
of  static  RAM  and  a  64K  firmware  module. 

Two  upgrades  purchased  with  the  system  were  a  16-bit  A/0,  D/A  conversion 
card  which  includes  8  programmable  gains  from  1  to  128  and  a  low-level 
+0.1V,  differential  analog  input  module  which  provides  up  to  4  isolated 
simultaneous  inputs. 

This  system  has  a  resolution  of  approximately  3.05  V  per  bit  with  a  gain 
of  1  and  has  a  maximum  sampling  rate  of  approximately  48000/(number  of 
channels)  samples  per  second. 

3.7.3  Electrode  Selection.  As  discussed  in  our  interim  report  (Madni, 
Chu,  Otsubo  and  Purcell,  1984),  the  Motion  Control  electrode  packages  were 
originally  selected  for  use  based  on  the  high-performance  characteristics 
of  their  integral  signal  pre-amplifier.  However,  these  electrodes  were 
physically  too  large  to  allow  for  placement  on  some  potentially  useful 
sites  (e.g.,  forehead).  Additionally,  the  sensitive  CMOS  circuitry  used  in 
the  design  of  the  pre-amplifier  was  extremely  susceptible  to  static 
discharge  induced  failure  with  no  simple  means  of  failure  indication  short 
of  experiencing  signal  degradation. 

It  was  decided  that  si  1 ver/si 1 ver  chloride  "button"  type  electrodes  along 
with  an  external  high-performance  pre-ampl if ier  would  be  more  suitable  for 
our  purposes.  The  characteristics  of  the  VIVO  Metric  si  1 ver/si lver 
chloride  electrodes  in  summarized  in  Table  3-2. 


TABLE  3-2 

CHARACTERISTICS  OF  VIVO  METRIC 
SILVER/SILVER  CHLORIDE  ELECTRODES 


-  Sensor  Dia. : 

4mm 

-  Housing  Dia.: 

7.2mm 

-  Overall  Height: 

6mm 

-  Electrode  Cavity: 

1mm  deep 

-  Lead: 

lm  shielded  cable 

3.7.4  Electromagnetic  Interference.  Our  preliminary  study  had  indicated 
that  60  Hz  interference  was  contaminating  the  myoelectric  signal.  To  this 
end,  we  constructed  a  Faraday  cage  (a  grounded,  wire-mesh  cage)  within  the 
room  in  which  a  subject  could  sit  while  performing  a  task,  essentially 
shielding  him  from  extraneous  R-F  interference.  Any  remaining  EMI  was 
removed  by  using  differential  electrodes  and  by  grounding  the  subject. 


4.  EXPERIMENTAL  STUDY 


4.1  Overview 

The  phase  one  series  of  experiments  were  designed  to  assess  the  "goodness" 
of  stochastic  model-derived  features  pertaining  to  operator  alertness  and 
workload  levels.  In  that  study,  behavior  variability  and  hardware 
limitations  prevented  totally  reliable  measures  to  be  gathered  from 
subjects.  To  alleviate  these  problems,  more  stringent  experimental  design 
and  methodologies  were  developed  to  improve  experimental  control  and 
circumvent  the  hardware  limitations  observed  during  the  initial  series  of 
experiments. 

This  section  describes  the  new  experimental  procedures  and  the  tasks 
involved  in  deriving  workload  correlates.  The  rationale  for  changing  the 
tasks,  variables  and  procedures  from  the  original  series  of  experiments, 
was  so  that  workload  performance  and  perceived  stress  could  be  measured 
more  reliably  with  model-derived  MES  features. 

4.2  Experimental  Hypotheses  and  Test  Procedure 

The  major  effort  of  this  program  was  to  be  able  to  predict  perceived 
workload  levels  from  model-derived  MES  features.  To  predict  workload 
levels,  these  MES  features  have  to  be  reliable  and  diagnostical  to  infer 
appropriate  operator  states.  Two  hypotheses,  which  are  characterized  as 
"reliability"  and  "diagnosticity,"  were  investigated  in  this  study.  These 
hypotheses  are: 

o  Reliability  -  A  minimum  of  one  set  of  values  in  the  ARIMA  model 
provides  invariant  or  near-invariant  pattern  values  for  each 
subject  within  a  predetermined  underlying  level  of  mental  load 
for  a  given  task  category,  thus  achieving  "reliability". 


o  Diagnosticity  -  Each  level  of  task  difficulty  has  its  own 
distinct  pattern,  providing  a  "diagnostic"  feature  for  predicting 
operator  wokload/altertness  with  model-derived  MES  features. 

To  ascertain  the  reliability  and  diagnosticity  of  MES-derived  features, 
several  issues  were  investigated.  The  first  involved  detecting  at  least 
one  set  of  values  in  the  ARIMA  model  for  each  underlying  loading/alertness 
level,  subject,  muscle  site,  and  task  level.  To  determine  this,  model 
coefficients  and  their  related  pi-weights  were  examined  for  invariance 
across  multiple  samples  within  trials  and  across  multiple  trials.  If 
values  were  found  that  satisfy  the  necessary  invariant  conditions,  they 
were  then  considered  "reliable".  These  reliable  features  must  also  have 
been  sufficiently  different  in  magnitude  for  each  level  of  the  task,  to 
exhibit  the  attribute  of  "diagnosticity." 

4.3  Experimental  Tasks 

The  original  series  of  experiments  were  based  on  simulating  a  subset  of  the 
Criterion  task  Set  (CTS)  workload  test  battery  that  was  developed  at  AFAMRL 
(Shingledecker,  1983).  These  experiments  were  based  on  the  degree  to  which 
they  satisfied  the  requirements  of:  (I)  validity  and  reliability,  (2)  flex¬ 
ibility  and  quantifiabil ity,  (3)  memory,  (4)  mental  mathematics/reasoning, 
and  (5)  choice  reaction  time.  Unfortunately,  data  from  only  one  subject 
was  used  to  validate  task  parameters.  Thus,  the  tasks  used  in  the  original 
experiments  were  not  stringently  pretested  and  validated  to  correlate  to 
the  desired  levels  of  "low"  and  "high"  workload  levels;  the  low  difficulty 
level  tasks  may  have  been  too  easy  for  the  subjects,  generating  boredom, 
and  the  high  levels  may  not  have  actually  represented  a  higher  task-imposed 
of  "stress".  There  were  observed  differences  in  performance  and  subjective 
ratings  between  the  high  and  low  levels  of  task  difficulty,  but  all 
subjects  seem  to  exhibit  boredom  and  fatigue  due  to  the  long  experimental 
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sessions  and  the  repetitive  task  performance  required  of  them.  Thus,  in 
the  pilot  set  of  experiments,  instead  of  experiencing  different  levels  of 
stress,  subjects  were  experiencing  boredom  and  fatigue. 

Additionally  the  original  tasks  were  discrete,  independent  task  components, 
rather  than  integrated  task  situations  to  match  the  Criterion  Task  Set 
(CTS)  under  development  at  AFAMRL.  This  minimized  the  external  validity  of 
the  tasks. 

To  alleviate  these  problems,  a  new  set  of  tasks  was  developed  to: 

o  Minimize  the  fatigue  and  boredom  that  resulted  from  the  first 
series  of  experiments; 

o  Achieve  a  higher  level  of  external  validity; 
o  Represent  validated  task  levels  of  low  and  high  stress. 

The  new  set  of  tasks  was  developed,  based  on  the  objective  that  when  given 
a  primary  task  at  a  comfortable  pace,  adding  a  secondary  task  will  increase 
the  level  of  difficulty  of  the  total  task.  Thus,  the  perceived  stress 
experienced  by  the  subject  should  theoretically  increase.  We  first 
analyzed  the  basic  CTS  and  selected  a  subset  to  employ  as  our  experimental 
tasks. 

The  overall  CTS  included  the  following  task  components: 

(1)  Perceptual  tasks 

o  Probability  monitoring  task 
o  Auditory  monitoring  task 
o  Visual  target  search  task 
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(2)  Central  Processing  Tasks 


o  Memory  tasks  -  memory  update,  memory  recall 
o  Manipulation  '■nd  comparison  tasks  -  linguistic  processing, 
mathematical  computation,  spatial  pattern  identification 
o  Reasoning  tasks  -  analogical  reasoning  and  grammar 
o  Planning  and  scheduling  -  flight  assessment  and  supervisory 
control 

(3)  Motor  tasks 

o  Critical  tracking  task 

A  subset  of  the  CTS  was  selected  from  existing  tasks  previously  researched 
and  validated  at  AFAMRL.  These  task  dimensions  included  a  motor  task 
(critical  tracking)  and  a  perceptual  task  (auditory  monitoring). 

The  two  tasks  used  in  this  series  of  experiments  are  described  in  the 
following  paragraphs. 

Critical  Tracking  Task 

The  critical  Tracking  Task  was  developed  as  the  primary  task.  In  this  task 
the  subjects  were  required  to  track  a  target  on  a  computer  display  screen 
by  controlling  a  "site"  object  with  a  joystick.  The  subject  was  directed 
to  track  and  surround  the  target  with  the  site,  as  the  target  randomly 
moved  around  the  screen.  Performance  was  based  on  the  average  distance 
between  the  site  and  target  over  the  entire  session.  The  more  control  the 
suuject  maintained  over  the  tracking  task,  the  lower  the  average  distance 
recorded.  This  task  was  developed  to  represent  the  desired  level  of  "low" 
stress. 


Critical  Tracking  Task  With  Added  Auditory  Monitorin' 


II 

The  same  Critical  Tracking  Task  was  then  combined  with  a  secondary  auditory 
monitoring  task.  This  second  task  required  the  subject  to  respond  to  an 
auditory  tone  by  pressing  a  button  on  top  of  the  joystick.  The  tones  were 
intermixed  so  that  one  tone  represented  signals  and  a  different  tone 
represented  "noise".  The  subject  was  required  to  distinguish  the  signals 
from  the  "noise"  by  pressing  the  button  only  when  the  signal  was  heard. 
The  subject  was  required  to  respond  within  1/2  second  from  the  time  the 
signal  was  presented.  If  no  response  was  made  within  that  time,  an 
omission  error  (a  "miss")  was  recorded.  If  the  subject  pressed  the  button 
when  no  tone  was  presented,  a  commission  error  (a  "false  alarm")  was 
recorded.  This  monitoring  task  was  presented  to  the  subject  with  the 
primary  critical  tracking  task,  thus  representing  a  more  stressful  task 
than  the  primary  task  alone. 

4.4  Experimental  Variables 

* 

The  experimental  design  utilized  a  within-subject  repeated  measures  design. 
The  following  independent  variable  was  tested: 

o  Task  loadinq  -  two  levels  of  difficulty  or  "stress" 

(1)  Critical  tracking  task  only  ("low"  stress) 

(2)  Critical  tracking  task  +  auditory  monitoring  tasks  ("high" 
stress) 

The  dependent  variables  gathered  from  the  varying  task  loading  situations 


were: 


_  o  Performance  measures 

a 

(1)  average  distance  between  the  site  and  target  for  the  critical 
tracking  task 

(2)  auditory  monitoring  responses,  (errors  of  omission  and 

'  ccmission) 

o  Subjective  measures  of  effort  and  stress  -  gathered  from 
post-test  questionnaires 

o  MES  data  gathered  through  the  data  acquisition  system 

In  the  initial  series  of  experiments  there  were  two  task  categories  and  two 
levels  of  task  difficulty. 

o  Task  category  -  two  types 

(1)  Central  processing  tasks 

®  (2)  Motor  tracking  tasks 

o  Task  loading  -  two  levels  of  difficulty 

f  (1)  Low  load 

(2)  High  load 

We  attempted  to  stay  within  the  same  criteria  established  in  the  first 
series,  i.e.,  CTS  implementation  utilizing  two  levels  of  stress.  The  only 
major  change  was  the  integration  of  the  task  categories  of  a  perceptual 
(auditory)  and  a  motor  tracking  task  together,  to  achieve  a  higher  level  of 
external  validity.  We  also  employed  the  use  of  secondary  task  loading  to 
validate  the  assumption  that  presenting  an  additional  task  places  higher 
cognitive  load  on  the  subject,  thus  increasing  his  perceived  "stress." 
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This  probe-reaction-time  technique  has  been  utilized  with  notable  success 

< 

in  uncovering  resource  demands  for  various  perceptual /motor  and  cognitive 
tasks  (Ells,  1973;  Posner  &  Bois,  1971;  Posner  &  Keele,  1969). 

To  validate  the  two  levels  of  "stress",  pilot  studies  were  performed. 
Initially,  five  test  subjects  were  presented  seven  varying  levels  of 
difficulty  for  the  critical  tracking  task.  They  were  given  30  seconds  to 
perform  each  level,  from  the  easiest  (Level  1)  to  the  most  difficult  (Level 
7).  Each  test  subject  practiced  and  became  familiar  with  all  the  levels. 
Each  was  then  given  the  opportunity  to  move  the  level  of  difficulty  up  or 
down  by  pressing  the  up-cursor  (increases  the  difficulty  level)  and 
down-cursor  (decreases  the  difficulty  level)  on  the  computer  keyboard. 
They  were  told  to  reach  the  level  at  which  they  felt  most  comfortable,  i.e. 

g 

not  too  easy  and  not  too  difficult.  We  assumed  that  this  level  would 
represent  the  "low"  stress  level.  Three  of  the  test  subjects  chose  level 
4,  while  two  subjects  chose  Level  5,  thus  we  established  Level  5  as  the 
level  of  "low"  stress,  to  minimize  the  problem  of  boredom  and  fatigue  that 

« 

was  found  to  occur  in  the  first  series  of  experiments.  We  again  tested  the 
pilot  subjects  on  the  task  only  at  Level  5,  and  all  agreed  that  it  was 
comfortable  level.  Those  who  initially  selected  Level  4  as  the  comfortable 
level  could  not  tell  the  difference  once  given  the  task  at  Level  5,  thus, 

< 

we  maintained  Level  5  as  the  established  level  of  "low"  stress. 

To  find  the  appropriate  level  for  the  secondary  task  of  auditory 
monitoring,  we  again  conducted  pilot  studies  in  which  the  same  test 
subjects  performed  the  same  critical  tracking  task  (Level  5)  with  varying 
levels  of  auditory  monitoring  (Levels  1-7).  Test  subjects  again  were  able 
to  practice  the  critical  tracking  plus  auditory  monitoring  task  to  become 
familiar  with  all  the  levels.  Next  they  were  given  the  opportunity  to 
control  the  levels  of  auditory  monitoring  by  pressing  the  up-cursor 
(increasing  frequency  of  signal  and  noise  presentation)  and  the  down-cursor 
(decreasing  frequency  of  signal  and  noise  presentation).  Subjects  were 

« 


4-7 


told  to  maintain  the  level  at  which  they  felt  comfortable,  i.e.  not  too 
easy  and  not  too  difficult.  All  subjects  chose  Level  4.  This  task  was  to 
represent  the  level  of  "high"  stress,  thus  we  increased  the  actual  auditory 
monitoring  task  to  Level  5.  We  attempted  to  increase  the  difficulty  by  two 
levels  (to  Level  6),  but  discovered  that  the  signals  were  presented  in  a 
predictable  pattern.  this  made  it  easier  instead  of  more  difficult  to 
respond.  Consequently,  we  maintained  the  auditory  presentation  of  Level  5 
as  the  established  level  of  "high"  stress. 

We  then  presented  both  tasks  to  all  test-subjects,  gathered  performance  and 
subjective  data,  and  found  that  there  was  an  observable  difference  between 
the  "low"  stress  and  "high"  stress  task  performances  based  on  the  average 
distance  from  site  to  target  on  the  critical  tracking  tasks,  and  on  the 
subjective  responses  to  the  amount  of  effort  exerted  and  stress  perceived. 

4.5  Muscle  Site  Selection 


Exist. ng  literature  (Madni ,  Chu,  Otsubo,  and  Purcell,  1984)  suggests  a 
number  of  muscles  (frontalis,  trapezius,  splenius,  temporalis,  masseter, 
brachioradial is)  which  undergo  sporadic,  potentially  diagnostic  activity 
during  high  workload  situations.  Our  new  data  acquisition  system  gave  us 
the  capability  of  acquiring  data  from  multiple  muscle  sites,  and  our  new 
electrodes  were  of  a  small  enough  size  to  permit  judicious  placement. 
However,  the  additional  processing  time  involved  in  modeling  the  data  from 
each  additional  muscle  restricted  us  at  this  time  to  a  maximum  of  two 
sites.  The  frontalis  and  trapezius  muscles  were  chosen. 

4.6  Subjects  and  Procedures 

In  this  phase  of  the  program,  12  male  subjects  were  recruited  from  the 
local  university,  and  from  the  available  Perceptronics  pool  of  volunteer 
subjects.  Subject  ages  ranged  from  19-35  years  old.  All  had  at  least  a 
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high  school  diploma  and  some  experience  with  using  computers.  A  2X2 
within-subject  repeated-measures  design  was  employed.  Each  subject 
performed  the  two  tasks  twice.  After  they  completed  each  task,  they  were 
given  a  questionnaire  (See  Appendix  C)  to  complete  regarding  their 
subjective  evaluations  of  the  effort  expended  and  the  difficulty  of  the 
tasks. 

In  this  phase  of  the  program,  several  major  changes  were  implemented. 
Before  performing  the  actual  experiment,  several  pilot  studies  were 
conducted  to  investigate  the  validity  of  the  tasks  to  be  presented.  Upon 
examining  the  results  of  the  pilot  studies,  we  were  able  to  substantiate 
the  expected  differences  between  the  two  levels  of  task  difficulty  based  on 
task  performance  and  perceived  stress. 

The  original  series  of  experiments  attempted  to  represent  low  and  high 
levels  of  difficulty  by  isolated,  discrete  tasks.  By  presenting  subjects 
with  different  tasks,  there  was  no  way  of  validating  whether  the  tasks  used 
were  actually  of  varying  degrees  of  difficulty  or  whether  they  were  just 
different  in  task  type,  that  is,  it  was  more  clear  whether  the  tasks 
really  fell  as  different  points  on  "level  of  difficulty"  continuum.  To 
correct  for  this  undefined  task  variability,  the  method  of  secondary  task 
loading  was  used;  an  auditory  monitoring  task  was  added  to  the  primary 
critical  tracking  task.  This  procedure  insured  that  if  there  was  a 
difference  between  the  primary  task  ("low"  stress  task)  and  the  primary 
plus  secondary  task  ("high"  stress  task),  then  it  must  be  due  to  the 
additional  cognitive  load  imposed  by  the  secondary  task  (see  Knowles, 
1963). 

We  also  recruited  12  subjects  participate  in  the  experiment,  compared  to 
only  3  subjects  in  the  initial  phase  so  that  the  statistical  power  of  the 
experiment  would  be  increased.  The  results  from  only  three  subjects  may 


have  been  due  to  chance.  Increasing  the  number  of  subjects  would  provide  a 
greater  probability  that  the  results  would  demonstrate  statistically 
significant  differences. 

The  total  time  required  of  each  subject  was  1  to  1-1/2  hours.  This  is 
vastly  different  from  the  two  to  three  sessions  required  of  them  in  the 
initial  series  of  experiments.  Reducing  the  amount  of  experimental  time 
minimized  the  probability  that  the  subjects  were  bored  and  fatigued  with 
the  experimental  procedure,  and  consequently,  reduced  the  probability  that 
physiological  measures  would  be  confounded. 

Each  subject  was  asked  to  read  the  instructions  (Appendix  A)  explaining  the 
experimental  procedure.  Subjects  were  encouraged  to  ask  questions  as  they 
read  the  instructions  to  ensure  that  they  understood  the  tasks.  They  were 
then  requested  to  read  and  fill  out  a  "personal  information  fact  sheet"  and 
"consent  to  act  as  an  experimental  subject"  form.  The  active  electrode 
assembly  was  then  attached  above  the  (frontalis  muscle),  and  to  the 
subject's  upper  back  (trapezius  muscle),  running  parallel  to  the  muscle 
fibers.  A  third  "ground"  electrode  was  positioned  near  the  elbow.  All 
electrodes  were  placed  on  the  subjects'  non-dominant  side  (right-handed 
subjects  had  the  electrodes  positioned  on  the  left  arm,  and  vice  versa). 
Subjects  were  given  an  orientation  and  practice  session  lasting  5-10 
minutes  to  minimize  the  learning  effect  that  ^ould  occur  during  the 
experiment  itself.  The  practice  session  was  concluded  when  the  subjects 
produced  comparable  scores  on  two  successive  trials  for  each  task,  the 
subjects  were  then  instructed  to  sit  comfortably  upright  in  the  chair  and 
cautioned  against  moving  the  side  of  their  bodies  to  which  the  electrodes 
were  attached.  At  this  time,  a  sample  time  series  of  data  was  recorded  for 
subjects  representing  their  "At  rest"  state.  At  the  end  of  each  task, 
subjects  were  asked  to  fill  out  a  questionnaire  of  subjective  ratings  and 
post-experimental  comments.  After  completing  this  form,  subjects  rested 
for  a  few  minutes  before  the  next  task  was  presented.  All  subjects  were 
given  two  trials  for  each  task  in  a  counterbalanced  order. 


Each  experimental  task  lasted  200  seconds.  During  each  trial,  data  sampled 
at  1  Khz  was  collected  in  1000  msec  windows  spaced  30  seconds  apart  in  the 
last  45  seconds  of  the  trial.  This  data  collection  scheme  allowed  us  to 
evaluate  feature  reliability  both  within  and  between  trials. 

4.7  Performance  Measures 

The  performance  measure  of  average  distance  between  site  and  target  was 
gathered  so  that  one  common  performance  measure  could  be  compared  between 
both  tasks.  In  the  original  series  of  experiments,  the  performance 
measures  for  low  and  high-stress  tasks  differed,  making  it  difficult  to 
reliably  determine  the  meaningfulness  of  performance  differences  observed. 

The  performance  measures  collected  during  the  experimental  trials  were: 

Critical  Tracking  Task  - 

o  Average  distance  between  site  and  target  -  the  average  distance 
measured  over  10  second  intervals  between  the  random  target  and 
the  site  controlled  by  the  subject 

Critical  Tracking  Task  +  Auditory  Monitoring  Task  - 

o  Average  distance  between  site  and  target  -  the  average  distance 
measured  over  10  second  intervals  between  the  random  target  and 
the  site  controlled  by  the  subject 

auditory  response  accuracy  -  the  correct  responses  (hits), 
omissions  (misses)  and  commissions  (false  alarms)  made  by 
subjects  to  the  auditory  "signals" 


These  performance  results  were  recorded  and  printed  for  each  task  (Appendix  D) 

Subjects  were  also  requested  to  fill  out  a  questionnaire  after  each  task 
was  completed,  the  questionnaire  required  that  subjects  rate  (on  a  sliding 
scale)  various  items  related  to:  (1)  their  perception  to  task  difficulty, 
and  (2)  their  perceived  level  of  effort  (Appendix  C). 

The  results  of  the  performance  measures  and  subjective  ratings  along  with 
the  statistical  and  pattern  analysis  are  discussed  in  Section  5.2. 


5.  EXPERIMENTAL  RESULTS  AND  DATA  ANALYSIS 


5.1  Overview 

This  section  discusses  the  subjective  responses,  objective  performance  data 
and  objective  MES  measurements  that  were  collectively  analyzed  in  order  to 
identify  ARIMA  features  for  varying  levels  of  stress  and/or  operator 
workload.  Included  is  a  discussion  of  the  trends  between  tasks  that  were 
identified  and  conclusions  that  were  drawn  from  the  ARIMA  model-based 
analysis. 

5.2  Task  Performance  and  Subjective  Ratings 

Analysis  of  task  performance  was  made  by  comparing  the  average  distance 
between  the  site  and  target  between  the  two  tasks.  These  results  indicate 
that  an  observed  difference  exists  between  performance  on  the  two  tasks 
(Figure  5-1).  Despite  the  fact  that  these  differences  were  statistically 
significant  (paired  5-test,  p  >  .05),  75%  of  the  subjects'  performance 
dropped  from  the  single  to  the  dual  task  situation.  With  an  increased 
number  of  subjects,  and  perhaps  a  greater  discrepancy  in  difficulty  levels 
of  the  single  and  dual  tasks,  it  is  reasonable  to  expect  that  statistical 
significance  would  have  supported  the  resulting  observed  performance 
differences.  The  observations  showed  that  9  subjects  out  of  the  12 
performed  worse  on  the  dual  task  (critical  tasking  +  auditory  monitoring 
task).  This  observed  difference  may  be  an  indication  that  the  dual  task 
may  have  been  perceived  as  more  stressful. 

Questionnaire  responses  were  then  analyzed  to  ascertain  differences  between 
the  "perceived  stress"  levels  on  the  two  tasks.  The  pertinent  questions 
relating  to  perceived  stress  were: 
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FIGURE  5-2(c). 

QUESTIONNAIRE  RESPONSES  TO  "HOW 
MUCH  STRESS  DID  YOU  EXPERIENCE? 


o  "What  level  of  effort  did  you  expend?" 
o  "How  challenging  was  the  session?" 

o  "How  much  strain  did  you  experience  during  the  session" 

Comparisons  between  the  two  tasks  exhibited  significant  differences  in 
response  for  these  questions  (paired  t-test,  p  <  .005;  Figure  5-2).  These 
differences  indicate  that  increased  "perceived  stress"  was  experienced 
during  the  dual  task,  even  though  there  were  no  statistically  significant 
differences  in  task  performance. 

5.3  Results  of  MES  Feature  Extraction 

MES  sampling  was  conducted  a  total  of  nine  times  for  each  subject,  once  at 
rest  before  any  of  the  tasks  began  and  twice  during  each  task  for  both 
muscles.  Each  file  of  MES  data  was  stored  and  subsequently  plotted  on  the 
CRT.  Samples  of  these  measurements  are  shown  in  Figures  5-3  and  5-4  for 
subjects  JOJ  and  AXC. 

The  ARIMA  analysis,  a  result  of  the  AUTOBJ  software  implementation,  was 
conducted  for  each  of  the  files  stored.  Sample  outputs  for  subject  JJJ  and 
AXC  are  shown  in  Appendix  0.  Modeling  results  included  univariate  model 
parameters,  residual  statistics  and  pi-weights  for  each  file.  The 
univariate  model  parameters  included  autoregressi ve,  moving  average  and 
trend  constant  parameters.  For  each  of  these  a  factor,  lag,  coefficient 
and  F-ratio  value  was  extracted.  The  residual  statistics  included  the  sum 
of  squares,  mean  square,  R-squared,  degrees  of  freedom  and  number  of 
residuals.  Univariate  models  were  expressed  as  pi-weights  (weighted  sum  of 
the  past  plus  a  random  shock).  The  general  form  of  the  model  expressed  by 
the  pi-weights  is: 

Z ( t )  =  Constant  +  P i ( 1 )  ’  Z(t-l)  +  Pi(2)  ’  Z ( t - 2 )  ...  +  Pi(n)  *  Z(t-n) 
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where: 


n  is  number  of  lags 
t  is  time 

Pi(n)  is  the  value  of  lag  at  n 
Z(t)  is  the  pi-weight  at  time 
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6.  DISCUSSION  OF  FINDINGS 
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Based  on  the  data  collected  from  the  Box  &  Jenkins  ARIMA  analysis,  we  were 
able  to  compare  the  MES  features  between  tasks,  between  trials,  between 
muscle  sites  and  between  the  rest  period  and  each  of  the  two  tasks. 

Trends  were  found  in  two  individuals  as  a  result  of  careful  analysis  of  the 
data.  Subject  JJJ  exhibited  a  definite  trend  with  the  frontalis  muscle 
between  tasks  1  and  2.  Table  6-1  shows  the  significant  pi-weights  for  this 
individual  with  respect  to  the  task,  window  and  trial.  A  conclusion  which 
can  be  drawn  from  this  information  is  that  the  pi-weights  are  generally 
significant  for  task  2  and  not  for  task  1.  Subject  AXC  exhibited  a  trend 
between  tasks  1  and  2  also  but  with  the  trapezius  muscle.  Table  6-2 
identifies  the  significant  pi-weights  with  respect  to  the  task,  window  and 
trial.  Pi-weights  for  this  subject  were  generally  significant  in  task  1 
and  not  in  task  2. 

The  subjective  data  indicate  a  statistically  significant  difference  in 
perceived  stress  between  the  single  and  dual-task  groups.  These  results 
are  supported  by  observed  differences  in  performance  for  9  of  the  12  (or 
75%  of)  the  subjects.  These  performance  differences,  although  in  the 
expected  direction,  did  not  reach  statistical  significance. 

The  first  obvious  explanation  for  the  lack  of  statistical  significance 
between  the  single  and  dual-task  performance  measures  is  the  relatively 
small  number  of  experimental  subjects  tested.  It  is  clear  that  the  greater 
the  number  of  subjects  tested,  the  greater  the  power  of  the  experimental 
manipulations,  and  consequently,  the  greater  the  probability  that  small 
differences  may  reach  statistically  significant  levels.  Since  we  knew  that 
a  large  number  of  subjects  would  not  be  accessible,  given  the  practical 
constraints  under  which  this  study  was  conducted,  we  employed  a  within- 
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TABLE  6-1 


N 


SIGNIFICANT  Pi -WEIGHTS  (JJJ) 


Subject:  JJJ  (Frontalis) 


Significant 

Pi-Weights 

Task 

Window 

Trial 
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1 

1 

1 

None 

1 

2 

1 

None 

1 

1 

2 

None 

1 

2 

2 

None 

2 

1 

1 

5,6 

2 

2 

1 

5, 6, 7, 8 

2 

1 

2 

None 

2 

2 

2 

5 

TABLE  6-2 

SIGNIFICANT  P-j-WEIGHTS  (AXC) 


Subject:  AXC  (Trapezius) 


Significant 
Pj -Weights 


Task 

Window 

Trial 

20-25 

1 

1 

1 

20,21,22,23 

1 

2 

1 

20,22,23,24,25 

1 

I 

2 

20,21,22,23 

1 

2 

2 

20,21,22,23 

2 

1 

1 

20,21,22,23,24,25 

2 

2 

1 

20,21,22,23,24,25 

2 

1 

2 

None 

2 

2 

2 

None 
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subjects  experimental  design,  which  generally  provides  a  more  powerful 
experiment  with  less  subjects  than  a  between-groups  design.  It  is 
apparent,  however,  that  even  under  those  conditions,  the  number  of  subjects 
tested  was  inadequate.  The  fact  that  three-fourths  of  the  subjects  tested 
did  show  performance  differences  between  the  single  and  dual-task 
situations  does  suggest,  however,  that  using  a  greater  number  of  subjects 
may  have  produced  results  that  were  statistically  significant. 

The  experimental  tasks  used  in  this  study  were  selected  because  they  are 
classical  tasks.  That  is,  motor  tracking  and  auditory  vigilance  tasks,  in 
combination,  are  used  throughout  the  skilled  performance  literature  to 
"load"  the  subject  and  increase  task  difficulty  levels.  It  is  well 
documented,  the  time-sharing  literature,  that  imposing  a  secondary  task  on 
a  user  increases  the  attentional  demands  on  subjects.  In  theory,  then,  the 
dual-task  is  more  difficult  to  perform  than  the  single  task  alone.  This 
has  been  noted  to  increase  the  level  of  stress  under  which  the  subject 
performs.  Because  the  tasks  used  were  chosen  based  on  their  successful  use 
in  similar  research  in  the  past,  the  lack  of  statistical  significance  of 
the  tracking  performance  measure  in  this  study  requires  an  alternative 
explanation. 

A  more  theoretical  but  certainly  plausible  explanation  for  this  lack  of 
statistical  significance  may  be  advanced  based  on  a  cognitive  model  of 
allocation  of  attention  (see  Kahneman,  1973).  While  performing  the  dual¬ 
task,  the  subjects  could,  in  theory,  have  allocated  their  attention  to  the 
primary  tracking  task,  even  though  they  were  instructed  that  maintaining 
performance  on  both  primary  and  secondary  tasks  was  critical.  Because  the 
primary  (tracking)  task  was  inherently  more  "entertaining"  (more  like  a 
video  game)  and  hence  perhaps  more  motivating  than  the  secondary  (signal 
detection  task),  it  is  distinctly  possible  that  subjects  ignored  the 
instructions  and  allocated  their  attention  to  the  primary  tracking  task. 
That  is,  subjects  may  well  have  allocated  just  enough  attention  to  the 


secondary  task  of  signal  detection  to  perform  to  their  own  levels  of 
perceived  adequacy,  rather  than  trying  to  perform  well  on  both  tasks.  If 
this  was,  in  fact,  the  case,  then  it  would  be  logical  that  some  decrease  in 
primary  task  performance,  and  hence  observed  differences,  would  be  the 
result,  but  that  those  observed  differences  would  not  have  been  profound 
enough  to  reach  statistical  significance.  This  may  also  explain  why  the 
differences  in  MES  readings  associated  with  the  dual-task  situation  also 
were  not  very  pronounced. 

To  determine  whether  subjects  did,  in  fact,  "satisfice"  (i.e.,  perform  to 
their  perception  of  adequacy)  in  performing  the  secondary  task,  one  could 
look  at  the  signal  detection  performance  of  subjects  in  the  dual -task 
situation.  Although  no  "control"  group  performing  signal  detection  only 
was  employed  in  this  study,  one  could  still  expect  that  if  subjects  were 
allocating  enough  attention  to  the  signal  detection  task,  performance 
should  have  been  above  chance.  Stated  differently,  signal  detection 
performance  at  or  below  change  levels  for  the  subjects  in  the  dual-task 
group  would  indicate  that  the  subjects  did  not  allocate  to  the  secondary 
task  the  attention  it  demanded. 

To  determine  the  viability  of  this  explanation,  a  post-hoc  analysis  was 
performed  on  the  raw  signal  detection  data  from  the  subjects  in  the  dual- 
task  situation.  The  results  indicate  that  all  subjects  performed  below 
chance  levels.  On  the  average,  subjects  signal  detection  responses  were 
only  correct  13%  of  the  time  (see  Table  6-3).  Three  subjects  had  more 
incorrect  than  correct  responses  (14-18%  more  incorrect  responses),  and 
only  two  subjects  managed  to  achieve  over  40%  in  correct  responses.  This 
may  suggest  that  the  subjects  performed  the  dual-task  at  their  level  of 
adequacy  by  allocating  less  attention  to  the  secondary  task  then  it  called 
for,  and  thus,  maintaining  a  consistent  level  of  performance  for  the 
primary  task. 


6-5 


TABLE  6-3 


SECONDARY  TASK  AVERAGE  CORRECT  AND  INCORRECT 
RESPONSES  FROM  TOTAL  SIGNALS  PRESENTED 


SUBJECTS 

TRIAL  1 

TRIAL  2 

AVERAGE 

WG 

14% 

23% 

19% 

(more  correct  responses) 

NI 

-10% 

-17% 

-14% 

(more  incorrect  responses) 

BN 

41% 

40% 

41.5% 

JM 

20% 

34% 

27% 

GN 

-19% 

-16% 

-18% 

RR 

16% 

24% 

20% 

DS 

50% 

41% 

46% 

AC 

-2% 

4% 

1% 

JJ 

-16% 

-21% 

-19% 

FH 

-37% 

-19% 

-28% 

DN 

33% 

31% 

32% 

JF 

12% 

18% 

15% 

Summary:  On  the  average,  subjects  achieved  13%  more  correct  responses  than 
incorrect  responses  (see  Section  6.0  for  a  discussion) 
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Appendix  A 


Software  Specifications 


Overview 


The  two  software  packages  involved  in  the  myoelectric  signal  acquisition 
consist  of  one  system  for  the  IBM  XT  and  one  for  the  Apple  lie.  The  IBM  XT 
is  used  to  control  the  newly  acquired  data  acquisition  system,  the 
experimental  parameters  and  the  Apple  lie.  The  Apple  lie  is  used  to 
present  the  tasks  to  the  subjects,  to  record  subject  responses  and  to  print 
a  subject  performance  record. 

Apple  lie  System 

There  are  two  tasks  present  on  the  Apple  He  which  are  both  written  in 
Apple  GraForth.  When  not  performing  either  task  or  printing  out  a 
performance  record,  the  Apple  waits  for  a  task  number  to  be  received  from 
the  IBM. 

Task  l  is  a  straightforward  tracking  task.  This  task  consists  of  a 
crosshair  which  moves  randomly  at  a  constant  rate  about  the  screen  along 
both  x  and  y  axes.  The  subject  has  joystick  control  over  a  circle  on  the 
screen.  The  object  is  to  keep  the  crosshair  as  close  to  the  center  of  the 
circle  as  possible  for  the  3.5  minutes  duration  of  the  task. 

Task  2  consists  of  task  1  along  with  a  simultaneous  auditory  task.  At 
random  intervals,  a  short  tone  is  emitted  from  the  Apple's  speaker.  Sixth 
per  cent  of  the  time,  this  tone  will  be  high  pitched.  At  this  point  the 
subject  has  exactly  0.5  seconds  to  respond  to  the  tone  by  pressing  the 
joystick  button.  The  system  then  responds  with  a  much  higher  tone  to 


acknowledge  the  button  press.  Forty  per  cent  of  the  time  however,  the  tone 
will  be  low  pitched.  Should  the  subject  respond  to  this  lower  tone,  it  is 
recorded  as  an  error. 


A  restriction  is  placed  on  the  high  tones.  Once  a  high  tone  is  emitted 
there  will  not  be  another  high  tone  for  at  least  0.5  seconds. 
Nevertheless,  any  number  of  low  tones  can  be  sounded  within  this  period. 
This  is  done  to  allow  the  subject  time  to  respond  to  the  tone. 

Also,  once  the  joystick  button  is  pressed,  it  is  not  polled  again  for  the 
next  0.5  second.  This  is  done  in  order  to  remove  the  possibility  of  a 
double  press  for  a  given  tone. 


Record inc 


Performance 


In  both  tasks,  a  running  average  of  the  distance  from  the  center  of  the 
civJe  to  the  center  of  the  crosshair  is  calculated.  Every  ten  seconds 
this  is  stored  in  a  buffer  and  the  average  distance  is  reset  to  zero. 


In  task  2,  responses  to  the  tones  are  also  recorded.  These  are  groupeu 
into  three  categories:  good  responses,  false  responses  and  no  responses. 
A  good  response  occurs  when  the  joystick  button  is  pressed  within  0.5 
seconds  after  a  high  tone  is  emitted.  A  false  response  is  recorded  when 
the  joystick  button  is  pressed  but  there  was  no  high  tone  emitted  within 
the  previous  0.5  seconds.  Finally,  a  no  response  occurs  when  0.5  seconds 
passes  after  j  high  tone  is  emitted  without  the  joystick  button  being 
pressed.  As  with  the  average  distance,  every  ten  seconds,  the  tone 
responses  are  stored  in  a  buffer  and  reset  to  zero. 


Upon  terminating  the  cask,  the  Apple  prints  its  response  buffers  and 
returns  to  its  initial  waiting  state. 


IBM  XT  System 

The  system  on  the  IBM  XT  is  completely  written  in  C  and  consists  of  four 
primary  functions. 

o  Control  the  Apple  lie. 

o  Interface  with  the  data  acquisition  system, 
o  Plot  the  data  files  received  from  the  data  acquisition  system, 
o  Interface  with  the  experimenter . 

Control  of  the  Apple  He 

The  IBM  communicates  with  the  Apple  He  over  an  RS-232  line  to  inform  it  to 
start  a  particular  task  and  to  find  out  when  this  task  is  finished.  When 
the  IBM  system  wants  to  start  a  particular  task,  it  sends  the  desired  task 
number  (I  or  2)  to  the  Apple.  The  Apple,  upon  receiving  this  number, 
prompts  the  subject  for  his  name  and  instructs  the  subject  to  press  the 
joystick  button  when  he  is  ready.  The  Apple  then  returns  the  task  number 
to  the  IBM  to  indicate  the  start  of  the  task.  Once  the  task  is  finished, 
the  task  number  is  once  again  sent  to  the  IBM  to  indicate  this  fact. 

Interface  with  the  Oata  Acquisition  System 

The  IBM  communicates  with  the  MDAS  7000  data  acquisition  system  over  an 
IEEE-488  (GPIB)  bus  to  set  up  the  sampling  parameters,  to  inform  the  MDAS 
7000  to  sample  data  and  to  download  files  to  the  IBM. 

The  only  sampling  parameters  involved  are  the  sampling  rate,  the  sampling 
duration  and  the  channe1  selection.  The  MDAS  7000  is  instructed  by  the  IBM 
to  sample  two  channels  specified  for  the  frontalis  and  trapezius  muscle 
sites  for  one  second  at  a  rate  of  1000  samples  per  second. 
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Once,  at  the  beginning  of  the  experiment  and  at  exactly  2.75  and  at  3.25 
minutes  into  each  task,  the  MOAS  7000  is  instructed  to  sample  the  two 
channels  and  to  send  the  two  files  to  the  IBM  for  storage. 

File  Format  Convention 


The  data  files  consist  of  a  100  character  descriptive  header  followed  by 
1000  integers  in  ASCII  format.  The  descriptive  header  is  in  the  following 
format: 


Name/Data/T i me/T ask#  /W i ndow# /T rial # /Mu sc  1 e  Site 

The  blanks  are  filled  in  with  the  appropriate  values.  The  windows  are 
numbered  as  follows: 

0  -  At  rest 

1  -  2.75  minutes  into  the  task 

2  -  3.25  minutes  into  the  task 

The  name  of  each  file  begins  with  the  three  initials  of  the  subject,  where 
the  letter  "X"  is  used  if  the  subject  has  no  middle  name.  This  is  followed 
by  a  letter  representing  the  muscle  site:  'a'  for  frontalis,  ' c '  for 
trapezius.  If  the  file  represents  a  sample  at  rest  then  the  string 
'atrest'  is  appended  to  the  filename.  Otherwise,  the  task  number,  the 
window  number  and  the  trial  number  are  appended  to  the  filename  in  that 
orour  along  with  a  '.dat1  suffix. 
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Plot  the  Data  Files 


The  data  plotting  package  on  the  IBM  is  further  explained  in  section  3.5. 
Whenever  the  IBM  receives  a  set  of  two  data  files  from  the  MDAS  7000,  it 
plots  the  first  300  points  of  each  file  simultaneously  with  a  horizontal 
magnification  of  2  and  a  vertical  reduction  of  31.  This  remains  on  the 
screen  until  a  new  set  data  is  received,  or  a  new  task  begins. 

Interface  with  Experimenter 

When  the  program,  'OPERA',  is  executed  on  the  IBM,  the  screen  clears  and 
the  system  prompts  for  the  subject's  name.  Once  a  valid  name  with  three 
initials  is  typed,  the  subject  is  sampled  at  rest,  and  the  two  files  are 
plotted  and  stored  and  the  following  appears: 

Practice  first?  (y/n) 

Should  the  experimenter  wish  the  subject  to  go  through  a  practice  round 
before  any  data  is  taken,  the  experimenter  responds  with  ' y ' .  Otherwise, 
the  actual  experiment  will  begin. 

If,  however,  a  practice  round  is  desired,  the  following  prompt  appears: 
Practice  task  # 

Valid  responses  to  this  are  a  '1'  and  a  '2'.  Any  other  responses  will  be 
rejected.  The  task  number  is  then  send  to  the  Apple  He  which  then 
presents  the  appropriate  task  to  the  subject.  Once  the  subject  has  entered 
his  name  and  pressed  the  joystick  button,  the  task  number  is  sent  back  to 
the  IBM.  The  I8M  acknowledges  with: 

Task  #  started 
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where  the  blank  is  filled  in  with  the  appropriate  task  number.  The  IBM 
then  sits  in  a  waiting  state  until  the  task  number  is  once  again  received 
from  the  Apple  indicating  that  the  task  has  terminated.  Upon  receiving  the 
character,  the  IBM  types: 

Continue  practice?  (y/n) 

at  which  point  it  will  repeat  the  entire  practice  sequence  if  the  response 
is  affirmative.  Otherwise,  the  actual  experiment  begins  and  the  following 
appears: 

Trial  # 

The  experimenter  replies  with  either  a  ' 1 '  or  a  '2'.  The  system  responds 
with: 

Task  (1,2) 

The  experimenter  replies  with  either  a  '1'  or  a  '2'  depending  on  the  task 
difficulty  desired.  The  task  number  is  then  send  to  the  Apple  and  the  IBM 
types: 

Waiting  for  task  to  start  on  Apple  He... 

Upon  receiving  the  task  number  from  the  Apple,  indicating  the  start  of  the 
task,  a  timer  is  started  and  the  following  appears: 

Task  #  _ started 

where  the  blank  is  filled  in  with  the  appropriate  task  number. 


At  exactly  2.75  and  3.25  minutes  into  the  task,  the  data  acquisition  system 
receives  a  message  from  the  IBM  to  sample  the  two  channels  specified  for 
the  frontalis  and  trapezius  muscles  for  one  second  at  a  rate  of  1000 
samples  per  second.  The  acquired  data  is  immediately  sent  to  the  IBM. 
Once  received,  the  two  sets  of  data  are  plotted  simultaneously  for  the 
experimenter  to  check  the  signal  integrity. 

When  the  task  is  terminated,  the  Apple  indicates  this,  as  in  the  practice 
session,  by  sending  the  task  number  to  the  IBM.  The  IBM  acknowledges  this 
by  typing: 

Continue?  (y/n) 

If  the  response  is  affirmative  the  entire  process  after  the  practice 
session  is  repeated.  Otherwise  the  program  terminates. 
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APPENDIX  B 

INSTRUCTIONS  FOR  EXPERIMENTAL  SUBJECTS 
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This  experiment  is  part  of  a  program  of  continuing  research  at 
Perceptronics  in  human  (pilot)  performance  and  decision  making.  The 
purpose  of  this  particular  experiment  is  to  analyze  ways  in  which  human 
operators'  myoelectric  (muscle)  signals  respond  and  how  a  computer  might 
help  to  determine  an  operator's  mental  state  via  these  signals.  You  are  an 
integral  part  of  this  research  since  your  performance  provides  the  baseline 
data  for  predicting  operator  performance,  and  estimating  the  effectiveness 
of  computer-based  analysis  techniques. 

Tasks  Overview 


There  are  two  types  of  tasks  that  you  will  be  asked  to  perform  in  the 
experiment.  They  are:  (1)  motor  (control)  tracking  task,  and  (2)  motor 
plus  perceptual  (control  +  auditory)  task.  Each  task  will  be  approximately 
3-1/2  minutes  each.  Please  concentrate  on  the  task,  as  your  responses  and 
performance  will  be  closely  monitored  and  scored.  There  will  be  a 
predetermined  pay-scale  for  performance  levels,  thus  the  better  you  do,  the 
more  money  you  will  make.  At  the  end  of  each  session,  you  will  be  given  a 
questionnaire  to  fill  out.  Information  on  these  questionnaires  will  not  be 
used  to  rate  your  score,  so  please  use  your  unbiased  judgment  to  answer 
those  questions. 

The  following  paragraphs  describe  the  two  types  of  tasks. 

Motor  (Control)  Tracking  Task.  In  this  task,  you  will  see  an  object  called 
the  "target"  which  moves  freely  around  the  screen.  You  will  be  asked  to 
track  the  target  by  controlling  another  object,  called  the  "site"  with  a 
joystick.  The  objective  is  to  track  and  surround  the  target  with  the  site 
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as  the  target  continually  moves  around  the  screen.  Performance  will  be 
measured  as  the  average  distance  between  the  site  and  target  over  the 
entire  task. 

Motor  +  Perceptrual  (Control  +  Auditory)  Task.  This  task  incorporates  the 
first  motor  tracking  task  with  a  perceptual,  auditory  task.  In  addition  to 
the  first  motor  tracking  task,  you  will  be  required  to  respond  to  high 
tones  by  pressing  the  orange  button  on  the  top  of  the  joystick  whenever  a 
high  tone  is  heard.  There  will  be  low  tones  and  high  tones--you  must 
discriminate  between  the  two  tones  and  only  respond  to  the  high  tones. 
Overall  performance  will  be  measured  as  (1)  average  distance  between  the 
site  and  target  and  (2)  good  responses,  false  responses  and  missed 
responses  to  the  high  tones. 
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answer  on,  then  you  may  just  guess 

l  How  simple  was  the  task? 
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simple 
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2  What  level  of  effort  did  you  expend7 
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very 
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low 
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3  Hew  much  did  you  enjoy  the  session7 
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APPENDIX  D 

SAMPLE  PRINTOUT  AND  ANALYSIS 

WITH  PERFORMANCE  MEASURE  PRINTOUT 
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PERFORMANCE  MEASURE  PRINTOUT 


PERFORMANCE  MEASURE  PRINTOUT 


SUBJECT:  JAY  MARTIN 


LEVEL 

AVERAGE 

DISTANCE 

GOOD 

RESPONSE 

NO 

RESPONSE 

FALSE 

RESPONSE 

1 

21 

0 

3 

3 

1 

31 

2 

0 

0 

1 

20 

1 

0 

0 

2 

20 

2 

2 

2 

2 

15 

3 

1 

0 

2 

17 

1 

2 

2 

3 

17 

5 

2 

2 

3 

15 

7 

2 

2 

3 

19 

8 

0 

0 

4 

19 

9 

2 

2 

4 

22 

6 

0 

0 

4 

12 

11 

2 

2 

5 

19 

8 

0 

0 

5 

21 

7 

3 

2 

5 

21 

3 

0 

2 

Appendix  D  (Cont'd) 

AUTOBJ  SAMPLE  PRINTOUT  (SUBJECT:  AXC i 


MODELING  RESULTS  FuR  TIME  SERIES  AXCAATRE 


»■»*********+***«♦***•»«*♦♦+♦♦+*+♦+♦♦**+***♦*+*****+****++*.*+♦+++♦♦ 

DAT  A  :  2  =  AXCAATRE  10000.  1000  OBSERVATIONS 


DIFFERENCING  ON 


NONE 


BACKCAS~ING  :  OFF 

****♦*************+♦*****  **♦+****■***■*■■**•*+**  ****-***  *♦*♦*■*■ -*■■*•♦*♦•#■■*■■*•* 
JN  I  VAR  I  ATE  '"ODE:-  PARAMETERS 

**■»*■*■**■«•**♦********+*****+*********+******+*********************  + 

-ACTOR  LAG  COEFFICIENT  T  RATIO 

**♦•**♦♦****♦********++**************  ******♦****************♦.**+*» 


MEAN 

AUTOREGRESSIVE  1 
AUTOREGRESSIVE  1 
AUTOREGRESS  I  VE  £ 

autoregressive  £ 

MOVING  AVERAGE  1 
MOVING  AVERAGE  2 


1 2557E+01 
17241E+00 
£01 23E+00 
47839E+00 
£27 1 6E+00 
43090E+00 
14074E+00 


4.  94 
6.  £0 
16.  06 
7.  99 
14.  64 
-4.  3£ 

*  *************************** 


T-E  RESIDUAL  STATISTICS 


iu'JHrttb 


MEAN  SQUARE 
R  SQUARED 


. 4 1 064E-0 1 
. 42796E-04 
. 51239E+00 


DEGREES  OF  -REEDuM 
NUMBER  OF  RcSIDUAlS 


THE  PI  WEIGHTS 


.'MVARIhTE  BuX-JEnkInS  MODEuS  CAN  BE  EXPRESSED  AS  a  WEIGHTED  SuM  OF  t_=;  zc 
;--S  A  RANDOM  SHOCK.  DEFINED  AS  T-E  a  I  WEIGHTS,  MODEuS  EXPRESSED  IN  T-ESS 
'ERMS  ARE  USEFUL  FOR  PRECASTING  AnD/OR  COMPARISON  TO  OTHER  TYPES  0-  MUDS 


G=.n-RA'_  -u4iM  Or  TrE  mGDEl  EXPRESSED  BY  T’-E  PI  wEIGH“S  IS  : 


lons~  ant  +  Pi(l;  *  Z(t-1>  +  Pi(£)  *  Z<t-£)  ...  +  Pi(r,)  *  Z(t- 


r? 

>  , 


■  - 

1 " 

■  . 

» 

U  » 

1  & 


n 

i*-. 

s; 


4 


v". 

u. 

< 


1 

LAG 

1  T4E  PI  Lt I G- i 

Constant 

— . 386402* Ol 

Pi  (  i ) 

4609uE*00 

Pi  (  3) 

— .  40467E-*-oo 

Pi  <  3) 

39584£*00 

Pi  (  4 ) 

190362*00 

Pi  (  5) 

— . 9 i 543E— 0 1 

Pi  <  6) 

44023E-01 

Pi  (  7) 

21 170E-01 

Pi  <  a ) 

lOiaiE-o: 

Pi<  10) 

46075E*00 

Pi (  11) 

231  192*00 

Pi (  12) 

2622 iE—0 1 

P i (  13) 

. S2S97E-01 

Pi<  14) 

. 39769E-01 

Pi<  15) 

. 19125E-01 

Pi (  30) 

— • 237 1 6E*00 

Pi  (  31) 

1 1405E*00 

Pi (  32) 

13720E-01 

Pi (  33) 

. 41 127E-01 

Pi <  34) 

. 19778E-01 

Pi<  70) 

.  140742*00 

Pi (  71) 

. 67684E-01 

Pi (  72) 

. 569562-01 

Pi<  73) 

. 55712E-01 

Pi (  74) 

. 26792E-01 

Pi<  75) 

.  12864E-01 

Pi<  ao> 

. 67863E— 0 1 

pi <  ai) 

. 32539E-01 

Pi <  83) 

116392-0’ 

Pi < 100) 

. 33380E-01 

Pi  (  101 ) 

. 16052E-01 

Pi (140) 

19809E-01 

_  _ _ — - _ * _ _ _ _ _ 

mCCE-InG  “I£S.l*5  -OR  T I  "E  SERIES  AXCCATRE 


****************  *****4*****+-*  *■*-**-*■**■*■■*•**-***+*-»•■*■*••*-*■***■*■*-«■  ♦■*•♦-*•-*•*-*■♦-•• 

DATA  :  Z  =  AXCCATRE  10000.  1000  OBSERVATIONS 

D I - rERENC I NG  ON  Z  :  NONE 
BACKCAS~ I NG  :  OFF 
UNIVARIATE  MODEL  PARAMETERS 

************** *4*********4***4*44*4****44+***44*****4*4*4+ 444**44 

FACTOR  LAG  COEFFICIENT  T  *ATIQ 


1 


*********************************  *********  ***♦*♦ ♦*♦***  *  *4*44**44  • 


1 

mEAn 

2* 

4- 

AUTOREGRESSIVE 

1 

1 

3 

AUTOREGRESSIVE 

1 

£ 

4 

AUTOREGRESSIVE 

1 

5 

AUTOREGRESSIVE 

£ 

30 

£ 

MOVING  AVERAGE 

1 

2 

7 

MOVING  AVERAGE 

Zt 

50 

*********4**************4* 


18903E+01 

.  l7£56E-*-00  5.34 

. 561 37E+00  .0. 83 

.  17522E-*-00  4. 24 

.  1 i SaaE-00  3.  37 

.275102*00  4.40 

l744<*E+00  -5.  34 


TnE  RtS-DOAL  STAiIbTIuS 


SL"1  OF  SQUARES 
"SAn  SQUARE 
R  SQUARED 


12556E+00 
1308 IE— 03 
622Q7E+00 


DESRuiS  OF  >-  *<i“L 
NUMBER  OF  R“3IDu 


THE  PI  UtlGnTS 


NIVARIATE  BOX— JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  wS 
uuS  A  RANDOM  Sr»OCK.  DEFINED  AS  T«E  PI  WEIGHTS,  (MQDEu 
ERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  C"-ER 


THE  GENERAL  FORM  OF  Tr-E  MODtu  iXPRiSSSD  BY  T>iE  PI  wEI 
Z(t>  =  Constant  -*■  Pi(l)  *  Z(t-l)  Pi  (2)  *  Z  (t-2;  ...  + 


1  • 

uAG 

Tt£  PI  Hr  I  G~T  : 

Constant 

*“«  4£6ci0t^’0 1 

!  Pi  ( 

1) 

17256E+00  i 

:  Pi  < 

2) 

— . 83646E+00  i 

!  Pi  < 

3) 

-.  £££6'5E-K>0  ! 

i  Pi  < 

4) 

— . 2301 1E+00  ! 

!  Pi  < 

5) 

—  .6 1 26  i  E— 0 1  1 

l  Pi  ( 

6) 

— . 63302E— 01  ! 

I  Pi  ( 

7) 

16853E-01  ! 

;  Pi  ( 

8) 

-. 174146—01  : 

!  Pi  < 

30) 

H5Q2E+00  ! 

!  Pi  ( 

31  ) 

. 1 9986E— 01  1 

!  Pi  ( 

32) 

. 33157E-01  ! 

!  Pi  ( 

uo  ) 

. 25793E-01  ! 

:  Pi  ( 

50) 

.  i744^E-»-00 

:  Pi  ( 

51  ) 

. 30101E-01 

;  Pi  ( 

52) 

.  145S2E-K>0 

1  Pi  < 

53) 

. 38Q47E— 0 1 

i  Pi  ( 

54) 

. 40141E-01  : 

:  Pi  < 

55) 

. 1 0687E— 0 1  ! 

!  Pi  ( 

56) 

.  1  i042E-01  : 

:  Pi  ( 

80) 

. 202056—0 1  : 

!  Pi ( 100) 

-. 3043.E-01 

:  p i ( i 02 ) 

-.  S5454E-0 1 

in  »h 


y* C 3 n  —  I N *3  ^ b. S  U  L.  3  ^GR  *•" I M h.  SERIES  PXC^lli* 


**  +  ***»  +  ***+****#*+*4'****4f*****  +  *  +  '»*+4+*'l»,  +  **^***  +  +**'*  +  *  +  *^  +  *^  +  ' 

DAT P  :  Z  =  PXCP111.  1000.  1000  CBSERV A~ 1 1  ’ 


D I - - i"<£.NC ING  On 


Z  :  NONE 


e-*****-*++*  +  +  +  *  +  *i*-  +  -«-**  +  - 


3ACKCAS” I NG  :  QR r 

*  **************  ************  *i 

UNIVARIATE  ^ODEi_  PARAME“ERS 

FACTOR  LAG  COEF-ICIENT  T  RATIO 


♦  *■»•*■****•*•»****•**■»***■***■*•»**■»■******■»***■»■**•»•**+*••■*********■*■*  *•****♦♦ 


m- AN 

-.  27796E+01 

ALTCREG 

RtiSSIV-  1 

20 

.  1 1512= -*-00 

MOVING 

AVERAGE  1 

1 

.  340E5E+00 

MOV  I NG 

AVERAGE  1 

3 

9 1 043E— 0 1 

MOV  ING 

AVERAGE  E 

40 

1 3533E-00 

^  i  •  — 1_ 

-3.  00 
-4.  El 

*****++*+******+*****^********.**  **♦•*•*  ♦•if'*********-*-*'-*’****-***-*- •*-*■■»•*■' 


T fit  RtSIDuAu  5T A  i  ISTIU3 


SUM  O'-  SQUARES  :  . 279 1 1 t^O 1  DEGREES  OF  F  Rt.- DOT  : 
m£AN  SQUARE  :  . 286262-02  NUMBER  0-  RSSIDUA-S  : 
R  SQUARED  :  .  12517E+00 

T~E  PI  wEIG— T3 


■jNIVARIA”  BOX-jENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  4E I  G— 72 
P*_JS  a  random  SmCCK.  DEFINED  AS  T-E  pi  UEIGATS,  MCDEi_5  EaP 
-R^S  PR-  .Ss'l'l  -OR  RCRECASTING  AnD/OR  COMPARISON  Ti j  Q~  — E-r  TVCji  r-  « 


-  GEnEFAl.  -CRM  OF  T-t  MCD-u  £XP>-<tSS2D  BY  T-E  PI  wEIG--3  IS  : 
Constant  *  Pi(i)  *  Z(t-l)  +  Pi  <2>  *  Z  <t-2)  ...  +  o !  ;  r,  ■  *  z 


!  LAG 

THE  PI  WEIG-” 

!  Constant 

—  . 

362S5E+0 1 

Pi  ( 

i  ) 

-. 

34025E—00 

Pi  ( 

2) 

1 1577E+00 

Pi  ( 

3) 

5 l 652E-0 1 

Pi  ( 

4) 

435522—0 1 

Pi  ( 

5) 

27060E-0 1 

Pi  ( 

20) 

1  121  2- -*-00 

Pi  < 

2.  ) 

33 1 46E-0 1 

Pi  ( 

22) 

-. 

1 29 Q02-0 1 

Pi  ( 

40  ) 

1 3533E-I-00 

Pi  ( 

41  ) 

460^*62— ‘j  1 

Pi  < 

42) 

1  US&7E— 1 

Pi  < 

60 ) 

151736-01 

P  i  ( 

80 ) 

*  • 

133 1 «2-  1 1 

it  ir 


DATA  :  Z  =  AXCC 1 1 i .  10000. 

DIFFERENCING  ON  Z  :  NONE 
BACKCASTING  :  OFF 


h*****»+*+++*********4****+i 

1000  observh-: 


UNIVARIATE  MODEL  PARAMETERS 

#•»■»***•»*■*■»*•**■»***■»*-»*■************■***■»*-•*****•»**•»•♦***■<•■*•**■**■****■♦*  + 
FACTOR  LAG  COEFFICIENT  T  RATIO 


1  MEAN 

a  AUTOREGRESSIVE  1 

3  AUTOREGRESSIVE  1 

4  AUTOREGRESSIVE  1 

5  MOVING  AVERAGE  1 


*♦♦+♦+■♦■♦**•**♦**♦***+*+*♦*♦»+♦+*♦♦**-»**-»*♦♦ 

19609E+01 

1  . 27079E+00  a. &a 

a  .  34385E+00  1 1 . £7 

3  .  1 9783E-*-00  6.35 

£0  -.  l68l£E+00  -5.61 

****-»*■»■■»■••****•*+**+***♦**•*■*+-*■•■*■»•♦-»■+*  +  -»■>■»■■•■■*•* 


**-»-**  •*•*-»♦-»  •***♦*♦  ■*-«■-*•  **■*•*  *•*•**-*  ■**■*-*•*•**  *■**■*■*•#*■**■#•***■*■■*•■*♦  **♦-«■+** 

Tr.E  RESIDUAL  STATISTICS 


SU'T  Or-  SQUARES 
*EPN  SQUARE 
R  SQUARED 


. Q7317E-02 
. aaoais-05 
. 56219E+00 


DEGREES  OF  FREEDOM 
NUMBER  OF  RESIDUALS 


THE  PI  WEIGHTS 


nj  I  vAR  I A-E  BOX-JENKINS  MGDeLS  CAN  Be  EXPRESSED  AS  A  wEIG-~ED 

_ S  A  RANDOM  SHOCK.  DEFINED  AS  T~E  PI  WEIG^”S,  MuDEeS  EaPRE 

ER^S  ARE  USEFUL  FOR  FORECASTING  A-mD/OR  COMPARISON  tq  GtmE-<  Tv  PEE 


Tee.  GeNeRAL  FORM  qf  ThE  MGDeL  eXPReeSED  BY  ~  — E  3  I  wEIGe  S  IS 
c<  t)  =  Constant  *  Pi(l)  »  Z(t-l)  Pi(£)  *  Z(t— £)  ...  Pi  in)  * 


in  in 


;« * 


- 

AG 

-E  PI  we  I  ip-  1 

Con 

st  ant 

-.  259  1  1 

pi 

( 

l ) 

i^7<.)7'cE-*-U0 

pi 

( 

2) 

— . 343S5E-00 

pi 

( 

3) 

-.  l97QSE-f00 

pi 

( 

20) 

. 1 38 1 22+00 

pi 

( 

2 1 ) 

. 50940E-01 

pi 

< 

22) 

. £46332-0 1 

Pi 

( 

23) 

.  372 1  **E— 0 1 

Pi 

( 

40) 

-.  353SSE-0 1 

P  . 

£ 

42) 

-.  1 2 1  66E— 0 1 

MGD£i_InG  RESULTS  FOR  TIME  SERIES  AXCA121. 


DATA  :  Z  =  AXCA121 .  1000. 

DIFFERENCING  ON  Z  :  NONE 
BACKCASTING  :  OFF 
UNIVARIATE  MODEL  parameters 

FACTOR  LAG  COEFFICIENT 


1000  OBSERVA-:  0!\  ; 


T  RATIO 


B  p 


1  MEAN  28824E+01 

2  AUTOREGRESSIVE  1  10  .  20174E-*00  £.  £7 

3  MOVING  AVERAGE  1  1  .  E7375E-t-00  3.65 

4  MOVING  AVERAGE  1  2  -. 12206E+00  -3.35 

5  MOVING  AVERAGE  2  <+0  -.  39747E-0.  -2.72 


28824E+01 
10  . 20174E-+-00 

1  .  27375E+00 

2  -. 1 2206E+00 

40  -.  S9747E-0. 


T^E  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


. 1 465SE+0 1 
. 1 4880E-02 
.111 39E+00 


DEGREES  OF  FREtD'JM 
NUMBER  0^  R^SIDuA-S 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JE.NKINS  MODELS  CAN  Bt  EXPRESSED  AS  A  WEIG-iTED 
PL.JS  A  RANDOM  Shock.  defined  hS  T-E  PI  WEIGHTS,  MODELS  EaPR: 
TE  RMS  ORE  USEFUL.  FOR  FQR£Cpsjjfg0  A:\D  /  OR  COM  PAR  ISON  TO  C^HEk 


t-E  GENERAL  FQ^tvi  qp  r-E  mQDEl  EXPRESSED  BV  ~-E  p  I 


Z  ( t  ) 


Constant  +•  Pi(l)  *  Z<t  — 1)  +  Pi(2u  *  Z  i 


1  l  n 


♦  L 


m  w  h 


WF  I  G  —  i 


Const  ant 

Pi 

( 

i ) 

Pi 

( 

2) 

Pi 

( 

3) 

Pi 

( 

10) 

Pi 

( 

1 1 ) 

Pi 

< 

40) 

Pi 

< 

41 ) 

Pi 

< 

50) 

b747(.)b*U  1 
2737rrE-*-oO 
47121 E-Ol 
463 1 3E— 0 1 
20175E-*-00 


89747E— 0 1 
2*56aE-01 
18105E— 01 


MODELING  RESULTS  FOR  7I*G  SERIES  AXCC121 


DATA  :  Z  =  AXCC121.  10000. 


iooo  cbs ~ rva~ i c 


DI>-f ERfcNClNG  ON 


NONE 


BACKC AST  I NG  :  OkF 

********■»+•»♦•*****-**■»*■*  ■**•»*■*-*■»■*■*•*+■*■***■*■+■»•■**++*+*■»•■*•**■•*■»■*-«•+-»■♦•*■♦ -n 

UNIVARIATE  MODEL  PARAMETERS 

FACTOR  i_AG  COEFFICIENT  T  RATIO 


MEAN 

AUTOREGRESSIVE  1 
AUTOREGRESSIVE  1 
MOVING  AVERAGE  1 
MOVING  AVERAGE  1 
MOVING  AVERAGE  2 


194S1E+01 
65 102E+00 
20266E+00 
23 592E+00 
523 1 i E+00 
2004lE-«-00 


**+♦-*-**■»*•»■»♦■♦*♦♦-*♦•*■»*♦*+*■*•*■••**■•■*+♦*♦***♦*♦■» 


9.  57 
7.  40 
7.  46 
6.  72 
-6.  15 

*+*+*+t+H 


THt  RESIDUAL  S i A  1  I S 1  I CS 


Sum  OF  SQUARES 
^EAn  SQUARE 
R  SQUARED 


. 75131 E-01 
. 75813E-04 
.  1 1171  E-*-00 


DEGREES  0“ 
NUMBER  OF 


P  R-  -D0|Y1 
<E5  _  D:_'A_S 


THE  PI  WEIGHTS 


-N  I  VAR  I  ATE  BOX -JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  wEIGr-TED  5-M 
PLuS  A  RANDOM  SHOCK.  DEFINED  AS  t~E  PI  WEIGHTS,  MljE_5  E*PRE5Si 
TERMS  ARE  USEFUL  FOR  FORECASTING  AnD/uR  COMPARISON  TO  G-,-,E'<  Tva 


T«t  G'tN^RAl.  rQRw  OF  T  — £  MGDti_  tjARuSSED 


1  a  ^  n  .  — 


Z  ( t  > 


Constant  ai(l)  *  Z(t-l) 


p  l  (  2  >  *  Z  '  t  -  _ 


4*  wJ 


♦  z  k- 


n  ui 


Const  arit 

Pi 

< 

l ) 

Pi  ( 

2) 

Pi 

( 

3 ) 

Pi  < 

4) 

Pi 

( 

5) 

Pi 

( 

6) 

Pi 

( 

7) 

Pi 

( 

8) 

Pi 

( 

3) 

Pi 

< 

10) 

Pi 

( 

11) 

Pi 

( 

12) 

Pi 

( 

13) 

Pi 

( 

14) 

Pi 

( 

15) 

Pi  < 

16) 

Pi  < 

17) 

Pi 

< 

18) 

Pi 

( 

19) 

Pi 

( 

20) 

Pi 

< 

22) 

Pi  < 

23) 

Pi 

< 

24) 

Pi 

< 

25) 

Pi 

( 

26) 

Pi 

< 

27) 

Pi 

( 

£8) 

Pi  ( 

29) 

Pi 

( 

30) 

Pi 

( 

31  ) 

Pi 

( 

32) 

Pi 

( 

33 ) 

Pi 

< 

34) 

Pi 

< 

35) 

Pi 

( 

36) 

Pi 

< 

37) 

Pi 

< 

38 ) 

Pi 

< 

40) 

Pi 

( 

42) 

Pi 

( 

43) 

Pi 

( 

44) 

Pi 

( 

45) 

Pi 

< 

46) 

Pi 

< 

47) 

Pi 

( 

46) 

l£484c-*-08 
£3592E-K>0 
— .  1  8298E+0  1 
— .  6i62iEH"00 
— . 78S70E+00 
508 422 -*-00 
532522+00 
39159E+00 
37095c -*-00 
£9£36E-*-00 
— .  ££30£E-*-00 
2J.499E+00 
1883 IE +00 
— .  156 89c +00 
— . 13552E+00 
1  i  40+E+00 
97797E-01 
82729E-01 
— . 70676E— 01 
— . 5995 i E-0 l 
. 149302+00 
.  209462+00 
. 920588-01 
.131 30E+00 
. 79132E-01 
. 87354E-01 
. 68003E-0 1 
. 6032+E-01 
. 46666E-01 
. 42565E-01 
. 34454E— 01 
. 30395E-01 
. 25194E-01 
. 2i84hc-01 
. 18332E— 01 
• l 575cc  Ol 
. 1 3306c— 0 1 
. 1 1379E-01 
-.31 936c— 01 
4344iE-0l 
-. 19691c— 01 
27370E-01 
1 67582—0 1 
1827 iE— 01 
-. 1 30772—01 
186432-01 


MODELING  RESULTS  FOR  TIME  SERIES  AXCA211 


**•*■**•*■****•**-»**■»'*■»*■*■»  ^ **■*•*■■#•■** *******♦******♦***•*•♦**■****+♦+♦+*■♦♦♦♦♦ 

DATA  :  Z  =  AXCA21I.  1000.  1000  OBSERVA" jN  = 

DIFFERENCING  ON  Z  1  NONE 


BACKCASTING  :  OFF 


UNIVARIATE  MODEL  PARAMETERS 


FACTOR  LAG  COEFFICIENT 


T  RA i 10 


♦♦■a-************-*****-* -fr******-***-*-*  ■*•****•*  ********+♦*-*■-*■■*•-*•-»•■»•■*•■*•-*■-*■-*■-*  ♦-*•->■ 


1 

MEAN 

33753E-*-0l 

o 

La 

AUTOREGRESSIVE 

1 

1 

— .  54S06E-*"00 

17.  E& 

AUTOREGRESSIVE 

1 

2 

-.  24740E+00 

-7.  86 

4 

AUTOREGRESSIVE 

2 

10 

.  14733E-«-00 

4.  £0 

5 

AUTOREGRESSIVE 

2 

20 

.  16534E+00 

5.  27 

6 

MOVING  AVERAGE 

1 

30 

-.  24072E-**00 

-7.  34 

THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
mEPn  SQUARE 
R  SQUARED 


.  3273SE-K51 
. 33681E-02 
. 345S2E+00 


DEGREES  0-  FREEDOM  : 
NUMBER  OF  RE5ID!_'A_S  : 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SwM  0-  ~-E 
aLuS  A  RANDOM  SHOCK.  DEFINED  AS  T-E  PI  WEIGHTS,  MCDEuS  EXPRESSED  In  ~-i 
TERMS  ARE  USEFUL  FOR  FORECASTING  An  D/OR  COMPARISON  TO  0~HE  ■<  TYPES  Q-  mi 


T«->E  GENERAL  FORM  OF  THE  MODti_  EXPRESSED  BY  T|-'E  PI  wEIGk'S  lb  : 

t - 2 )  ...  Pi  (•'•■>  *  Z 


Z  (t ) 


Constant  Pi<l)  *  Z<t  — 1>  Pi(2/  *  Z< 


s_AG  I  i  "it  s  I  w-.IG-iT 


Constant 

-.  20677E+0  1 

Pi  ( 

1) 

.  5*606E-*-00 

Pi  < 

2) 

. 247402+00 

Pi  < 

10) 

14738E+00 

Pi  < 

ir> 

804782—0 1 

Pi  < 

12) 

— . 3646 i E— 0 1 

Pi  < 

20) 

-.  1 6534S-r-00 

Pi  < 

2i> 

S0286E-01 

Pi  < 

22) 

-. 40S05E-01 

Pi  ( 

30) 

.  24072E+00 

Pi  ( 

31) 

13145E+00 

Pi  ( 

32) 

-. 59553E-01 

Pi  < 

40) 

. 35477E— 0 1 

Pi  ( 

41) 

. 19373E-01 

Pi  < 

50) 

. 398012-01 

Pi  < 

51) 

.21  73+E— 0 1 

Pi  ( 

60) 

-. 57947E-01 

Pi  ( 

61) 

. 3i642E-0l 

Pi  ( 

£2) 

. 1 43362—01 

Pi  ( 

SO) 

.  13949E-01 

MGD£LING  RESULTS  FOR  TIME  SERIES  AXCC211. 


* -**-**■* -*-***■*  •»■**■*  •*■***-*•*-*■+*•**•#  ■*■*■*■*■**•**■#*-*  •*•**+■**■*•*■*  **■ ■*•*■***•***■*■+  ♦  ♦ 

DATA  :  Z  ■  AXCC211.  lOOOO.  1000  uBSERVA"M'NS 

DIFFERENCING  ON  Z  :  NONE 


BACKCAS'ING  :  OFF 

*****•***■»**•****■»•»■»■*■»•***■»■*•»•*•#*■*■»■*#■*#****■»■******■****■••♦**♦+*+*•*♦■•■■*♦♦ 
UNIVARIATE  MODEL  PARAMETERS 

**#**•»*■*  •***■»*•»*****■»•»■»•*-»■»•*-»  •**■*•*•»***•»******•»■  ■*•»■»*****■«■•»■*•***♦+■•■+♦■•■+■*• 
FACTOR  LAG  COEFFICIENT  T  RATIO 


*#■*#•»**■»•****■»**•*■»■»*-*•*■»■»■»*■»**■**-*******■»*-»■»******+■*■»■•*■»•*+■«■+•«•  +  ■*•*♦■*••*•*♦ 


1 

O 

w 

MEAN 

AUTOREGRESSIVE 

1 

2 

-. 15278E+01 
. £91182+00 

3 

AUTOREGRESSIVE 

1 

3 

. 25295E+00 

4 

AUTOREGRESSIVE 

2 

10 

. 32819E+00 

5 

MOVING  AVERAGE 

1 

1 

. 2.918E+00 

6 

MOVING  AVERAGE 

1 

2 

. 62209E+00 

7 

MOVING  AVERAGE 

2 

60 

.  1  4956E+00 

17.  s: 

9.  24 
10.  £9 
6.  89 

13.  34 
a  =;  3 

.  AW- 


**•*****•*•»****•*•*****•*■*■*******■*■**■****■****■*■•■■»****■**  +  ■••  +  ♦■••♦♦*■*•♦*♦♦♦♦»♦ 


TiE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


47244E-01 

48208E-04 

22098E+00 


DEGREES  Or-  •- RsEDuM  : 
NUMBER  OF  R2SID\-A_S  : 


T HE  P  I  Hi  I  Gn  l  S 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  Wt-IGnTED 
=LL-S  A  RANDOM  ShOCK.  DEFINED  AS  ThE  hi  WEIGHTS.  '"uDEuS  E^HRi 
'ERMS  ARE  USEFUL  "OR  FORECASTING  AnD/OR  COMPARISON  TO  C'“hFR 


THE  GENERAL  FORM  OF  THE  MODEL.  EXPRESSED  BY  T-E  A I  wEIG 
Z(t)  =  Constant  -*-  Pi(l)  *  Z(t-1>  -*•  Pi  (2)  *  Z  <t-2)  ...  3 


LAG  !  THE  PI  WEIGhT 


Constant 

-.  2221  IE *-02 

Pi  ( 

l ) 

2i91BE+00 

Pi  < 

£) 

— . 13613E+01 

Pi  ( 

3) 

— .  687S6E-I-00 

Pi  ( 

4) 

99758E-f00 

Pi  < 

5) 

6464W-E+00 

Pi  < 

"6) 

-. 76227E+00 

Pi  < 

7) 

-.  569E2E-I-00 

Pi  ( 

8) 

-. 59S97E+00 

Pi  ( 

9) 

-. 48539E+00 

Pi  < 

10) 

80710E+00 

Pi  ( 

11) 

-.  478S5E+00 

Pi  ( 

IE) 

-.  3S0£7E-*-00 

Pi  ( 

13) 

-. 298E5E+00 

Pi  < 

14) 

30193E+00 

Pi  ( 

15) 

£517  j.E-K'0 

Pi  < 

16) 

E4300E+00 

Pi  < 

17) 

-. 20985E+00 

Pi  < 

18) 

— .  1 97 1 6E-*-00 

Pi  ( 

19) 

-.  17376E+00 

Pi  ( 

SO) 

-.  16074E-*-00 

Pi  ( 

Si) 

-. 1433EE+00 

Pi  < 

SS) 

-.  13141  £-*-00 

Pi  < 

23) 

— .  1  l79SE-*-00 

Pi  ( 

24) 

— .  10760E+00 

Pi  < 

25) 

-. 96966E- 0 1 

Pi  < 

SS) 

-.  8819EE-01 

Pi  < 

27) 

-. 79653E-01 

Pi  ( 

S3) 

- . 723EEE— 01 

Pi  < 

29) 

-. 65403E-01 

Pi  ( 

30) 

— . 59326E— 01 

Pi  ( 

3 1 ) 

53690E-01 

Pi  < 

32) 

48674E-01 

Pi  ( 

33) 

-. 4h063E-01 

Pi  ( 

34) 

39938E— 01 

Pi  ( 

35) 

-. 36166E-01 

Pi  ( 

) 

— . 3E773E— 0 1 

Pi  ( 

37) 

-. E9S83E— 0 1 

Pi  ( 

3Q ) 

— . 2S893E— 0 1 

Pi  ( 

39) 

E**3£0E-0l 

Pi  ( 

40) 

-. EE069E— 0 1 

-YP-S 

H  ~S  - 

i  ( n ; 


in  rn 


Pi (  4i> 
P i <  42) 
Pi<  43) 
Pit  44) 
Pi<  45) 
Pi (  46) 
P  l  (  47J 
P i <  48) 
Pi (  60) 
Pit  61) 
Pit  66) 
Pit  63) 
Pi (  64) 
Pit  65) 
Pi (  66) 
Pit  67) 
Pit  68) 
Pit  69) 
Pit  70) 
Pit  71) 
Pit  72) 
Pit  73) 
Pit  74) 
Pit  75) 
Pit  76) 
Pit  77) 
Pit  78) 
Pit  79) 
Pit  80) 
Pit  81) 
Pit  82) 
Pit  83) 
Pit  84) 
Pit  85) 
Pit  86) 
Pit  87) 
Pit  88) 
P i < 1 20  > 
Pi <122) 
Pi t 123) 
Pi ( 124) 
Pi  (125) 
Pi  (126) 
Pi (127) 
Pi  ( 126) 
Pi t 129) 
Pi ( 130) 
Pi ( 131 ) 


1999. E— 01 
-.1811  IE— 01 
16406E— 01 
1 48662— 0 1 
13464E— 01 
-. 12197E— 01 
1104SE-01 
-. 10009E— 01 
— .  15262E+00 
-. 35550E-01 

2061  it -*-00 

-. 1051 26+00 
15126E+00 
98549E— 01 
1  i570E+00 
— . 66665E— 0 1 
90970E— 01 
— . 73852E— 01 
12185E+00 
-. 72650E-01 
-. 57808E-01 
-. 45453E-0 1 
45924E-01 
-. 38341E-01 
36973E— 0 1 
- .319565-01 
30O04E-01 
— . 2S456E— 0 1 
24464E-0 1 
-.218205-01 

2000 iE-01 
-. 1 7958E— 01 
16379E— 01 
— . 1 476 1 E— 01 
13424E-01 
-. 121255-01 
—  •  liOO'aE— 01 
-. 22835E-01 
30833E— 01 
-. 157292-01 
-. 22629E-01 
1 47442—0 1 
17309E-01 
-. 1 2966E— 01 
13610E-01 
-. 1 10492-01 
— .  16227E— 01 
-. 1 0869E— O 1 


mODEi_ ING  RESULTS  FOR  TIME  SERIES  AXCA22 1. 


********■»**♦■»******■»*********•»■»**+*•♦*****+*++***+»■»♦♦**+*♦♦■»*■*+*■>' 


DATA  :  Z  =  AXCA221.  1000. 


1000  GBSERVA  :  I  u  \G 


DIFFERENCING  ON 


BACKCASTING  :  OFF 


NONE 


UNIVARIATE  MODEL  PARAMETERS 


FACTOR  LAG  COEFFICIENT 


i  *<ATIO 


1  MEAN 

2  AUTOREGRESSIVE  1  1  3C622E+00  -10.06 

3  AUTOREGRESSIVE  1  3  .  13107E+00  4. 76 

4  AUTOREGRESSIVE  1  4  .  13222E+00  4.  15 

5  AUTOREGRESSIVE  2  10  . 34040E+00  .0.75 

6  AUTOREGRESS I VE  2  20  . 207Q9E+00  6.  36 

7  MOVING  AVERAGE  1  40  -.167296+00  -4.66 

***-*■***•*•*-*■*•**■*■*-*  *********+*•***•#*+*-*■+■**■*-*••*•*•*•*+-*■•*••»•-*■-*♦■*+++  +  +  +  ♦■*+ 

TmE  RESIDUAL  STATISTICS 


-. 33S2SE+0 1 
1  -. 3C622E+00 

3  . 15107E+00 

4  . 1 3222E+00 


3uM  uF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


. 23 1 6 1E+01 
. 300S4E-02 
. 35080E+00 


DEGREES  u1"  fRi-tDC11” 
NUMBER  uR  Rr 5  I DuAlS 


THE  PI  WEIGHTS 


nIVQRIA  s.  BOX— JENKINS  MODELS  CAN  Bfc.  EXPRESSED  AS  A  wEIG-T-D  5_l',,  0  - 
uuS  A  RANDOM  ShOCK.  DEFINED  AS  TmE  PI  HEIGHTS,  MODEuS  EXPRESSED  In 
ER"*S  ARE  USEFUL  FOR  FORECASTING  AnD/OR  COMPARISON  TO  C“hE-<  T*  P£=  C- 


Tht  GENtRAL  RQRM  Oh  ThE  MODfcL.  EXPR^SSn-D  BY  ThE  u  I  wt.'j"  3  .  b 
*  Constant  +  Pi(l)  *  Z(t-l)  +  Pi(£)  *  Z(t-£)  ...  +  R  l  ( n  >  * 


ill  ill 


T'-it  PI  wtlG^T 


Const  ant 

Pi 

< 

1) 

Pi 

( 

3) 

Pi 

( 

4) 

Pi 

< 

10) 

Pi 

< 

1  1) 

Pi 

< 

13) 

Pi 

( 

14) 

Pi 

< 

20) 

Pi 

( 

£i> 

Pi 

( 

23) 

Pi 

< 

£4) 

Pi 

( 

40) 

Pi 

( 

41) 

Pi 

( 

43) 

Pi 

( 

44) 

Pi 

< 

50) 

Pi 

( 

51) 

Pi 

( 

60) 

Pi 

< 

6i> 

Pi 

< 

80) 

44675E+0! 

.  30822E+0 
15107E+00 
13222E+0 
34040E+00 
1 0498E+00 
. 5 142SE— 0 1 
. 45007E-01 
20789E-‘-00 
64076E-01 
. 31407E-01 
. 274S7E— 0 1 
.  167£9E-t"00 
— .  5i  56 i E— 0 1 
. 25273E— 0 1 
. £81 1SE-01 
. 56946E-0 I 
. 17552E— 01 
. 34778E— 01 
.  10719E-01 
27986E-01 


MODELING  RfcSULTS  FOR  7 1  Ms.  SERIsS  AXCC221. 


DATA  :  Z  =  AXCC221.  10000.  1000  Q3SERVA7 M nS 

DIFFERENCING  ON  Z  :  NONE 
BACKCASTING  :  OFF 

* ****  *********************  *******************  *******it**-*"^***^- -*■-*■ -^  ■*■*  *■  -+ 

UNIVARIATE  MODEL  PARAMETERS 

FACTOR  LAG  COEFFICIEN~  T  ^A“Ij 


1 

wtH(N 

14343c.  *01 

2 

AUTOREGRESSI VE 

1 

10 

.  £334SE-<-00 

a.  96 

3 

moving 

AVERAGE 

1 

1 

.  3S625E+00 

_  u .  6  o 

4 

MOVING 

AVERAGE 

1 

3 

-.21 290E-*-00 

—  7.  32 

5 

MOVING 

AVERAGE 

2 

20 

15848E+00 

-4.  77 

#*-»****  +  ***********-*+**-)HH(-  +  *-*********  +  **-»  +  *-»-**  +  *»»»-nn 

the  residual  s~atis~ic3 


SUM  Of-  SQUARES 
MEAN  SQUARE 
R  SQUARED 


.  13406E+00 
. 136102-03 
.  £27  1  SE-i-00 


DEGREES  QR  r REuD 
NUMBER  OF  R-SIDu* 


THE  PI  WEIGHTS 


_iNl  VARIATE  BOX— JENK I  NS  MODELS  CAN  BE  EXPRESSED  AS  A  wEIGhTxJ  5i~M 
3L'_'S  A  RANDOM  Shock.  DEFINED  AS  T~E  pi  WEIi5rtT5,  MODEuS  Ea^RESSs. 
“ER^S  are  useful  for  forecasting  AnD/QR  comparison  TO  0~HE-<  Ty  p 


THE  GENERAL  FORM  OF  ThE  MODEL  EXPRtSSsD.  5V  "h£  PI  wE-G-  3 
Z(t>  =  Constant  +  Pi(l)  *  Z<t-D  +  Pi(£)  *  Z  <t-£)  ...  +  P:  in) 


i  LAG 

The.  PI  WEIGhT 

!  Constant 

-. 19453E+01 

Pi  ( 

l) 

-. 39625E+00 

Pi  ( 

2) 

-. 15702E+00 

Pi  ( 

3) 

.  15066E+00 

Pi  ( 

A) 

.  1AA07E+00 

Pi  ( 

5) 

. 90516E-01 

Pi  ( 

7) 

29171E-01 

Pi  ( 

6) 

-. 30830E— 01 

Pi  ( 

9) 

-. 13023E— 01 

Pi  ( 

10) 

£82 AGE +00 

Pi  < 

11)  . 

— . 1 0535E+00 

Pi  ( 

12) 

-. 3897GE-01 

Pi  ( 

13) 

. AA686E— 0 1 

Pi  < 

1A) 

. A01 36E-01 

Pi  ( 

15) 

. 2A201E— 01 

Pi  ( 

20) 

. 15897E+00 

Pi  ( 

21 ) 

. SA803E— 01 

Pi  < 

22) 

. 26A0iE— 01 

Pi  < 

23) 

— . 23382E— 0 1 

Pi  ( 

2A) 

— . 23063E— 01 

Pi  < 

25) 

-. 1A759E-01 

Pi  ( 

30) 

. AAS07E-01 

Pi  < 

31 ) 

. 1670  IE— 01 

Pi  < 

AO) 

-. 25 1 942-0 1 

Pi  < 

41  ) 

-. 10272- -01 

n  iii 


'CDEuInG  ^sS'JL 1 5  -'jR  “I^e  SEkIeS  AXuhI 


DATA  :  Z  =  AXCA112.  1000. 


000  CSSiRVh  " 


DIEi-E.Re.NC ING  ON  Z  :  NONe 


SACKCASTING  :  OFF 
UNIVARIATE  MODEL  PARAMETERS 

**•*•*■#■#*•*■*  ■**•***■***■*■*  ■*•******•**■*■*'***•*#•*  •**■*•***•#***  ***■**■**•***♦  +  ■*■*•■*  ■**■••* 


FACTOR  LAG  COEFFICIENT- 


T  RATIO 


**■**■**■**■*•*•**•*'*■***•**< 

1  MEAN 

2  AUTOREGRESSIVE  1 

3  AUTOREGRESSIVE  1 

4  AUTOREGRESSIVE  1 

5  AuTOREGRESS IVE  2 

6  AUTOREGRESS I VE  2 

7  MOVING  AVERAGE  1 


21326E+01 
£01 1 4E-+-00 
.  1  i07SE-*-00 
.  151£2E-*-00 
.  £S5£7E-*-00 
. 19493E+00 
84166E-01 


-6.  23 

3.  41 

4.  72 
9.  63 
6.  £4 

—2.  49 


r  ■*■*■**■*■*•***■*■  **♦■***•*■«■*■*•***•*•■*•  ■*•-*♦■*♦ 


THE  RESIDUAL  STATISTICS 


Su-M  IJI-  SOuARES 
*EAN  SQUARE 
R  SQUARED 


. 1 1 556E*0 1 
.  12040E-02 
. 22893E+00 


DEGREeS  OF  PREEDu^ 
NUMBER  OF  RES  I  Du:- 


THE  P I  wE I GH i S 


uNIVARIATE  BOX— JENKINS  MQDEuS  CAN  Be  EXAReSSED  AS  A  weIG  — —  eD  S  _  'Y*  _■ 
OluS  A  RANDOM  SHOCK.  DEFINED  AS  T’-E  3 1  wET3H“3,  M0DE..5  E,'  :i::£5S-D  In 
“ERMS  ARE  USEFUu  -qp  FORECASTING  AnD/QR  COMPARISON  TO  C_mE-(  "v  ^Ej  O' 


The  GENERAL.  FGhM  Or-  The  MuDeu  eXAReSSeD  By  t— e  2I  wtl  3-  ‘5  IS 
Z  1 1  )  =  Const  ant  ■+•  3i<1)  *  Z  <  t  —  1  >  +  P  l  (  2  )  *  Z  >.  z  -2 )  ...  -*•  P  i  \  n  .*  * 


u  AG 

T-.j.  o  ‘  wr  I  o“ 

Constant 

- .  ESSESE-^.’’  1 

Pi  <  1 ) 

.  201  l*E-*-00 

Pi  <  2) 

— .  1  1 076S-*-00 

Pi  (  3) 

l5l22E-*-00 

Pi(  10) 

29527E-00 

Pi <  11) 

593S2E— 0 1 

Pi (  12) 

. 32705E-0 1 

Pi  (  13) 

. 4^652E-0l 

P l <  30 ) 

19493E+00 

Pi <  3l ) 

3S206E-01 

Pi<  32) 

. 21590E-01 

Pi (  33) 

. 2S477E— 01 

P i <  40) 

. 84 1 682—0 I 

Pi  <  41) 

16930E— 01 

Pi<  43) 

. 12728E— 01 

Pi <  50) 

. 24858E-01 

Pi (  70) 

. 1 6407E— 0 1 

"GDE-ING  RESULTS  FOR  TIME  SERIES  AXCC112. 


DATA  :  Z  =  AXCC112.  iOOOO.  1000  OrSE  RVA-  I  G'mE 

DIFFERENCING  ON  Z  :  NONE 
SGCKCASTIMG  :  OFF 
UNIVARIATE  MQDEu  PARAMETERS 

*****+***++♦****■**«*******«**♦+♦**•»****«**+♦+♦*■»«+♦♦■*♦*♦  ****♦*■*•*♦ 

FACTOR  LAG  COE" j 1 T  RATIO 


1 

mean 

-.111 13E+0 1 

2 

AUTOREGRESSIVE 

1 

1 

. S5572E-01 

2.  10 

3 

AUTOREGRESSIVE 

2 

. 43537E+00 

.5.  80 

4 

AUTOREGRESSIVE 

4 

3 

.  267£4E-,-00 

8.  £3 

er 

-j 

AuTGREGRESS IVE 

Zj 

10 

.  17603E-*-00 

5 .  bb 

6 

M0VING  AVERAGE 

1 

20 

27£22E+00 

-a.  5: 

*  •*■*•*■#■* ■*■**■#•  »»♦♦+♦  ♦■*•*■•*  ♦-**♦-***-**  *•*♦♦•*  ♦♦-*•*■-»■  * 

THE  RES  I  DUAi_  S~A~IS“ICS 


SLM  OF  SQUARES  :  . 53127E-02 
wEAn  SQUARE  :  .  54156E— 0>5 
R  SQUARED  :  . 51480E+00 


D  r.  G  9  i  n.  5  'w  *“  -  n  Z  C  : 

Ng  '_’V'  ^  P  _  ^  t 


jNlVARIPiE  BOX— JE.NK  I NS  ^GDc.i_S  CHN  at  aXP^aSaaD  hs  P 
3i_'_5  P  FANDOM  St-C'CK.  DEFINED  AS  T-E  a  I  wEIG— _5.  !Y1|_DE_3  E1, 0"': 
^EF+S  ARE  JSEFL'u.  -OR  -GFECPSTInG  A"yD/CF  COMPARISON  tG  G  t-iE-< 


T-E  GENERAL  FORM  OF  T-E  MODEl.  EXPRESSED  BY  T 
Z  (t )  =  Constant  +  Pi(l)  *  2<t  —  i)  +  Pi(2.»  *  Z(t  — 


PG 

THE  P I  Wfc  I  6— 

Const  ant 

1 5SO96+0 1 

Pi  < 

i ) 

65572E-01 

Pi  < 

2) 

435372+00 

Pi  < 

3) 

26764E+00 

Pi  ( 

10) 

1 7Q03E+00 

Pi  ( 

1  1 ) 

1  i67*E~ 01 

Pi  < 

12) 

775092—0 1 

Pi  ( 

13) 

476462— Ol 

Pi  ( 

20) 

272222+0 0 

Pi  < 

21) 

17350E-01 

Pi  ( 

22) 

1 18526+00 

Pi  ( 

23) 

72858E-01 

Pi  ( 

30) 

43464E-0 1 

Pi  ( 

32) 

21 100E-01 

Pi  < 

33) 

1 297 1 s— 01 

Pi  ( 

40) 

74106E— 01 

Pi  < 

42) 

32263E— 0 1 

Pi  ( 

43) 

1983+E-01 

Pi  < 

50) 

13193E— 01 

Pi  ( 

60) 

20 1 73E-01 

"GDE-ING  RESULTS  FOR  TIME  S6R 


****** ***************** ***********  ***  +  ***-i 

QPTA  :  Z  =  AXCP122.  lOOO. 


1000  G35EFVX 


DIF  FERENC InG  On  Z  :  NGNt 
BACKCPS~ING  :  QF- 

***********************************************  ♦***-! 

ON  IVOR  I  PTE  MODEl.  PPRAME'ERS 

******************************************»*****  +  **-i 


FACTOR 

LAG 

COEFFICIENT 

:  RP  i  I  G 

1 

+  EPN 

-. 20727E+0 1 

d 

POTCREGRESSIVE 

1 

1  0 

. 4 1 583E+00 

13.3" 

3 

A'_ir0R2S-E5S  I  YE 

30 

. 260562+00 

9.  52 

4 

HOY  I NG  OVEFPGE 

1 

1 

. 53277E+00 

1  6.  3'- 

~Zj 

+GYING  PYEFPGE 

-j 

127  1 OE +00 

—  u .  2 

6 

HCV I nG  average 

c. 

20 

1 Q495E+00 

-  5 .  ^  r 

*  »  • 


III  If I 


A0-A16S  567  OPERATOR  ALERTNESS/HORKLOAD  ASSESSMENT  USING  STOCHASTIC  2/2 
MODEL-BASED  ANALV  <U>  PERCEPTRONICS  INCUUOOOLAND  HILLS 
CA  A  MADNI  ET  AL  NOV  85  PFTR-1126-85-11 
UNCLASSIFIED  AF0SR-TR-86-8325  F49628-83-C-8881  F/G  5/9  NL 


-I  111 


THE  RESIDUAL.  S~  AT I ST I  CS 


SUM  0-  SQUARES 
■"EAn  SQUARE 
R  SQUARED 


. 51580E+01 
. 53506E-02 
.  52418E+00 


D^GRE-S  u*“  f  : 

number  or  res i d^m^s  s 


THE  PI  WEIGHTS 


jN I  VAR  I ATE  BOX— JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  HE 
US  A  RANDOM  SHOCK.  -  DEFINED  AS  T«E  PI  wEIGh  S,  MODEL 
RMS  ARE  USEFUL  FOR  FORECASTING  A.nD/OR  COMPARISON  TO  OTHER  ~YP£s 


T-e  GEnERAl.  -uRm  OF  T>it  MuDtu  fcXPRESSEO  BY  ~-E  oi  wEIG~~S  I 
Z(t)  =  Constant  P  i  <  1  >  *  Z(t  — 1)  +  ai  (2>  *  Z  (t— 2/  ...  -*•  3i(n; 


lAG 

TnE  PI  Wt,IG^T 

Constant 

-.  49928E-*-01 

Pi  ( 

i> 

-. 53277E+00 

Pi  < 

2) 

— .  28384E-*-00 

Pi  < 

3) 

24120E-01 

Pi  < 

4) 

. 54364E-0: 

Pi  < 

5) 

. 65306E-01 

Pi  ( 

6) 

. 37858E— 01 

Pi  < 

7) 

. 13196E-01 

Pi  < 

10) 

42048E-*-00 

Pi  < 

1  1) 

— . 22386E+00 

Pi  ( 

12) 

-. 1 1857E+00 

Pi  < 

14) 

. 2327 iE— 01 

Pi  < 

15) 

. 27468E-0 1 

Pi  < 

16) 

.  15870E— 01 

Pi  < 

20) 

. 1 830 1E+00 

Pi  ( 

2  i ) 

. 9757*E-0l 

Pi  < 

22) 

. 52279E-01 

Pi  ( 

25) 

-.  1 1949E-01 

Pi  ( 

30) 

-.  2028 1  E-t-00 

Pi  ( 

3 i ) 

-.  10Q0bE-K*0 

Pi  ( 

32) 

57706E-01 

Pi  ( 

34) 

.  liOSOE-Ol 

Pi  ( 

35) 

. 13243E— 01 

Pi  < 

40) 

-.351 36E-01 

Pi  ( 

41  ) 

-.  1868SE-01 

Pi  ( 

50) 

. 37517E-01 

Pi  ( 

51 ) 

. 199Q8E-01 

Pi  < 

52) 

. 1 0670E— 0 1 

IH  1/1 


MODELING  RESULTS  FOR  TIME  SERIES  AXCC122. 


**•**■»**•*■»•*•******■*■»■**■*■****■»***■»■***■•'■•■*•»■»•■***■*••■***-•■**•*♦***■••*+-•■*♦♦*♦■•♦ 

DATA  :  Z  =  AXCC122.  10000.  1000  GBSSRVA”:0N5 

DIFFERENCING  ON  Z  :  NONE 


EfACKCAS'lNG  :  OFF 

***********************************  **+**-**+***+-*-**-«-* 

UNIVARIATE  MODEu  PARAMETERS 

**************************************************■*■■<»•-»-♦♦•»■*-*♦♦-••-»■♦-•■*• 

FACTOR  LAG  COEFFICIENT  T  RATIO 


*******♦**♦********♦*****♦************************♦•*•*♦-»•*  +  -••  +  -•■*-•■-»■-•■♦ 


1  MEAN 

2  AUTOREGRESSIVE  1 

3  AUTOREGRESSIVE  1 

4  AUTOREGRESSIVE  l 

5  MOVING  AVERAGE  1 

6  MOVING  AVERAGE  1 

7  MOVING  AVERAGE  2 


— .  1095lE-*-0l 

1  . 19957E+00 

2  .  10596E+00 

3  . 4a950c+00 

1  .  17681E+00 

2  1C’57E-01 

20  22387E+00 


4.  07 

2.  39 
17.  09 

3.  04 
-.  20 

-7.  00 


********* #♦**■»*******♦■♦•»♦♦**•»**■♦*******■****♦*********♦**+**  *****  * 


T-E  RESIDUAL  STATISTICS 


SUM  OF  SQuARtS 
*EAN  SQUARE 
R  SQuARED 


88485E-02 

89379--05 

4l752E-*-00 


DEGREES  OF  F^EtDo™ 
NUMBER  OR  RrSIDuAuS 


J  III 


THt  PI  wtIGi i S 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGnTED 
"  lS  A  RANDOM  SiiuCK.  DEFINED  AS  T>-*E  =1  wEISi-i_S,  MODE-5  EXPRE 
RTS  ARE  USE-UL  FOR  FORECASTING  AND/OR  COMPARISON  TO  G~HE  ■< 


ThE  GENERAL  FORT  OF  Th£  MODEl  EXPRESSED  BY  T|-tE  °I  hEIGh -S  IS  : 

Z(t)  =  Constant  Pi<l)  *  Z(t-l)  +  Pit£)  ■*  Z  (t-£)  ...  +  Pi  (n)  *  Z(t-n; 


:  lag 

THE  PI  WEIGhT 

i  Constant 

-.  1 9£6 1E+01 

:  Pi <  i) 

— .  37638E+00 

!  Pi (  £> 

-. 16177E+00 

!  Pi t  3) 

5i405E-K>0 

!  Pi (  4) 

Q9148E-01 

i  Pi t  5) 

-. 10£3£E— 01 

!  Pit  £0) 

. ££387E+00 

1  Pi (  £1) 

. Q4£59E— 01 

1  Pit  ££) 

. 36£15E-01 

!  Pit  £3) 

. 1 i506E+00 

!  Pit  £4) 

. 19957E-01 

1  Pit  40) 

-. 501 16E-01 

:  Pit  41) 

-. 18863E-01 

i  Pit  43) 

-. £576£E-01 

!  Pit  60) 

. 1 1E19E-01 

mODEl IMS  RESULTS  FOR  TIME  SERIES  AXCA£l £ 


DATA  :  Z  =  AXCAS1S 


1000  GBSERVA 


DIFFERENCING  ON  Z 


BACKCASTING  :  OFF 


UNIVARIAiE  MODEl  PARAMETERS 

##******♦*■* -#•*•*  ■*•*•*■*-**■*  ■*■**-***■**#*•***  +  +  #  +  *-*  +  +  *  +  **  +  +  -«--*  +  +  +**  +  *  +  ♦♦♦♦HU 

FACTOR  LAG  COEFFICIENT  T  RATIO 


*♦■*•#  •*■!*■  ■*•******■*  +  •#■*♦*■*-»'■****-*•**•*♦  +  ■*■*•*■**  +  *■*-*•*•*■*•*  +  +  +  ■*  +  +  •**■**•«■  *  +  -*  +  +  -l 


1  MEAN 

£  AUTOREGRESSIVE  1 

3  AUTOREGRESSIVE  1 

4  AUTOREGRESSIVE  l 

5  AUTOREGRESSIVE  £ 

6  MOVING  AVERAGE  1 

7  MOVING  AVERAGE  1 
Q  MOVING  AVERAGE  £ 


19370s. +01 
47S£8s+00 
35548E+00 
1 3653E+00 
1987£E-*-00 
53595E-t-00 
i  £0£83E+O0 
39£00t+00 


1.  5£ 
1.  40 
£.  1  £ 
6.  18 
1 . 70 
.  76 
-13.  09 


+*+♦♦  ++***++*♦♦+*+♦♦♦+*♦++♦♦*♦*+♦++♦♦♦♦♦♦♦♦*+»*■♦■»•♦♦*+♦**  +  *♦*+♦-»•** 


i  n  m  i- 


m  r 


THE  RES I DUAL  statistics 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


.  64C37E+00 
.  6541 i E-03 
.  63847E+00 


Di.GREr.5  O'-  “RErDQM  :  37 

NUMBER  UF  RESIDUALS  :  36 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  r- 
uS  A  RANDOM  SHOCK.  DEFINED  AS  T«E  °I  WEIGHTS,  MGDEi_5  EXPRESSED  IN  “ 
RMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHE^  TYPE'S  0“ 


-  3 


MvD- 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESStD  BY  T’-iE  PI  wEIG*-~S  IS  : 

Z  ( t )  »  Constant  >  Pi  (1)  *  Z(t-l)  Pi  (£)  *  Z  (t-£)  ...  +  Pi  <n)  *  I(r-n 


LAG  1  THE  PI  WfIGnT 


!  Constant 

105£0t+0£ 

Pi  < 

i) 

-.  101£6E+0l 

Pi  < 

£) 

-.  1 1010E+01 

Pi  < 

3) 

-.  93£0£E-t-00 

Pi  ( 

4) 

-. 72284E+00 

Pi  < 

5) 

57S45E+00 

Pi  < 

6) 

-. 45556E+00 

Pi  < 

7) 

-. 36106E+00 

Pi  < 

Q) 

-. £859£E+00 

Pi  < 

9) 

-.  ££646E-K>0 

Pi  < 

10) 

-.  37809E-t-00 

Pi  ( 

li) 

-. 15385E+00 

Pi  ( 

1£> 

-.  88507E— 01 

Pi  ( 

13) 

-. 5:51 :E-01 

Pi  ( 

14) 

45559E-01 

Pi  ( 

15) 

-.  3**866E-01 

Pi  < 

16) 

-. £79£7E— 01 

Pi  < 

17) 

-. ££039E— 01 

Pi  ( 

18) 

1 7477E— 0 1 

Pi  ( 

19) 

-.  13837E-01 

Pi  < 

£0) 

. 36104E+00 

Pi  < 

Si) 

. 388£6E-K>0 

Pi  ( 

£2) 

.  4£473E-t-00 

Pi  < 

£3) 

.  3599iE-*-00 

Pi  < 

£4) 

. 27904E+00 

Pi  ( 

£5) 

.  £££55E+00 

:  £9 
Pi<  30) 
Pi <  31) 
32 
Pi (  33) 
Pi  (  34) 
Pi  (  35) 
Pi (  36) 
Pi  (  40) 
Pi  (  41 ) 
Pi<  46) 
Pi <  43) 
Pi  <  44) 
Pi  (  45) 
Pi<  46) 
P i <  47) 
Pi (  46) 
Pi (  49) 
Pi (  50) 
Pi<  51) 
Pi C  52) 
Pi (  60) 
Pi  (  6l) 
Pi <  62) 
Pi <  63) 
Pi (  64) 
Pi (  65) 
Pi <  66) 
Pi <  67) 
Pi <  68) 
P i  (  69) 
Pi (  70) 
P i <  80) 
Pi (  81 ) 
°i(  82) 
Pi <  83) 
=i<  84) 
Pi (  85) 
D  i  <  86) 
Pi < 102) 


.  175886+00 
.  13940E+00 
•  1 10396+00 
. S7437E-01 
. 14715E+00 
. 59467E— 01 
. 34027E— 01 
. 19663E-01 
. 174406-01 
. 13336E— 01 
. 1 0685E— 01 
1 4947E+00 
15229E+00 
1&656E+00 
141146+00 
10943E+00 
87274E— 01 
68970E— 01 
54666E— 01 
43286E— 01 
34288E— 01 
57693E-01 
23319E— 01 
13345E— 01 
. 58594E— 01 
. 59697E-01 
. 65292E-01 
. 55326E— 01 
. 428956-01 
. 34218E— 01 
. 27036E— 01 
. 21429E-01 
. 16969E-01 
. 13441E-01 
. £26166-01 
. 22969E-01 
. £34016-01 
. 25595E-01 
. £16886-01 
.  168156-01 
. 1341 16-01 
.  105986-01 
. 100336-01 


MODELING 


3  " 


♦  *■#•*■*■#•♦•#•♦**♦♦»■*•******■*•****■***■***■*■»•*■*■**■*■*****-*•*•*■**■*•■»•*•*■**♦*■*•■*•****■*••*• 


data  :  Z  =  AXCC212.  10000. 


1000  OBSERVATIONS 


D I F F ERfc.N0 I NG  ON  Z  :  NONE 
BACKCASTING  :  OFF 
jMWR:ptE  MODEu  PARAMETERS 

♦  ♦•s-***********-***********************************************-*-*** 


FACTOR  LAG  COEFFICIENT 


T  -(AT  10 


*♦♦  **********•*■******■»■«■****■»■»**********•*■•■*♦*•»■**  *****•»******■»■■»■***♦  + 


WEAN 

AUTOREGRESSIVE  1 
AUTOREGRESSIVE  1 
•MOVING  AVERAGE  1 
"OV I n*G  AVERAGE  2 


12356E+0 1 
.  42857E-*-00 
. 24043E+00 
. 343922+00 
— .  26843E+00 


4.  93 
G«  l  1 
3.  88 
-8.  £5 


********************************  **************************  *♦***■»•■*•* 


T-iE  RESIDUAL  STATISTICS 

=S=SS=Z^==SS==3Z=S===== 


S'.v'  0-  SQUARES 
MEAN  SQUARE 
R  SQUARED 


. 4? l 33E-0 1 
. 47466E-04 
. 1 5904E+00 


DEGREES  OF  FREEDOM  s 
NUMBtR  OR  Ri-SIDUAlS  : 


THE  PI  WEIGHTS 


M  VAR  I ATE  BOX-JENKINS  MODElS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SuM  OR  T-E 
luS  A  RANDOM  SHOCK.  DEFINED  AS  ThE  PI  WEIGHTS,  MuDElS  EXPRESSED  IN  “i 
ERH5  AR£  lSEFOL  fqR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  0-  M 


■E  GENERAL  c0RM  OF  ThE  MODEl  EXPRESSED  BY  T^e  PI  wEIGr'S  IS  : 


Z  <t  > 


Constant  ♦  Pi(l)  *  Z(t-l)  +  Pi  (2)  *  Z  (t-£) 


+  Pi  (n)  *  Z  \  ~ 


THE  PI  WEIGHT 


Constant 

Pi 

< 

i> 

Pi 

( 

2) 

Pi 

( 

3) 

Pi 

< 

4) 

Pi 

< 

5) 

Pi 

< 

50) 

Px 

< 

5 1 ) 

Pi 

< 

52) 

Pi 

< 

53) 

Pi 

( 

54) 

Pi 

< : 

LOO) 

Pi 

(101  ) 

Pi 

c 

LOS) 

Pi 

(103) 

Pi 

c 

L50) 

Pi 

( 151 ) 

24784c. -*-01 
77E49E+00 
5061 lE+00 
17406E+00 
59864E— 01 
20588E— 01 
26843E+00 
20736E+00 
1 35Q6E-*-00 
46724E— 01 
1 6069E-01 
72057E— 01 
55664E- 01 
36469E— 01 
12542E-01 
19343E-01 
14942E-01 


_ING  RES JL“S  FOR  TIME  SERIES  AXCA22 


*********-*****************************************************"*** 

DA'rA  .  z  =  AXCA222.  1000.  1000  OBSERVATI OnS 

DIFFERENCING  ON  Z  :  NONE 
BACrtCASTING  !  OFF 

******************************************-*********************** 

'_'N  I  variate  MODEl  parameters 

**■*•*•*•***  *•****•*■***■*■*■*  ****************************  ■*■*■*■*■#■•«•■»(•*  ***■»***** 
FACTOR  LAG  COEFFICIENT  T  RATIO 

**■»■**  •»*•»•»**•»•*************  ********  ***************  ***********  **•#•**■* 


MEAN 

AUTOREGRESSIVE 
AUTOREGRESSIVE 
AUTOREGRESSIVE 
AUTOREGRESSIVE 
MOVING  AVERAGE 


18880E+01 
. 13944E+00 
. 97 1 OxE— 01 
. 24565E+00 
.  £324»E-*-00 
1618  1  E-*-00 


4.  41 
3.  06 
7.  81 
7.  37 
-4.  41 


*****•»■**♦*♦*♦**  *♦**♦*■*■****♦****♦■»■*****•***♦*♦***♦  **♦*♦+**♦*♦*****♦ 


THE  RcSIDUAl  STATISTICS 


SUM  of  SQUARES 
m  E  A  \  S  0  l  A  R  E 
R  SQUARED 


. 772042+00 
. 795 1 0E-03 
. 25175E+00 


DEGREES  O1-  FREEDOM 
number  uf  residuals 


T«£  RESIDUAL  STATISTICS 


5l'y1  OF  SQuhRslS  :  . 27636E-01  DEGREES  OF  FREEDuM  : 

MEAN  SQUARE  :  . 28142E-04  NUMBER  OR  RtSIQUAcS  : 

R  SQUARED  :  .  29660E-I-00 


THE  PI  WEIGHTS 


NiIs/PRIA;-  BGX  — JENKINS  MGDiLS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OR  T  —  E 
i__'5  A  RANDuhn  SruCn.  DEFINED  AS  Ti-iE  PI  WEIGHTS,  IMuDEuS  EXPRESSED  IN  T-i 
t  RMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  0-  * 


’  —E  GENERAc  rORM  OF  i  he  MuDEu  EXPRESSED  BY  T*ifc  PI  WEIGHS  IS  : 

Z  <t  >  =  Constant  +  Pi(l>  *  Z(t-l)  Pi  (2)  *  Z  <t-£)  ...  +  Pi  (n;  *  Z'.'-m 


LAG 

THE  PI  WEIGHT 

Constant 

-.  1 35452-t-Ol 

Pi  (  2) 

—. 32224E+00 

Pi  (  3) 

— .  20626E-*-00 

Oi(  10) 

.  14213E+00 

Pi (  12) 

-.  45799E-01 

Pi<  13) 

-. 23316E— 01 

Pi (  50) 

. 31596E+00 

Pi <  52) 

.  lOlSiE+OO 

Pi (  53) 

. 65170E— 01 

Pi <  60) 

-. 44907E-01 

Pi <  62) 

. 1447 IE— 01 

Pi (100) 

-. 9383iE— 01 

Pi < 102) 

-. 32169E— 01 

Pi (103) 

-. 20591E-01 

Pi ( 1 10) 

. 141Q9E— 01 

P l ( 150) 

. 3.543E-01 

P l < 1 52 ) 

.  101S4E— 01 

•*j  r-  •;  iii  lit 

'-0  oj  n  Hi  d 

IT*  in  111  u 


«  -  w  —  V  -  v  *  V  »  w 


V''*  L 


:i » 


Appendix  D  (Cont'd) 


AUTOBJ  SAMPLE  PRINTOUT  (SUBJECT:  JJJ) 


MODELING  RESULTS  FOR  TIME  SERIES  JJJAATRE 


DATA  s  Z  *  JJJAATRE  1000. 
DIFFERENCING  ON  Z  i  1 )  1  OF  ORDER  1 
BACKCASTING  :  OFF 


UNIVARIATE  MODEL  PARAMETERS 


1000  OBSERVATIONS 


x.  I  I  X.i 


.  X  X2XXXXXXI 


;■  T'  x  X  Z  X'1 


FACTOR  LAG  COEFFICIENT 


1  AUTOREGRESSIVE  1 

2  MOVING  AVERAGE  1 

3  MOVING  AVERAGE  1 

4  MOVING  AVERAGE  2 


. 19333E+00 
. 97338E+00 
661 48E—01 
17422E+00 


T  RATIO 


6.  07 
41.  12 
-2.  67 
-3.  36 


THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


46S52E+0 1 
49489E-02 
68648E+00 


DEGREES  OF  FREEDOM 
NUMBER  OF  RESIDUALS 


THE  PI  WEIGHTS 


UN I  VARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  T*E  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  Or  MODEuS 

THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

Z<t)  =  Constant  +  Pi(l)  *  Z<t-1)  +  Pi(£)  *  Z  (t-2)  ...  +  Px(n)  *  Z(t-n) 


‘""lag 

!  THE  PI  WEIGHT 

Constant 

. OOOOOE+OO 

Pi  ( 

1) 

19756E+01 

Pi  ( 

2) 

19273E+01 

Pi  < 

2) 

-..18141E+01 

Pi  ( 

4) 

1 639 1E+01 

Pi  < 

3) 

14716E+01 

Pi  < 

6) 

13lS7E-t-01 

Pi  < 

7) 

11751E+01 

Pi< 

8) 

10491E+01 

Pi  < 

9) 

93643E+00 

Pi  < 

10) 

83383E+00 

Pi  ( 

11) 

74602E+00 

Pi  ( 

12) 

663S6E+00 

Pi  ( 

13) 

59431E+00 

Pi  ( 

14) 

53045E+00 

Pi  ( 

IS) 

47343E+00 

Pi  < 

18) 

42258E+00 

Pi  < 

17) 

37717E+00 

Pi  C 

18) 

- . 33664E+00 

Pi  < 

19) 

30047E+00 

Pi  ( 

20) 

-.46151E+00 

Pi  < 

21) 

23464E+00 

Pi  C 

22) 

20904E+00 

Pi  ( 

23) 

17340E+00 

Pi  < 

24) 

1 5363E+00 

Pi  ( 

23) 

13607E+00 

Pi  < 

26) 

12127E+00 

Pi  ( 

27) 

108 15E+00 

Pi  < 

28) 

— .  96308E— 01 

Pi  ( 

29) 

86129E— 01 

Pi  ( 

30) 

— .  76872E— 01 

Pi  ( 

31) 

686 1  IE— 01 

Pi  ( 

32) 

6 1238E— 01 

Pi  ( 

33) 

54657E-01 

Pi  ( 

34) 

— . 48784E—01 

Pi  < 

33) 

43542E-01 

Pi  < 

36) 

38863E— 01 

Pi  < 

37) 

34687E-01 

Pi  ( 

38) 

30960E— 01 

Pi  < 

39) 

27633E-01 

Pi  ( 

40) 

24664E— 01 

Pi  < 

41) 

22O14E-01 

Pi  ( 

48) 

1  9648E  — 01 

Pi  ( 

43) 

17337E-01 

Pi  ( 

44) 

1 5632E— 01 

Pi  < 

43) 

13970E-01 

Pi  ( 

46) 

12469E-01 

Pi  < 

47) 

Ill  29E-01 

Pi  < 

30) 

. 16631E+00 

Pi  (  51) 
Pi (  52) 
Pi <  53) 
Pit  54) 
Pi <  55) 
Pit  56) 
Pit  57) 
Pit  52) 
Pi t  59) 
Pit  SO) 
Pit  61) 
Pit  62) 
Pi t  62) 
Pit  64) 
Pi t  65) 
Pit  66) 
Pi t  67) 
Pit  66) 
Pi f  69) 
Pit  70) 
Pi f  71) 
Pit  72) 
Pit  73) 
Pit  74) 
Pit  75) 
Pit  76) 
Pi  t  77) 
Pit  76) 
Pi t  79) 
Pit  80) 
Pi  t  81) 
Pit  82) 
Pi tlOO) 
Pi t 101 > 
Pi t 102) 
Pi tl03) 
Pi (104) 
Pi t 105) 
Pi (106) 
Pi (107) 
Pi (108) 
Pi ( 109) 
Pi (110) 
Pi  ( 111 ) 
Pi  (112) 
Pi ( 113) 
Pi (114) 
Pi (115) 
Pi (1 16) 
Pi (117) 
Pi  (118) 
Pi (120) 
Pi (151) 
Pi ( 152) 


. 33712E+00 
. 32948E+00 
. 31043E+00 
. 28055E+ 00 
. 25190E+00 
. 22522E+ 00 
. 201 16E+00 
. 1 7958E+00 
. 16030E+00 
. 14308E+00 
. 12771E+00 
. 11396E+00 
. 10174E+00 
. 90804E— 01 
. 81046E-01 
.72338E-01 
. 64565E-01 
.57627E-01 
.  51434E-01 
. 79589E— 01 
.40152E-01 
. 35769E—01 
.  29631E-01 
. 2625 IE— 01 
. 23244E-01 
.  20717E-01 
-*_liL474E-*0i 
.  16465E-01 
. 14713E-01 
. 13131E-01 
. 11720E-01 
. 10461E-01 
2900 IE-01 
56757E-01 
57423E-01 
.54102E-01 
..  48894E— 01 
. 43902E-01 
■.  39251E-01 
. 35058E— 01 
31298E-01 
. 27937E-01 
■.  24936E-01 
. 22257E— 01 
■.  19865E-01 
.  1 773  IE— 01 
.  15825E-01 
. 14125E— 01 
. 12607E-01 
. 1 1252E-01 
. 10043E-01 
.  1 3869E— 01 
.  10237E— 01 
. 10004E-01 


353S3  =  : 


MODELING  RESULTS  FOR  TIME  SERIES  JJJCATRE 


*4****#*HHHHHHHHHHHH»*HHHHHHHHHH»**-HHHt***#-»-**-***-»**-»**.**** 

DATA  :  Z  =  JJJCATRE  10000.  1000  OBSERVATI ONS 

DIFFERENCING  ON  2  s  1 )  1  OF  ORDER  10 
8ACKCASTING  s  OFF 
UNIVARIATE  MODEL  PARAMETERS 

FACTOR  LAG  COEFFICIENT  T  RATIO 

\ 


1 

AUTOREGRESSIVE 

1 

10 

-.  68971E+00 

-24. 07 

a 

MOVING 

AVERAGE 

1 

1 

. 2S697E+00 

9.  47 

3 

MOVING 

AVERAGE 

1 

3 

-.  17975E+00 

-5.  92 

4 

MOVING 

AVERAGE 

2 

20 

. 43790E+00 

12.  25 

5 

TREND  CONSTANT 

-.  81478E— 03 

-7.  61 

THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES  :  . 38776E-01  DEGREES  OF  FREEDOM  :  375 
MEAN  SQUARE  s  .  39770E-04  NUMBER  OF  RESIDUALS  :  960 
R  SQUARED  s  .  84003E+00 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  3A3~ 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  'MESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODELS 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 


Z(t) 


Constant  -*•  Pi(l)  *  Z(t-l)  +  Pi(2)  *  Z (t-2) 


+  P  l  (  n  )  *  Z  ( t  -  n  ) 


LAG 

THE  PI  WEIGHT 

Constant 

16236E— 02 

Pi  (  1) 

28697E+00 

Pi  (  2) 

82350E— 01 

Pi  <  3) 

. 13612E+00 

Pi  (  4) 

. 96383E-01 

Pi  <  5) 

. 42461E-01 

Pi  (  6) 

15877E— 01 

Pi  (  7) 

218S1E— 01 

Pi  (  a> 

13912E-01 

‘  Pi (  10) 

30668E+00 

Pi  <  11) 

85507E— 01 

Pi (  12) 

~ .  24333E— 01 

Pi (  13) 

. 48143E— 01 

Pi (  14) 

. 29185E-01 

Pi <  15) 

. 12749E-01 

Pi <  20) 

1 1265E+01 

Pi  (  21) 

32251E-MX) 

Pi (  22) 

92495E-01 

Pi (  23) 

. 1 7595E+00 

Pi <  24) 

. 10846E+00 

Pi (  25) 

. 4775 IE— 01 

Pi <  26) 

17923E— 01 

Pi <  27) 

— . 24639E—01 

Pi (  28) 

15654E-01 

Pi (  30) 

13181E+00 

Pi (  31) 

3501  lE—Ol 

Pi <  33) 

. 20875E— 01 

Pi <  34) 

. 12284E-01 

Pi (  40) 

49329E+00 

Pi (  41) 

14124E+00 

Pi (  42) 

40513E-01 

Pi <  43) 

.  77044E-01 

Pi (  44) 

. 47497E— 01 

Pi <  45) 

. 20912E-01 

Pi (  47) 

1 0790E— 01 

Pi (  50) 

577 19E-01 

Pi (  51) 

15331E-01 

Pi (  60) 

21600E>00 

Pi (  SI) 

61845E-01 

Pi (  62) 

17748E-01 

Pi <  63) 

. 33733E-01 

Pi (  64) 

. 20797E-01 

Pi (  70) 

— .  25275E— 01 

Pi <  80) 

94585E— 01 

Pi (  81) 

27082E— 01 

Pi (  83) 

.  14771E-01 

Pi (  90) 

1 1068E— 01 

Pi (100) 

41418E-01 

Pi (  101  ) 

1  1859E— 01 

Pi <120) 

18156E-01 

V 


MODELING  RESULTS  FOR  TIME  SERIES  JJJA111. 


DATA  s  Z  -  JJJA1 11.  1000. 


1000  OBSERVATIONS 


DIFFERENCING  ON  Z  s  NONE 
BACKCASTING  :  OFF 
UNIVARIATE  MODEL  PARAMETERS 

M*«**HH>**#**««*«««'»«*«'»'»44*** ***■»•»■»■»■**■»**  *** 

FACTOR  LAG  COEFFICIENT  T  RATIO 


1  MEAN 

2  AUTOREGRESSIVE  1 

3  AUTOREGRESSIVE  1 

4  AUTOREGRESSIVE  1 

5  AUTOREGRESSIVE  2 

6  MOVING  AVERAGE  1 


. 13158E+01 

1  50671E+00 

2  -.  266S7E+00 
4  . 60033E— 01 

10  .  39674E+00 

30  -. 13514E+00 


-16.  46 
-8.  67 
2.  1 1 
9.  38 
-4.  11 


THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


.  37846E+01 
. 390 29 E -02 
. 30725E+00 


DEGREES  OF  FREEDOM  : 
NUMBER  OF  RESIDUALS  : 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MGDE—S 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 


Z  <t  > 


Constant  +  Pi(l)  *  Z(t-l)  +  Pi(£)  *  Z(t-2) 


+  Pi (n)  *  Z(t-n; 


MODELING  RESULTS  FOR  TIME  SERIES  JJJC111 


DATA  :  2  =»  JJJC1U.  10000.  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  :  NONE 
BACKCASTING  s  OFF 
UNIVARIATE  MODEL  PARAMETERS 

FACTOR  LAS  COEFFICIENT  T  RATIO 


1  MEAN 

2  AUTOREGRESSIVE  1 

3  AUTOREGRESSIVE  1 

4  AUTOREGRESSIVE  1 

5  AUTOREGRESSIVE  2 


29Q26E+0 1 
. 84105E-01 
.  23309E+00 
. 17367E+00 
.  Ill 33E+00 


2.  £6 
8.  30 
5.  49 
3.  49 


THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


. 8975SE-0 1 
. 91401E-04 
.  13713E+00 


DEGREES  OF  FREEDOM 
NUMBER  OF  RESIDUALS 


3 

U 

i": 


I  K 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  ThESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODELS. 


9 


m  > 


r.~' 


V. 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

Z(t>  »  Constant  +  Pi<l)  *  Z  (t-1)  ♦  Pi(2>  *  Z  <t-2)  ...  ♦  Pi  (n>  *  Z  <t-n) 


LAG 

THE  PI  WEIGHT 

Constant 

-. 46740E+01 

Pi  ( 

1) 

.  -. 84105E-01 

Pi  < 

2) 

-. 2S509E+00 

Pi  ( 

3) 

-.  17367E+00 

Pi  < 

10) 

-.  1  U35E+0O 

Pi  < 

12) 

. 28403E-01 

Pi  < 

13) 

. 19339E-01 

MODELING  RESULTS  FOR  TIME  SERIES  JJJA121. 


Rsasssssasssssssnssssssssss 


DATA  :  Z  »  JJJA121.  1000.  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  :  NONE 
BACKCASTING  s  OFF 

♦♦***■*#■**•»•****■*****•**■****#**•»*•***-»**■*♦*•**•*•****•*-***■»♦♦**•»*•*■»**•♦**■» 
UNIVARIATE  MODEL  PARAMETERS 

FACTOR  LAG  COEFFICIENT  T  RATIO 


1  MEAN 

2  AUTOREGRESSIVE  1 

3  AUTOREGRESSIVE  1 

4  AUTOREGRESSIVE  2 
3  MOVING  AVERAGE  1 


. 12344E+01 

1  -. S4072E+00 

2  -.  27013E+00 

10  . 66730E+00 

20  -.  1S308E+00 


-17. 60 
-8.  73 
23.  10 
-4.  37 


THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


. 29880E+0 1 
. 30337E-02 
. S6930E+00 


DEGREES  OF  FREEDOM 
NUMBER  OF  RESIDUALS 


983 

988 


■  mA  i  mA  .  ^  ^  ^  <  >A,.  .i4  .  ^  I  . 


THE  PI  WEIGHTS 


UNIVARIATE  BOX “JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  P 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN^THES 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  Or  MODELS 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

Z(t)  *  Constant  Pi<l)  *  Z(t  —  1)  +  Pi(2)  *  Z<t-2)  ...  *  Pi(n)  *  Z(t-n) 


LAG 

% 

THE  PI  WEIGHT 

Constant 

. 14310E+01 

Pi  <  1) 

.  54O72E+00 

Pi  <  2) 

. 27013E+00 

Pi (  10) 

-. 66750E+00 

Pi  (  11) 

-.  36093E+00 

Pi <  12) 

-.  18031E+00 

Pi (  20) 

. 1650GE+00 

Pi (  21) 

-.  Q9262E— 01 

Pi <  22) 

—. 44592E— 01 

Pi (  30) 

.  11019E+00 

Pi <  31) 

. 59582E— 01 

Pi (  32) 

.  29765E-01 

Pi (  40) 

-. 2725 1E—01 

Pi (  41) 

. 14735E-01 

Pi <  50) 

-.  18190E-01 

MODELING  RESULTS  FOR  TIME  SERIES  JJJC121 


**-*•**•»*#♦■*•»•»*■***■*****■***■*♦*■*♦♦■**'»■**•******•**♦**#■»•**•******•»♦*•*•**■»*■* 
DATA  :  Z  -  JJJC121 .  10000.  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  s  NONE 

BACKCASTING  :  OFF 

UNIVARIATE  MODEL  PARAMETERS 

****************************  **■*■»•*-»■**■#■*••*-**•*•»•■*-*■**■*+*■*■»•*-**■*■*•*-**+■»  •*■*■»• 
FACTOR  LAG  COEFFICIENT  T  RATIO 

******************************************************************* 

1  MEAN  -. 29711E+01 

2  AUTOREGRESSIVE  1  2  .140432+00  4.45 

3  AUTOREGRESSIVE  1  3  .  1S761E+00  5.33 

4  AUTOREGRESSIVE  1  4  . 79626E-01  2.52 

5  AUTOREGRESSIVE  2  10  .  93544E-01  2.93 

6  AUTOREGRESSIVE  2  20  .  15126E+00  4.76 

7  MOVING  AVERAGE  1  30  -.  14436E+00  -4.44 

*************************  ******* *******************-»*■«.*+ *-***♦ 


cl  Id 


THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES  :  . 12944E+00 
MEAN  SQUARE  :  . 1333SE-03 
R  SQUARED  :  .  11919E+00 


DESREES  OF  FREEDOM  :  3£3 

NUMBER  OF  RESIDUALS  :  376 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  ThE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODELS. 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

Z(t)  »  Constant  ♦  Pi(l)  *  Z(t-l)  +  Pi<2)  *  Z(t-2>  ...  +  PiCn)  *  Z(t-n) 


LAG  i  THE  PI  WEIGHT 


!  Constant 

:  Pi  ( 

2) 

!  Pi  ( 

3) 

i  Pi  < 

4) 

1  Pi  < 

10) 

!  Pi  ( 

12> 

:  Pi  < 

13) 

!  Pi  ( 

20) 

!  Pi  < 

22) 

!  Pi  < 

23) 

:  Pi  < 

24) 

:  Pi  < 

30) 

!  Pi  < 

32) 

;  Pi  ( 

33) 

!  Pi  < 

34) 

:  Pi  ( 

40) 

:  Pi  < 

50) 

!  Pi  < 

60) 

-.  3392 1E+01 
14O4QE+00 
-.  l6761E*0O 
79626E-01 
-.  93344E-01 
. 13I41E-01 
. 1 5679E-01 
-.  15126E+00 
. 21249E— 01 
. 25332E-01 
. 12044E-01 
. 14436E+00 
. 202S0E— 01 
.  24195E-01 
. 1 1495E-01 
. 13504E-01 
. 21835E-01 
-.  20Q39E-01 


MODELING  RESULTS  FOR  TIME  SERIES  JJJA211. 


DATA  :  Z  =  JJJA21 1 , 


1000. 


1000  OBSERVATIONS 


DIFFERENCING  ON  Z  :  NONE 


BACKCASTING  :  OFF 

**  ♦♦****+*******************#****#**-»**4********+***********  ■**  +  ■** 

UNIVARIATE  MODEL  PARAMETERS 


FACTOR 


LAG  COEFFICIENT 


T  RATIO 


*****  *************************************************************  * 


1  MEAN 

2  AUTOREGRESSIVE  1 

3  AUTOREGRESSIVE  1 

4  MOVING  AVERAGE  1 

5  MOVING  AVERAGE  1 

6  MOVING  AVERAGE  1 

7  MOVING  AVERAGE  2 


. 131 54E+01 
. 2509QE+00 
.  17444E+00 
. 42Q55E+00 
.  7331  BE -01 
90160E— 01 
69039E-01 


7.  96 
5.  62 
13.  42 
2.  11 
-2.  80 
-2.  12 


»+*HHHHt***»**##***M*»*»**»****»**#*****«***t***«*4**»»********** 

THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


. 50976E+01 
. 52391 E-02 
. 27477E+00 


DEGREES  OF  FREEDOM 
NUMBER  OF  RESIDUALS  : 


THE  PI  WEIGHTS 


UNIVARIATE  BOX -JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  PAS 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MGDE^ 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 


Z  (t  > 


Constant  Pi(i>  *  Z(t-l)  Pi(2)  *  Z(t-2) 


+  Pi(n)  *  Z ( t - n) 


* 

I 

LAG 

THE  PI  WEIGHT 

Constant 

.  2S614E-f01 

Pi  <  1) 

42S55E+00 

Pi  <  2) 

25697E+00 

Pi  <  3) 

31387E— 0 1 

Pi  (  5) 

. 1844SE— 01 

Pi  (  6) 

. 12376E— 01 

Pi  (  10) 

-.2S130E+00 

Pi (  11) 

10737E+00 

Pi <  12) 

6468 1E-Q1 

Pi t  13) 

12971E-01 

Pi (  20) 

17456E+00 

Pi <  21) 

74656E— 01 

Pi  t  22) 

— . 4487SE— 01 

Pi t  40) 

. 69039E— 01 

Pit  41) 

.29567E-01 

Pit  42) 

. 17741E-01 

Pit  50) 

. 17359E-01 

Pit  60) 

. 12054E-01 

- - -  _  -  _  .  _  , - .  — 

MODELING  RESULTS  FOR  TIME  SERIES  JJJC2U. 


DATA  ;  Z  -  JJJC2U.  10000.  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  :  NONE 
BACKCASTING  :  OFF 

***♦*•*•*■*  ♦*■*■*■*•*•**■*•*■*•**-*•#-*■»■*  •*■*■*■*■*  ■»■*■*■*■*-*•»-******■»■**■»•)►*-»*■****■**•»  ■**■*■*** 
UNIVARIATE  MODEL  PARAMETERS 


FACTOR 

LAG 

COEFFICIENT 

T  RATIO 

#■*•#*♦ ■*•#•*■*■*•»•**•*■*•*■*■**■**■*•****■#•*■#*■*■*■*•*■*  *  ■***■■*•*■*-*+ ■*-*-*-*-**-*  **■*■ 

1  MEAN  21062E+01 

2 

AUTOREGRESSIVE 

1 

1 

. 10464E+00 

3.  32 

3 

AUTOREGRESSIVE 

1 

2 

. 20734E+00 

6.  71 

4 

AUTOREGRESSIVE 

1 

3 

. 1 7823E+00 

5.  £5 

3 

AUTOREGRESSIVE 

2 

10 

.  1433SE+00 

4.  46 

6 

AUTOREGRESSIVE 

2 

20 

. 12431E+00 

3.  86 

7 

MOVING  AVERAGE 

1 

30 

10870E+00 

-3.  29 

THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES  :  . 10776E+00  DEGREES  OF  FREEDOM  : 
MEAN  SQUARE  :  . 11109E-03  NUMBER  OF  RESIDUALS  : 
R  SQUARED  :  . 16005E+00 


THE  PI  WEIGHTS 


u'NIVPRIATE  SOX-JENKINS  MODELS  CAM  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  TiE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TG  OTHER  TYPES  OF  MGDE-S. 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

Ztt)  *  Constant  ♦  Pi(i)  *  Z(t-l)  *  Pi<2)  *  Z<t-2>  ...  ♦  Pi(n)  *  Z  (t-n) 


1 

1 

LAG 

THE  PI  WEIGHT 

Constant 

-.  309O3E+01 

Pi  <  1) 

10464E+00 

Pi  t  2> 

-.  20734E-*-00 

l  Pit  3) 

-.  17823E+00 

Pit  10) 

-. 14339E+00 

Pit  11) 

.  1S004E-01 

Pit  12)  • 

. 29739E-01 

Pit  13) 

.  2355BE— 01 

Pit  20) 

-.  12431E+00 

Pit  21) 

. 13008E-01 

Pit  22) 

. 23800E— 01 

Pit  23) 

. 22137E-01 

Pit  30) 

. 10870E+00 

Pi t  31) 

.  11374E-01 

Pit  32) 

. 22360E-01 

!  Pit  33) 

. 19374E-01 

Pit  40) 

. 15387E— 01 

Pit  SO) 

.  13313E-01 

Pit  60) 

1 181SE-01 

MODELING  RESULTS  FOR  TIME  SERIES  JJJA221. 


DATA  :  Z  *  JJJA221.  1000.  lOOO  OBSERVATIONS 

DIFFERENCING  ON  Z  s  NONE 
BACKCASTING  :  OFF 

*#-*■*-*■*♦-*■**■***•»•*♦■**■*■**■*■*■*■*-*■»■*■»****•*■*-*■***■**■»•»*■*•*♦•*■»*■»**•*  *•*■***•**■* 
UNIVARIATE  MODEL  PARAMETERS 

FACTOR  LAG  COEFFICIENT  T  RATIO 


1 

MEAN 

13395E+01 

2 

AUTOREGRESSIVE 

1 

2 

. 1 1936E+00 

63 

3 

AUTOREGRESSIVE 

1 

3 

.  29663E+00 

9.  73 

4 

AUTOREGRESSIVE 

1 

4 

. 2004 1E+00 

6.  58 

3 

AUTOREGRESSIVE 

2 

10 

. 4961 1E+00 

16.  34 

6 

MOVING  AVERAGE 

1 

1 

. 46162E+00 

14.  44 

7 

MOVING  AVERAGE 

2 

30 

-.  23880E+00 

-7.  47 

iAt.Uk 


THE  RESI DUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


. 58201E+0 1 
. 59449E-02 
. 56236E+00 


DEGREES  OF  FREEDOM 
NUMBER  OF  RESIDUALS 


THE  PI  WEIGHTS 


G73 


386 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODE-S. 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

.  2  ( t  >  =  Constant  +  Pi<l)  *  Z<t-1)  +  Pi<2>  *  Z <t-2>  ...  ♦  Pi(n)  *  Z<t-n) 


LAG 

THE  PI  WEIGHT 

Constant 

36289E+01 

Pi  <  1) 

46162E+00 

Pi  (  2) 

-. 33266E+00 

Pi  (  3) 

45019E+00 

Pi  (  4) 

-. 40822E+00 

Pi  (  5) 

-.  1B844E+00 

Pi  (  6) 

86990E— 01 

Pi  <  7) 

-.  40157E-01 

Pi  <  8) 

18537E-01 

Pi<  10) 

50006E+00 

Pi  <  11) 

-. 23084E+00 

Pi <  12) 

-.  47245E-01 

Pi(  13) 

. 12535E+00 

Pi  <  14) 

.  15729E+00 

Pi  <  15) 

. 72608E— 01 

Pi <  16) 

.  33517E-01 

Pi  (  17) 

. 15472E-01 

Pi (  30) 

. 23880E+00 

Pi  <  31) 

. 1 1023E+00 

Pi (  32) 

.  79437E-01 

Pi  <  33) 

. 10750E+00 

Pi  <  34) 

. 37482E-01 

Pi <  35) 

. 45000E-01 

P i (  36) 

. 20773E-01 

P  i  <  40 ) 

. 1 1941E+00 

Pi (  41) 

.  55124E-01 

Pi  <  42) 

. 1 1282E-01 

P  i  (  43) 

-.  29933E-01 

Pi<  44) 

37560E-01 

Pi<  45) 

-.  17339E-01 

Pi  <  60) 

57024E-01 

Pi  (  61 ) 

-.  26323E-01 

Pi<  62) 

18969E-01 

Pi<  63) 

-.  25671E-01 

Pi  (  64) 

-.  23278E-01 

Pi<  65) 

-.  10746E-01 

Pi<  70) 

-. 28516E-01 

Pi <  71 ) 

-.  13163E-01 

Pi<  SO) 

, 13617E-01 

MODELING  RESULTS  FOR  TIME  SERIES  JJJC221. 


DATA  :  2  »  JJJC221.  10000.  1000  OBSERVATIONS 

DIFFERENCING  ON  2  :  NONE 
SACKCASTING  ;  OFF 

#*-*•»•* *****•»■*■»■***■*•#*■»•***•»■»-»*■*■*-»•*■»+■*■»■*■»*•»■»*■*■»■»■*■»■»■*■»•*  ■*•*•»•***■**• *•*•**  *■*■* 
UNIVARIATE  MODEL  PARAMETERS 

#•#•#**•#*•#■**■*■»■*■*•*■*•#■»■*■*■»•*■**•*■**■*-»■»•*•*■#■*■#■*■#*•*■***-*■** •*■*■*■**•*•*•**•**•»•'# -»•»-»•»  -*■#•+ 
FACTOR  LAG  COEFFICIENT  T  RATIO 

mm*****M*Mm*#**m«*«#**#M*+*+M******m+*********4#** 


1 

MEAN 

19154E+01 

2 

AUTOREGRESSIVE 

1 

2 

.  14333E+00 

4.  58 

3 

AUTOREGRESSIVE 

1 

3 

.  13359E+00 

4.  27 

4 

AUTOREGRESSIVE 

2 

10 

. 71010E-01 

2.  23 

3 

MOVING  AVERAGE 

1 

70 

S6921E-01 

-2.  65 

THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES  :  .  12371E+00  DEGREES  OF  FREEDOM  :  S82 
MEAN  SQUARE  :  . 13209E-03  -  NUMBER  OF  RESIDUALS  :  987 
R  SQUARED  :  . 32182E-01 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  PA 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODE, 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

Z(t)  =»  Constant  +  Pi<l)  *  Z(t-l)  +  Pi(2)  *  Z<t-2)  ...  +  Pi  (n)  *  Z(t-n) 


.  . 

LAG 

THE  PI  WEIGHT 

Constant 

— . 23407E+01 

Pi  < 

2) 

14333E+00 

Pi  ( 

3) 

-. 13359E+00 

Pi  ( 

10) 

710 10E— 01 

Pi  ( 

12) 

. 10178E-01 

Pi  ( 

70) 

. 86921E-01 

Pi  < 

72) 

. 12459E-01 

Pi  < 

73) 

.  11612E-01 

_____ 

MODELING  RESULTS  FOR  TIME  SERIES  JJJA212 


•»  ■****•*•♦•»■*•*■*••*•*■*■*•*■»  ■*•*■** *■»■»+■**++■»■»■*•»■*+■»■»■»■»*+  ■*■*■*•*■*■*■**■*■#■*■*■*■*■*  ■**•*■*** 


DATA  :  Z  »  JJJA212.  1000. 

DIFFERENCING  ON  Z  :  NONE 


1000  OBSERVATIONS 


BACKCASTING  :  OFF 

#***-*•#*-*■******♦*•*  HH»-********-**********-IH» -*■*-*-**■*•****+*  **-****-)H»*-*-*  ■*■*•■* 

UNIVARIATE  MODEL  PARAMETERS 

*-*■*-*-#■*■#*•»•*■* ■*•*■*■**■»■*-#■* *■*-#■*■*-*•#■*■*•*-*■*•*■#■*•*■*•**■*■#■*■*■*•* **•*■*■*•■*+*■*■**-*•*•#. ■»*+ 
FACTOR  LAG  COEFFICIENT'  T  RATIO 

*•»■»•»*****•»*-***■)►*  ^ •*-*-»•*■*■#•**■*■#■*■»■*■* ***♦**■*•*•*  ■*•*•*■**•*■**-*■»*■*  ■#■■**•**-*  *■*-*•  ***** 


1  MEAN 

2  AUTOREGRESSIVE  1  1 

3  AUTOREGRESSIVE  1  2 

4  AUTOREGRESSIVE  2  10 

5  MOVING  AVERAGE  1  10 


53939E+0 1 
sassoE+oo 
37125E+00 
.89172E+00 
.  49908E+00 


-19. 76 
-12.  54 
45.  31 
13.  08 


♦♦♦♦♦♦♦♦♦♦♦♦♦♦-a-#*********-*****-#***-***************-**  ■*-#■-*-***•*  ■*-*■*-#•■*■*•* 


THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


. 5SS62E+0 1 
. 59880E-02 
. 66509E+00 


DEGREES  OF  FREEDOM  :  9 32 

NUMBER  OF  RESIDUALS  :  9S8 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SLIM  OF  THE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MCDE-S. 

THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

Z  ( t )  .»  Constant  >  Pi(l)  *  Z(t-l)  +  Pi<2)  *  Z(t-2)  ...  *  Pi  (n)  *  Z(t-n) 


I 


LAG 

1  THE  PI  WEIGHT 

Constant 

1 9254E+02 

Pi  (  1) 

.  38350E-KJ0 

-Pit  2) 

. 37123E+00 

Pit  10) 

13908E+01 

Pit  11) 

229S3E+00 

Pit  12) 

14377E+00 

Pit  20) 

-.69411E+00 

Pit  21) 

11473E-MX) 

Pit  22) 

-.72749E-01 

Pit  30) 

34S41E+00 

Pit  31) 

57260E—01 

Pit  32) 

36307E-01 

Pit  40) 

172S9E+00 

Pit  41) 

2B377E— 01 

Pit  42) 

1S120E— 01 

Pit  30) 

862S3E— 01 

Pit  51 ) 

142S2E— 01 

Pi  t  SO) 

43062E-01 

Pit  70) 

— . 21491E— 01 

Pit  SO) 

10726E-01 

MODELING  RESULTS  FOR  TIME  SERIES  JJJC212. 


~  #♦**■*•**■»*•**•»*#■»■*■*#*■*■»•»*•»•»■»♦•»#•»■*♦■*■»*■»■»*■»**■»•**•»■»•»****•*•»■**■••+*■***■**■*■» 
DATA  s  Z  ■  JJJC212.  100.  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  s  NONE 

BACKCASTING  :  OFF 

UNIVARIATE  MODEL  PARAMETERS 

*■»**■*■**■»■******#*■*#*■*•»*■»**■»•**■*■*  *■**•**-*-*■ -»*-**•* -*-*-*•  *•*-*•*-*•#•■*■»•■*•*•*■*-**  *-*+-*■-*+ 
FACTOR  LAG  COEFFICIENT  T  RATIO 


1 

MEAN 

60968E+01 

2 

AUTOREGRESSIVE 

1 

1 

.  10810E+00 

3.  37 

3 

AUTOREGRESSIVE 

1 

2 

.  25623E+00 

8.  32 

4 

AUTOREGRESSIVE 

1 

3 

.  14199E+00 

4.  49 

3 

AUTOREGRESSIVE 

2 

10 

.  24391E+00 

7.  75 

6 

AUTOREGRESSIVE 

2 

20 

.  13294E+00 

4.  18 

7 

MOVING  AVERAGE 

1 

50 

S4S24E-01 

-2.  53 

THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


. 73506E+03 
. 7S779E+00 
.  1343SE+00 


DEGREES  OF  FREEDOM 
NUMBER  OF  RESIDUALS 


370 

377 


I 


f.- 

tv' 

►  ' 


✓  . 


THE  PI  WEIGHTS 


UN  I  VAR  I ATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  PS  A  WEIGHTED  SUM  Or  TnE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MCDE-S. 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

ZCt)  =  Constant  *  Pi(l)  *  Z<t-1>  +  Pi(2)  *  Z  ft-2)  -  +  Pi  (n)  *  Z(t-n) 


It 


LAG 

THE  PI  WEIGHT 

Constant 

-.  95185E+01 

Pi  (  1) 

-.  10810E+00 

Pi  (  2) 

-. 25623E+00 

Pi  <  3) 

-.  14199E+00 

Pi (  10) 

-.24591E+00 

Pi <  11) 

.  2S5S3E-01 

Pi <  12) 

.6301  IE-01 

Pi  <  13) 

. 34917E-01 

Pi<  20) 

-.  13294E+00 

Pi (  21) 

.  14370E-01 

Pi  (  22) 

. 34063E-01 

Pi  <  23) 

.  18876E— 01 

Pi <  50) 

.84624E-01 

Pi (  52) 

. 21664E-01 

Pi (  S3) 

. 12016E-01 

Pi (  60) 

.20810E-01 

Pi  (  70) 

. 1 1250E-01 

MODELING  RESULTS  FOR  TIME  SERIES  JJJA222. 


*****************•****+**■)«*+*■****+************♦**■*•*■*■*■***•**-*•***■#•■*•»• 

DATA  :  Z  =  J JJA222.  1000.  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  :  NONE 
BACKCASTING  :  OFF 

********■**■#■-*■*■***■***■*■**■***■*•»•■*■***■**  -tut*** *■***■**■»•■*  ♦-»■**•»■■*■***-*  ■*■***■»•»■ -*•* * 
UNIVARIATE  MODEL  PARAMETERS 

***»*********HHf*****«***#*********»+4***«*1HHHHH»*********+ **■**•*  +•* 

FACTOR  LAG  COEFFICIENT  T  RATIO 

**-*•  +**********************♦*+*******■* -*#***-*  **■#•*•»•***■*♦♦■**-**•*  •*■**-#■*  * 


1 

MEAN 

-.  5B224E+01 

2 

AUTOREGRESSIVE 

1 

1 

. 58105E-01 

1.  86 

n 

w 

AUTOREGRESSIVE 

1 

3 

. 20355E+00 

6.  47 

4 

AUTOREGRESSIVE 

2 

10 

. 87782E+00 

45.  61 

5 

MOVING  AVERAGE 

1 

2 

-.  13798E+00 

-4.  36 

S 

MOVING  AVERAGE 

o 

10 

. 29425E+00 

7.  61 

**t********+**********+*«**************4*****+*********+****t,**t 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


981 

987 


THE  RESIDUAL  STATISTICS 


. 49844E+00 
. 50810E-03 
. 82802E+00 


DEGREES  OF  FREEDOM  : 
NUMBER  OF  RESIDUALS  : 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODELS 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

Z<t)  =  Constant  +  Pi(l)  *  Z(t-l)  +  Pi(2)  *  Z(t-2)  ...  +  Pi(n)  *  Z  (t-n) 


!  THE  PI  WEIGHT 


LAG 


Constant 

Pi  ( 

1) 

Pi  < 

2) 

Pi  < 

3) 

Pi  < 

4) 

Pi  ( 

5) 

Pi  < 

10) 

Pi  < 

11) 

Pi  ( 

12) 

Pi  ( 

13) 

Pi  < 

14) 

Pi  ( 

15) 

Pi  < 

20) 

Pi  ( 

21) 

Pi  < 

22) 

Pi  < 

23) 

Pi  ( 

30) 

Pi  ( 

32) 

Pi  ( 

40) 

-. 13845E+02 
38105E-01 
. 13798E+00 
-.  19553E+00 
-. 1 9038E— 01 
.  26979E-01 
-. 1 1720E+01 
.  33837E-01 
. 1S171E+00 
. 11411E+00 
-. 22313E-01 
13745E-01 
-. 344S3E+00 
.  10019E-01 
.47578E-01 
. 33369E-01 
-. 10147E+00 
. 14000E-01 
29862E— 01 


MODELING  RESULTS  FOR  TIME  SERIES  JJJC222. 


**•»-»**■»•*•****-»*•»■»•».*******♦♦  *-**■*+*+ «Ht-*******-»-»  ******•**♦■***■*•»■•* 

DATA  :  2  =  JJJC222.  1000.  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  s  NONE 
BACKCASTING  s  OFF 

**•»•*•»■**■»■*■*•»•»***•»■»***•»•»■»*•»***■»*■*■»•».***•*■******♦♦■»-**♦♦**+*■♦*****  ***** 
UNIVARIATE  MODEL  PARAMETERS 

************************«*HMt******^****«-»-»*-» *-»*****•»•****  ***♦*■■*  ** 

FACTOR  LAG  COEFFICIENT  T  RATIO 


*-*■***■ ♦**•»*******•*•»•»■#**•»#**•»•***#**#■»****♦■»*  ***•»**•**•»■»**  *********** 


1 

MEAN 

.  731S3E+00 

2 

AUTOREGRESSIVE 

1 

1 

-.  97088E— 01 

-3.  05 

W 

AUTOREGRESSIVE 

1 

2 

.  18436E+00 

5.  87 

4 

AUTOREGRESSIVE 

1 

3 

. 1 76OOE+O0 

5.  53 

5 

AUTOREGRESSIVE 

2 

10 

.  21626E+00 

6.  75 

6 

AUTOREGRESSIVE 

2 

30 

. 15128E+00 

4.  74 

7 

MOVING  AVERAGE 

l 

20 

13290E+0O 

-4.62 

♦♦St-*****************************************-* -t*.**-***-*****-*  •*■*•#•****■# 


THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


111 33E+02 
11597E-01 
16219E+00 


DEGREES  OF  FREEDOM 
NUMBER  OF  RESIDUALS 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED.  AS  A  WEIGHTED  SUM  OF ' THE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODE-S 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

2  ( t )  =  Constant  +  Pi<l)  *  Z<t-1>  +  Pi  (2)  *  Z  Ct-2)  ...  +  Pi(n)  *  2  <t-n) 


LAG  1 

THE  PI  WEIGHT 

Constant 

. 97376E+00  . 

Pi  < 

1) 

. 97088E— 01 

Pi  ( 

2) 

-. 1S436E+00 

Pi  < 

3) 

-.  17600E+00 

Pi  ( 

10) 

—.21 626E+00 

Pi  ( 

11) 

20996E— 01 

Pi  ( 

12) 

. 39869E— 01 

Pi  ( 

13) 

. 38060E— 01 

Pi  ( 

20) 

. 15290E+00 

Pi  ( 

21) 

-.  14845E-01 

Pi  < 

22) 

. 28188E— 01 

Pi  ( 

23) 

.  26910E-01 

Pi  ( 

30) 

1 1821E+00 

Pi  < 

31) 

11477E-01 

Pi  ( 

32) 

. 21794E— 01 

Pi  < 

33) 

.  20805E— 01 

Pi  < 

40) 

— . 2337 8E— 01 

Pi  ( 

50) 

. 18075E-01 

MODELING  RESULTS 


TIME  SERIES  JJJA112. 


*■*•*•*■**•*■**■* *■*•* *************-*******'*****-********-*-*** **■**-****•****  *** 
DATA  :  Z  =  JJJA112.  1000.  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  s  NONE 

BACKCASTING  :  OFF 

#****************************************^  •**-*■***■**■*■**■*■■***■***  •**■**♦ 
UNIVARIATE  MODEL  PARAMETERS 

FACTOR  LAG  COEFFICIENT  T  RATIO 

**■***•»*#•**■****  •»**•*•******•»  •»*#■»*■»***•»■»******■»■»***•»  *♦*■»**•»*■»■******** 

1  MEAN  -.  67906E+01 

2  AUTOREGRESSIVE  1  1  -. 26564E+00  -8.68 

3  AUTOREGRESSIVE  2  10  .  52048E+00  .  17.31 

4  MOVING  AVERAGE  1  20  -.  1 721  1E+00  -4.  96 

*********#*****•********-****************♦*****************♦•**■*■**•** 


SUM  OF  SQUARES 
MEAN  SQUARE 
R  SQUARED 


THE  RESIDUAL  STATISTICS 


985 

989 


. 19409E+01 
. 19705E-02 
. 36941E+00 


DEGREES  OF  FREEDOM  : 
NUMBER  OF  RESIDUALS  : 


THE  PI  WEIGHTS 


UNIVARIATE  BQX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODELS. 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

Z(t)  *  Constant  ♦  Pi(l)  *  Z  (t-l>  +  Pi(2>  *  Z(t-£>  ...  +  Pi(n)  *  Z(t-n) 


LAG 

"  "the  pi  weight 

Constant 

-7ao709i+or 

!  Pi  < 

1 ) 

.  26564E+00 

:  Pi  < 

10) 

-. 52048E+00 

:  Pi  < 

11) 

-.  I3826E+00 

Pi  < 

20) 

. 17211E+O0 

:  Pi  < 

21) 

-. 45719E-01 

Pi  < 

30) 

. 89580E-01 

:  Pi  ( 

31) 

. 23796E-01 

Pi  ( 

40) 

29622E-01 

;  Pi  ( 

50) 

-.  15418E— 01 

. 

.  —  m  .  r 

« 


I  1 


MODELING  RESULTS  FOR  TIME  SERIES  JJJC112. 


*4*#**»»*****«*****#«4«*****4»*«******«*#***«***«*********««***** 

DATA  2  =  JJJC112.  1000.  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  s  NONE 
BACKCASTING  :  OFF 

***************************************************************** 
UNIVARIATE  MODEL  PARAMETERS 

***************************************************************** 
FACTOR  LAS  COEFFICIENT  T  RATIO 

***************************************************************** 


1 

MEAN 

.  631 13E+01 

2 

AUTOREGRESSIVE 

1 

1 

.  79309E—01 

2.  48 

3 

AUTOREGRESSIVE 

1 

2 

.  21900E+00 

7.  04 

4 

AUTOREGRESSIVE 

1 

3 

. 16608E+00 

5.  19 

5 

AUTOREGRESSIVE 

2 

10 

.  26382E+00 

8.  15 

6 

AUTOREGRESSIVE 

2 

20 

. 20279E+00 

6.  29 

7 

MOVING  AVERAGE 

1 

30 

83247E-01 

-2.  44 

***************************************************************** 

THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES  s  . 74O56E+01  DEGREES  OF  FREEDOM  :  370 
MEAN  SQUARE  :  .  76347E-02  NUMBER  OF  RESIDUALS  :  977 
R  SQUARED  s  . 22761 E+00 


THE  PI  WEIGHTS  . 

saasasaassaaai; 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODELS . 


2  (t) 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

*  Constant  +  Pi(l)  *  2<t-l)  *  Pi(2)  *  Z  <t-2>  ...  +  Pi  <n>  *  Z  (t-n> 


LAG  : 

THE  PI  WEIGHT 

Constant  i 

. 9saaiE+oi 

Pi  ( 

1)  : 

79309E-01 

Pi  < 

2) 

21900E+00 

Pi  < 

3)  ! 

16608E+00 

Pi  ( 

10)  : 

2S382E+00 

Pi  ( 

11)  : 

.  20S23E-Q1 

Pi  ( 

12)  : 

. 37773E-01 

Pi  ( 

13)  : 

. 43816E-01 

Pi  < 

20)  ! 

20279E+00 

Pi  ( 

21)  : 

.  160Q3E— 01 

Pi  < 

£2)  i 

. 44410E-01 

Pi  < 

23)  i 

. 33SS0E-01 

Pi  ( 

30) 

. 83247E-01 

Pi  < 

32)  : 

. 1823 IE— 01 

Pi  < 

33)  : 

. 1 3826E— 01 

Pi  < 

40)  : 

.  21962E-01 

Pi  ( 

so)  : 

. laaaiE-oi 

MODELING  RESULTS  FOR  TIME  SERIES  JJJA122. 


********************************************************  *♦*#****♦ 

DATA  :  Z  -  JJJA122.  lOOO.  '  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  j  NONE 
BACKCASTING  :  OFF 

**  ******************  ***************** ****** ******** *********  ***** 

UNIVARIATE  MODEL  PARAMETERS 

***************************************************************** 
FACTOR  LAG  COEFFICIENT  T  RATIO 

***************************************************************** 


1 

o 

MEAN 

AUTOREGRESSIVE 

1 

10 

S8009E+01 
.  S4461E+00 

26.  11 

3 

MOVING 

AVERAGE 

1 

1 

.  26997E+00 

8.  78 

4 

MOV  I NG 

AVERAGE 

2 

30 

.  1 1061E+00 

3.  34 

************************  »********♦*»♦*♦  ************************** 

THE  RESIDUAL  STATISTICS 


SUM  OF  SQUARES  : 
MEAN  SQUARE 


.  13O79E+01 
. 1 8336E-02 


DEGREES  OF  FREEDOM 
NUMBER  OF  RESIDUALS 


THE  PI  WEIGHTS 


UNIVARIATE  BQX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  TnE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODELS 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

2(t)  *  Constant  Pi(l)  *  Z(t-l)  +  Pi(2)  *  Z(t-2)  ...  +  Pi(n)  *  Z  (t-n) 


LAG 

THE  PI  WEIGHT! 

Constant 

17226E+02  ! 

Pi  (  1) 

26997E+00  i 

Pi  (  2) 

-. 72886E-01  ! 

Pi  (  3) 

19677E-01  ! 

Pi <  10) 

64461 E+00  ! 

Pi <  11) 

17403E+00  ! 

Pi (  12) 

46983E-01 

Pi <  13) 

12664E-01  ! 

Pi (  50) 

11061E+00  : 

Pi (  51) 

-.  29S61E-01  ! 

Pi  <  60) 

7129QE— 01  ! 

Pi <  61) 

19248E-01  ! 

Pi (100) 

-. 12234E-01  ! 

MODELING  RESULTS  FOR  TIME  SERIES  JJJC122. 


*■***-*■»**•*■»•*•»**■*****■*•********♦***■»*•»***»****•»******•»**•***♦*♦**•»*** 

DATA  :  Z  -  JJJC122.  1000.  1000  OBSERVATIONS 

DIFFERENCING  ON  Z  :  NONE 
BACKCASTING  :  OFF 
UNIVARIATE  MODEL  PARAMETERS 

FACTOR  LAG  COEFFICIENT  T  RATIO 


»*»-»***'»-»**-»*'»»**********'»-»***#***-»*'»-»-»****«****-»*«"»-»-»-»#-K-»*-»-»-»*-N-» 


1 

MEAN 

.  66319E+01 

2 

AUTOREGRESSIVE 

1 

1 

-.  12195E+00 

-3.  aa 

w 

AUTOREGRESSIVE 

1 

2 

.  15636E+00 

4.  99 

4 

autoregressive 

1 

3 

.  l89iaE+00 

6.  01 

5 

AUTOREGRESSIVE 

2 

10 

.  13376E+00 

4.  20 

6 

AUTOREGRESSIVE 

2 

20 

.  15959E+00 

5.  00 

7 

MOVING  AVERAGE 

1 

30 

-.  12349E+00 

-3.  77 

**■»•**■***#*•»-»****•»**■***■»#•»■»*****■»**♦**♦*■»*****•*****  +  ■******-»***  +  *** 


THE  RESIDUAL  STATISTICS 


SUM  OP  SQUARES 
MEAN  SQUARE 
R  SQUARED 


. 14183E+02 
. 14622E-01 
. 12S77E+00 


DEGREES  OF  FREEDOM  :  970 

NUMBER  OF  RESIDUALS  :  977 


THE  PI  WEIGHTS 


UNIVARIATE  BOX-JENKINS  MODELS  CAN  BE  EXPRESSED  AS  A  WEIGHTED  SUM  OF  THE  PAST 
PLUS  A  RANDOM  SHOCK.  DEFINED  AS  'THE  PI  WEIGHTS,  MODELS  EXPRESSED  IN  THESE 
TERMS  ARE  USEFUL  FOR  FORECASTING  AND/OR  COMPARISON  TO  OTHER  TYPES  OF  MODE-S. 


THE  GENERAL  FORM  OF  THE  MODEL  EXPRESSED  BY  THE  PI  WEIGHTS  IS  : 

Z(t)  *  Constant  +  Pt(l)  *  Z(t-l)  ♦  Pi(2>  *  Z(t-2)  ...  +  Pi  (n)  *  Z(t-n) 


LAG 

!  THE  PI  WEIGHT 

Constant 

. 85672E+01 

Pi  < 

1) 

.  121 95E+00 

Pi  < 

2) 

“  •  1 56 3 Sc  •♦‘00 

Pi  < 

3) 

1S918E+00 

Pi  < 

10) 

13376E+00 

Pi  ( 

11) 

-.  16312E-01 

Pi  < 

12) 

. 20915E— 01 

Pi  < 

13) 

. 25304E-01 

Pi  < 

20) 

1 5959E+00 

Pi  ( 

21) 

1 9462E-0 1 

Pi  ( 

22) 

.24S53E-01 

Pi  ( 

23) 

. 30190E-01 

Pi  ( 

30) 

. 12349E+00 

Pi  ( 

31) 

15060E-01 

Pi  ( 

32) 

. 1S309E-01 

Pi  ( 

33) 

. 23362E-01 

Pi  ( 

40) 

. 1S518E-01 

Pi  < 

50) 

. 1S707E-01 

Pi  ( 

60) 

1 5250E-01 

